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EXECUTIVE  SUMMARY 

The  Accessible  Antenna  Package  (APACX)  is  a  versatile,  automated, 
combined- antenna- and- propagation  set  of  computer  subroutines  for  predicting 
electric  far- field  strength,  directive  gains,  power  gains,  and  transmission 
loss  for  16  types  of  linear  antennas  that  are  typically  used  at  frequencies 
between  150  kHz  and  500  MHz.  The  types  of  antennas  considered  are: 
horizontal  dipole,  vertical  monopole,  vertical  monopole  with  radial-wire 
ground  screen,  elevated  vertical  dipole,  inverted-L,  arbitrarily  tilted 
dipole,  sloping  long-wire,  terminated  sloping  V,  terminated  sloping  rhombic, 
terminated  horizontal  rhombic,  side- loaded  vertical  half-rhombic,  horizontal 
Yagi-Uda  array,  horizontally  polarized  log-periodic  dipole  array,  vertically 
polarized  log-periodic  dipole  array,  curtain  array,  and  sloping  double 
rhomboid . 

Expressions  for  fields  from  these  antennas  for  direct  and  ground- 
reflected  waves  are  available  in  the  literature  and  one  form  is  anployed  in 
the  ECAC  computerized  SKY WAVE  model.  These  expressions  do  not  include  radio 
surface  wave  terms  and,  therefore,  are  not  applicable  to  ground-wave 
analyses.  Expressions  for  the  surface  waves  were  developed,  combined  with  the 
existing  expressions,  and  the  combinations  included  in  APACK.  Previously 
available  ground-wave  propagation  models  implicitly  assume  that  the  antenna  is 
an  infinitesimal  point  source,  i.e.,  a  Hertzian  dipole,  whereas  APACK 
explicitly  accounts  for  the  antenna  structure. 

APACK  accounts  for  the  effects  of  lossy  ground  under  the  antennas  in 
three  far- field  regions:  planar  earth  within  the  radio  horizon  (for  low 
antennas) ,  spherical  earth  within  the  radio  horizon  (for  high  antennas),  and 
the  diffraction  region  beyond  the  radio  horizon.  Automated  switching  criteria 
are  used  to  determine  which  of  the  three  regions  is  appropriate  for  the 
analysis. 

The  current  distribution  for  a  resonant  antenna  is  assumed  to  be 
sinusoidal,  and  the  current  distribution  for  a  traveling-wave  antenna  is 
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assumed  to  be  exponential.  Comparisons  of  APACK  predictions  with  various 
types  of  data  are  presented  indicating  the  versatility  and  reasonableness  of 
APACK.  In  particular,  for  an  electrical  large,  sloping,  double  rhomboid, 
APACK  predicts  behavior  that  is  similar  to  that  predicted  by  the  Numerical 
Electromagnetic  Cbde  (NEC) .  However,  one  significant  advantage  of  APACK  is 
that  it  requires  only  1/1000th  of  the  computer  run  time  of  the  NEC  analysis. 
It  is  found  that  APACK  is  applicable  except  when  the  antenna  employs  resonant 
elements  of  length  very  close  to  integral  multiples  of  a  wavelength  since  the 
sinusoidal  current  distribution  assumption  is  not  appropriate  in  this  case. 
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PREFACE 


The  Electromagentic  Compatibility  Analysis  Center  (ECAC)  is  a  Department 
of  Defense  facility,  established  to  provide  advice  and  assistance  on 
electromagnetic  compatibility  matters  to  the  Secretary  of  Defense,  the  Joint 
Chiefs  of  Staff,  the  military  departments  and  other  DoD  components.  Hie 
center,  located  at  North  Severn,  Annapolis,  Maryland  21402,  is  under  the 
policy  control  of  the  Assistant  Secretary  of  Defense  for  Communication, 
Command,  Control,  and  Intelligence  and  the  Chairman,  Joint  Chiefs  of  Staff,  or 
their  designees,  who  jointly  provide  policy  guidance,  assign  projects,  and 
establish  priorities.  ECAC  functions  under  the  executive  direction  of  the 
Secretary  of  the  Air  Force  and  the  management  and  technical  direction  of  the 
Center  are  provided  by  military  and  civil  service  personnel.  The  technical 
support  function  is  provided  through  an  Air  Force-sponsored  contract  with  the 
IIT  Research  Institute  (IITRI). 

To  the  extent  possible,  all  abbreviations  and  symbols  used  in  this  report 
are  taken  from  American  National  Standard  ANSI  Y10.19  (1969)  "Letter  Symbols 
for  Units  Used  in  Science  and  Technology"  issued  by  the  American  National 
Standards  Institute,  Inc. 

Users  of  this  report  are  invited  to  submit  comments  that  would  be  useful 
in  revising  or  adding  to  this  material  to  the  Director,  ECAC,  North  Severn, 
Annapolis,  Maryland  21402,  Attention;  XM. 
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Section  1 


SECTION  1 
INTRODUCTION 


BACKGROUND 


The  EL ectromagnetic  Compatibility  Analysis  Oenter  (ECAC)  is  a  Department 
of  Defense  facility,  established  to  provide  advice,  assistance,  and  analyses 
on  electromagnetic  compatibility  (EMC)  matters.  Many  analyses  have  been  and 
are  being  performed  in  the  frequency  region  between  150  kHz  and  500  MHz  for 
both  desired  and  interfering  signals  that  propagate  by  ground  wave  as  well  as 
sky  wave. 

The  antenna-gain  models  currently  used  by  ECAC  for  ionospheric  sky-wave 
analyses  are  those  developed  by  the  Institute  for  Telecommunication  Sciences 
(ITS)  as  part  of  a  model  for  predicting  ionospheric  propagation  (H  FMUFES-4 ) . 1 
Because  the  antenna-gain  models  developed  by  ITS  account  for  the  contributions 
of  the  direct  and  grounds reflected  waves  of  linear  antennas  mounted  over  lossy 
ground  but  do  not  consider  the  contribution  of  the  surface  wave,  the  models 
are  not  suitable  for  use  in  ground-wave  analyses. 

R>r  sky-wave  calculations,  it  is  proper  to  neglect  the  contribution  of 
the  surface  wave.  Also,  it  is  true  that  the  contribution  of  the  surface  wave 
is  small  for  the  case  of  a  horizontally  polarized  wave.  However,  the 
contribution  of  the  surface  wave  is  significant  for  ground-wave  calculations 
involving  vertically  polarized  waves. 

The  surface  wave  is  guided  along  the  surface  of  the  earth,  much  as  an 
electromagnetic  wave  is  guided  by  a  transmission  line.  Since  the  surface  wave 


1Barghausen,  A.  F. ,  Finney,  J.  W. ,  Proctor,  L.L. ,  and  Schultz,  L. D. , 
Predicting  Long-Term  Operational  Parameters  of  High-Frequency  SKYVAVE 
Telecommunication  Systans,  ESSA  Technical  Report  ERL  1 1 0 -ITS-78 , 
Institute  for  Telecommunication  Science,  Boulder,  00,  May  1969. 
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is  affected  by  losses  in  the  earth,  its  attenuation  is  directly  affected  by 
the  permittivity  and  conductivity  of  the  earth. 

ttien  both  the  transmitting  and  receiving  antennas  are  located  at  the 
surface  of  the  earth,  the  direct  and  ground- reflected  waves  cancel  each 
other.  In  this  case,  the  transmitted  wave  reaches  the  receiving  antenna 
entirely  by  means  of  the  surface  wave,  assisting  that  there  is  no  sky-wave  or 
tropospheric-scatter  propagation. 

ECAC  has  an  operational  method- of- moments  computer  program,  the  Numerical 
Electromagnetic  Oode  (NEC),  that  predicts  the  gain  pattern  of  an  antenna  in 
the  presence  of  a  lossy  ground  and  includes  the  surface  wave  term.  NEC 
requires  as  input  detailed  description  of  the  antenna  dimensions  and  location 
and  sizes  of  conducting  obstacles  (such  as  guy  wires) .  In  addition,  NEC 
requires  considerable  computer  time  to  make  its  predictions. 

Other  models  in  use  at  ECAC3'4  predict  basic  transmission  loss  without 
calculating  field  strength  for  a  Hertzian  dipole  but  do  not  account  for  the 
actual  antenna  configuration.  Therefore,  an  automated  model  was  needed  to 
calculate  the  gain  and  field  strengths  for  commonly  used  linear  antennas . 

This  model  should  be  capable  of  predicting  rapidly  directive  and  power  gains 
suitable  for  both  ground-wave  and  sky-wave  analyses  in  addition  to  predicting 
ground-wave  transmission  loss  in  the  150  kHz  to  500  MHz  region.  The  model  was 


2 

Burke,  G.  and  Foggio,  A.,  ttamerical  Electromagnetic  Cbde  (NEC)  —  Method  of 
Moments,  Part  I:  Program  Description  -  Theory,  Ifert  II:  Program  Description 
Code,  and  Part  III:  User's  Guide,  Technical  Docunent  116,  Naval  Ocean 
Systems  Center,  San  Diego,  CA,  18  July  1977  (revised  2  January  1980). 

3Meidenbauer ,  R. ,  Qiang,  S. ,  and  Duncan,  M. ,  A  Status  Report  on  the 
Integrated  Propagation  System  (IPS),  ECAC-TO-78-023,  Electromagnetic 
Compatibility  Analysis  Center,  Annapolis,  MD,  October  1978. 

“Slaiuzzo,  M.A.  and  Brazier,  W. E. ,  A  Theoretical  Ground  VAve  Propagation 
Model  -  NX  Model,  ESD-TR-68-315,  Electromagnetic  Compatibility  Analysis 
Center,  Annapolis,  MD,  December  1968. 
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designed  to  be  modular  in  the  sense  that  it  would  be  a  collection  of 
subroutines  to  be  called  by  the  various  other  ECAC  models  as  needed.  This 
collection  of  subroutines  is  termed  APACK. 

OBJECTIVES 


Hie  objectives  of  this  report  are: 

1.  Ib  docvment  the  equations  that  comprise  APACK  and  provide  an  overview 
of  the  methods  used  to  develop  the  equations 

2.  lb  compare  APACK  predictions  for  antenna  gain  and  transmission  loss 
with  other  available  data . 

APPROACH 


The  Hertzian  Dipole 

The  Hertzian  dipole  is  used  as  a  building  block  in  the  analysis  of  linear 
antennas,  because  the  fields  of  linear  antennas  are  given  in  terms  of 
integrals  of  the  current  distribution.  Since  the  current  distribution  of  the 
Hertzian  dipole  is  assumed  to  be  uniform,  the  integrals  describing  the  fields 
of  a  Hertzian  dipole  are  simplified.  Actual  linear  antennas  with  nonuniform 
current  distributions  can  then  be  simplified  for  the  purposes  of  analysis  by 
considering  them  to  be  comprised  of  a  superposition  of  many  Hertzian  dipoles, 
each  with  a  uniform  current  distribution. 

The  equations  of  a  Hertzian  dipole  over  lossy  planar  earth  are  well 
Known. 6  Hue  actual  antenna  structure  was  considered  as  a  superposition  of 
Hertzian  dipoles  giving  a  sinusoidal  current  distribution  for  the  antenna. 


5Banos,  A.,  Jr.,  Dipole  Radiation  in  the  Presence  of  a  Conducting  Half-Space, 
Pergamon  Press,  Oxford,  England,  1966. 

®Weeks,  W.L. ,  Antenna  Engineering,  McGraw-Hill,  New  York,  NY,  1968. 
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Appropriate  formulations  were  applied  to  extend  the  equations  for  the  fields 
over  planar  earth  to  spherical  earth  within  the  radio  horizon  and  to  the 
diffraction  region  beyond  the  radio  horizon. 

Even  the  analysis  of  the  fields  of  a  Hertzian  dipole  is  not  simple  when 
it  is  located  over  lossy  planar  earth.'0'3  The  analysis  herein  is  based  on 
that  of  Norton.10  Norton  provides  practical  formulas  to  calculate  the  field 
intensity  of  electric  and  magnetic  dipoles  over  lossy  ground.  Norton  also 
separated  his  solutions  into  a  space  wave  and  a  surface  wave.  It  has  been 
shown  that  Norton's  formulation  gives  close  agreement  with  the  numerical 
solution  of  Sommerfeld’ s  equations.11 

Linear  Antennas  Over  Lossy  Planar  Earth 

APACK  predicts  basic  transmission  loss  from  calculated  values  of  electric 
field  strength.  The  APACK  field  strengths  account  for  the  actual  antenna 
structure  and  include  the  effects  of  the  contributions  of  the  direct,  ground- 
reflected,  and  surface  waves. 


7 Sommerfeld,  A. ,  "Uber  die  Ausbreitung  der  Wellen  in  der  drahtlosen 
Telegraphic,"  Ann.  Physik,  Vol.  28,  1909,  pp.  665-736. 


O 

Sommerfeld,  A.,  "Uber  die  Ausbreitung  der  Wellen  in  der  drahtlosen 
Telegraphic,"  Ann.  Physik,  Vol.  81,  1926,  pp.  1135-1153. 


g 

Sommerfeld,  A.,  Partial  Differential  Equations,  Academic  Press,  New  York, 
NY,  1949. 


10Norton,  K. A.,  "The  Propagation  of  Radio  Waves  Over  the  Surface  of  the 

Earth  and  in  the  Upper  Atmosphere":  Part  I,  Proc.  IRE,  Vol.  24,  October  1936, 
pp.  1369-1389;  Part  II,  Proc .  IRE ,  Vol.  25,  September  1937,  pp.  1203-1236. 

^Kuebler,  W.  and  Snyder,  S. ,  The  Sommerfeld  Integral,  Its  Computer 
Evaluation  and  Application  to  Near  Field  Problems,  ECAC-TN-75-002, 
Electromagnetic  Compatibility  Analysis  Center,  Annapolis,  MD, 

February  1975. 
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The  equations  developed  for  the  contributions  of  the  direct  and  ground- 

reflected  waves  from  the  16  linear  antennas  listed  in  TABLE  1  have  been 

derived  previously  by  others.  Laitenen  derived  formulas  for  the  fields  of 

linear  antennas  assvxning  a  sinusoidal  current  distribution.  He  included 

only  the  direct  and  ground-reflected  waves  in  his  calculation.  and 

1 3 

Walters  derived  similar  formulas,  and  Ma  later  extended  the  results  using  a 
more  accurate  three-term  current  distribution.14  However,  Laitenen,  Ma,  and 
Valters  (see  References  12,  13,  and  14)  do  not  consider  the  surface  wave. 

The  surface  wave  formulas  used  in  APACK  were  formulated  in  terms  of  the 
fields  of  a  current  element  over  lossy  planar  earth  employing  Morton's 
equations.  The  expressions  for  directive  gain  and  power  gain  were  derived 
from  the  expressions  for  the  far  fields. 

Thus,  the  initial  steps  in  the  development  of  the  APACK  equations  were 
the  formulation  of  the  expressions  for  directive  and  power  gains  in  terms  of 
the  far  fields  and  the  formulation  of  the  expressions  for  the  far  fields  of  a 
current  element  over  lossy  planar  earth.  These  formulations  are  presented  in 
Section  2. 

Assigned  Antenna  Current  Distribution 


The  determination  of  the  fields  of  the  current  element  requires  a 
knowledge  of  the  current  distribution  on  the  element  which,  in  turn,  requires 
a  knowledge  of  the  current  distribution  on  the  actual  antenna  structure.  The 


1  2 

Laitenen,  P. ,  Linear  Cbmmunication  Antennas,  Technical  Report  No.  7,  u.  S. 

Army  Signal  Radio  Propagation  Agency,  Rort  Monmouth,  NJ,  1959. 

'^Ma,  M.T.  and  V&lters,  L.  C. ,  Power  Gains  for  Antennas  Over  lossy  Plane  Ground, 
Technical  Report  ERL  104-ITS  74,  Institute  for  Telecommunication  Sciences, 
Boulder,  00,  1969. 

14Ma ,  M.T. ,  Theory  and  Application  of  Antenna  Arrays,  Wiley  Interscience, 

New  York,  NY,  1974. 
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equations  derived  for  APACK  assume  that  the  current  distribution  for  a 
resonant  antenna  is  sinusoidal  and  that  the  current  distribution  for  a 
traveling-wave  antenna  is  exponential . 

TABLE  1 

TYPES  OF  LINEAR  ANTENNAS  PRESENTLY  INCLUDED  IN  APACK 


1.  Horizontal  dipole 

2.  Vertical  monopole 

3.  Vertical  raonopole  with  radial-wire  ground  screen 

4.  Elevated  vertical  dipole 

5.  Inverted-L 

6.  Arbitrarily  tilted  dipole 

7.  Sloping  long-wire 

8.  Terminated  sloping-V 

9.  Terminated  sloping  rhombic 

10.  Terminated  horizontal  rhombic 

11.  Side- loaded  vertical  half-rhombic 

12.  Horizontal  Yagi-Uda  array 

13.  Horizontally  polarized  log-periodic  dipole  array 

14.  Vertically  polarized  log-periodic  dipole  array 

15.  Qirtain  array 

16.  Sloping  double  rhomboid 

Sinusoidal  current  distribution  was  first  treated  by  Pocklington. 1 5 
However,  it  is  well  known  that  sinusoidal  current  distribution  is  only  an 
approximation.  Schelkunoff  and  Ftiis^®  made  the  following  statements  on  the 
relationship  between  the  current  distribution  and  the  radiation  pattern  as 
well  as  the  radiated  power. 


1 5Pocklington,  H.E.,  "Electrical  Oscillations  in  Wires,"  Qjmb.  Phil.  Soc. 
Proc . ,  25  October  1897,  pp.  324-332. 

^Schelkunoff,  S.A.  and  Ftiis,  H.T. ,  Antennas,  Theory  and  Practice, 

John  Wiley  and  Sons,  New  York,  NY,  1952. 
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The  radiation  pattern  and  power  radiated  by  the  antenna  are 
insensitive  to  errors  in  the  assumed  form  of  current 
distribution.  The  antenna  current  must  be  known  more 
accurately  if  we  are  interested  in  the  minima  in  the  radiation 
pattern;  but  it  need  not  be  known  accurately  otherwise.  .  .  . 

A  primary  consideration  in  the  design  of  APACK  was  that  the  resulting 
model  execute  in  minimal  time  for  each  required  value  of  gain  or  transmission 
loss.  This  consideration  is  due  to  the  fact  that  hundreds  or  thousands  of 
values  of  gain  may  be  necessary  for  an  analysis  of  one  circuit  throughout  the 
HF  band  or  thousands  of  values  of  transmission  loss  may  be  required  for 
analyzing  a  ground-wave  circuit  over  a  wide  range  of  frequencies.  The  need 
for  large  numbers  of  predictions  involving  many  frequencies  precludes  the  use 
of  method-of-moments  models  that  use  matrix  techniques  to  compute  the  antenna 
current  distribution  at  each  frequency. 

Also,  because  it  was  important  that  the  model  execute  in  minimal  time,  a 
sinusoidal  current  distribution  was  assumed  for  resonant  antennas.  This 
distribution  is  not  as  accurate  or  elegant  as  the  three-term  current 
distribution  used  by  Ma  (see  Reference  14),  but  it  does  provide  reasonable 
results  except  when  the  lengths  of  the  resonant  elements  are  very  close  to 
integral  multiples  of  a  wavelength. 

The  restrictions  associated  with  the  sinusoidal  current  distribution  for 
resonant  antennas  do  not  arise  with  traveling-wave  (i.e.,  nonresonant) 
antennas  when  an  exponential  current  distribution  is  assumed.  The  use  of  the 
exponential  current  distribution  was,  therefore,  believed  to  be  reasonable 
without  restrictions. 

Bctension  of  the  Rjrmulation  to  Spherical  Earth 

then  the  antenna  is  located  close  to  the  surface  of  the  earth,  the  earth 
can  be  considered  as  planar.  However,  when  the  feed  point  of  the  antenna  is 
located  several  wavelengths  or  more  above  the  surface,  the  earth  can  no  longer 
be  considered  as  planar  for  calculations  within  the  radio  horizon. 
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Therefore,  the  formulation  of  the  fields  of  a  current  element  over  lossy 
planar  earth  (presented  in  Section  2)  was  extended  to  account  for  the 
curvature  of  the  earth.  Although  the  Bremmer  formulation17  can  be  used  within 
the  radio  horizon,  many  terms  of  the  series  for  the  fields  are  required  to 
make  the  series  converge  in  this  region. 

Norton  provided  approximate  formulas  that  account  for  the  curvature 
without  resorting  to  the  rigorous  Bremmer  techniques  and  are  computationally 
efficient.  Norton's  formulas  were  thus  used  to  extend  the  APACK  equations  to 
account  for  the  curvature  of  the  earth  within  the  radio  horizon.  This 
extension  is  presented  in  Section  3. 

Extension  of  the  Fbrmulation  to  the  Diffraction  Region 

In  the  diffraction  region  beyond  the  radio  horizon,  the  basic  formulation 
presented  in  Section  2  for  the  fields  of  a  current  element  over  lossy  planar 
earth  was  modified  by  making  use  of  the  radiation  vector  and  the  Bremmer 
formulation  (see  Reference  17).  The  radiation  vector  is  similar  to  the 
formulation  presented  in  Section  2  and  accounts  for  the  geometry  of  the 
antenna  structure.  The  Bremmer  secondary  factor  accounts  for  the  geometry  of 
the  path. 

The  resulting  electric  far-field  components,  presented  in  Section  4,  are 
in  terms  of  the  product  of  the  radiation  vector  and  the  Bremmer  secondary 
factor.  The  advantage  of  this  formulation  for  APACK  calculations  in  the 
diffraction  region  was  that  the  radiation  vector  and  Bremmer  secondary  factor 
could  be  calculated  independently,  and  thus  the  two  routines  could  be  modular. 


17Brenuner,  H.,  Terrestrial  Radio  Vfeves,  Elsevier  Publishing  00., 

New  York,  NY,  1949. 

18No  rton,  K.A.,  "The  Calculation  of  Ground  wave  Field  Intensity  Over  a 
Finitely  Conducting  Spherical  Earth,"  Proc.  IRE,  Vol.  29,  No.  12, 
December  1941,  pp.  623-639. 
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Resulting  Formulas  for  Field  Strength  and  Gain 

The  formulation  presented  in  Section  2  with  the  extensions  presented  in 
Sections  3  and  4  was  used  to  derive  equations  for  the  electric  far-field 
components,  directive  gain,  and  power  gain  for  the  16  types  of  linear  antennas 
considered  (see  TABLE  1).  Sinusoidal  current  distribution  was  assumed  for 
resonant  antennas,  and  exponential  current  distribution  was  assumed  for 
travel ing-wave  (i.e.,  nonresonant)  antennas. 

lb  simplify  the  use  of  this  report  for  reference  purposes,  key  equations 
for  determining  the  field  strengths  and  directive  gains  of  the  antennas  are 
listed  in  TABLE  2  in  Section  5.  TABLE  2  refers  to  appropriate  equations  in 
APPENDIXES  A  and  D  in  which  the  mathematical  details  for  each  antenna  are 
provided.  In  addition  to  TABLE  2,  Section  5  includes  a  figure  showing  the 
geometry  of  each  antenna  and  a  brief  introduction  to  each  antenna  for  the 
uninitiated  reader. 

Transmission  loss 


Because  the  basic  equations  derived  for  APACK  are  for  the  electric  far 
fields  of  the  antennas,  calculations  of  ground-wave  transmission  loss  are  also 
straightforward.  The  general  equations  presented  in  Section  6  were  used  to 
derive  the  transmission  loss  relative  to  free-space  loss  in  terms  of  the  power 
gains  of  the  antennas  and  the  ratio  of  the  actual  disturbed  field  at  the 
observation  point  to  the  free-space  field  at  the  observation  point . 

Automated  criteria,  also  presented  in  Section  6,  were  used  to  determine 
in  which  of  the  three  regions  the  far-field  observation  point  lies:  planar 
earth  (for  low  antennas),  spherical  earth  (for  high  antennas),  or  the 
diffraction  region  beyond  the  radio  horizon.  These  criteria  are  based  simply 
on  the  path  length,  operating  frequency,  and  heights  of  the  transmitting  and 
receiving  antenna  feed  points  above  ground.  The  criteria  were  also  used  in 
conjunction  with  the  equations  listed  in  TABLE  2  of  Section  5  so  that 
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calculations  of  field  strength  and  gain,  in  addition  to  transmission  loss, 
would  automatically  account  for  the  location  of  the  observation  point. 

Comparisons  Between  APACK  Predictions  and  Other  Data 

Gains  predicted  by  APACK  were  compared  with  gains  obtained  from  other 
sources.  These  comparisons,  presented  in  Section  7,  include  an  example  of 
each  of  the  16  types  of  antennas  and  various  ground  constants  (i.e., 
conductivities  and  permittivities) .  fehile  the  comparisons  are  not  exhaustive, 
they  do  indicate  that  APACK  gain  predictions  are  reasonable. 

Transmission- loss  predictions  made  by  APACK  were  compared  with 
transmission-loss  predictions  obtained  from  other  sources  and  are  presented  in 
Section  8.  Various  ground  constants,  typical  of  soil  and  sea  water,  were  used 
to  demonstrate  the  versatility  of  the  model . 
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SECTION  2 

GENERAL  CONSIDERATIONS  FOR  CALCULATING  GAIN  AND 
FIELD  INTENSITIES  OF  A  LINEAR  ANTENNA 


INTRODUCTORY  REMARKS 


The  directive  and  power  gains  calculated  by  APACK  use  the  definitions  of 
gains  in  terms  of  the  electric  far  fields  of  the  antenna  being  considered. 
Electric  far  field  intensities  are  calculated  from  those  of  a  current  element 
located  above  planar  earth.  Fresnel  reflection  coefficients  are  used  to 
account  for  the  presence  of  ground  both  in  formulating  the  electric  far  fields 
and  in  formulating  the  radiation  resistance. 

This  section  presents  general  expressions  for  calculating  directive  and 
power  gain.  It  also  presents  the  expressions  for  the  electric  far  fields  of  a 
current  element  located  above  planar  earth.  The  formulation  for  the  current 
element  is  extended  in  Section  3  to  include  the  effects  of  spherical  earth 
within  the  radio  horizon  by  using  the  divergence  factor  with  the  Fresnel 
reflection  coefficients.  The  radiation  vector  and  Bremmer  secondary  factor 
are  used  in  Section  4  to  include  the  diffraction  region  beyond  the  radio 
horizon . 

GAIN  CALCULATIONS 


The  directive  gain  ( gd)  of  an  antenna  in  a  given  direction  is  defined  as 
the  ratio  of  the  radiation  intensity  in  that  direction  to  the  average  power 
radiated  per  unit  solid  angle.  Thus: 
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9,5(9/$)  =  directive  gain  in  the  direction  specified  by  the  spherical- 


$(9,<f>)  = 


coordinate  angles  0  and  $  (numerical  ratio) 

radiation  intensity  in  the  direction  specified  by  the  spherical- 
coordinate  angles  9  and  $,  in  watts/steradian 


$av  =  average  power  radiated  per  unit  solid  angle,  in  watts/steradian 


Wr  =  total  power  radiated  by  the  antenna ,  in  watts , 


The  radiation  intensity  in  a  given  direction  can  also  be  defined  in  terms 
of  electric  field  intensity  by  : 


♦<ef$) 


where 


r2 1  E  ( 9  ,4>)  I  2 


(2-2) 


r  =  spherical- coordinate  radial  distance,  in  meters 

E(9,<f>)  =  electric  far  field  in  the  direction  specified  by  the  spherical- 
coordinate  angles  9  and  $,  in  volts/meter. 

From  Equations  2-1  and  2-2: 
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where 

=  current  at  the  antenna  feed  point  (base)  ,  in  amperes 

Ryk  =  radiation  resistance  referred  to  the  antenna  feed  point  (base) , 
in  ohms. 

Hie  maximum  value  of  the  directive  gain  is  called  the  directivity.  Hie 
directivity  is  sometimes  loosely  referred  to  as  the  "gain,"  but  this  usage  is 
depracated. 

Hie  power  gain  ( gp )  of  an  antenna  in  a  given  direction  is  defined  by: 


gp(9,<t>) 


4 it  $(9  ,4>) 

w. 

in 


(2-4) 


where 


gp (£,4>)  =  power  gain  in  the  direction  specified  by  the  spherical-coordinate 
angles  9  and  4>  (numerical  ratio) 


W^n  =  power  input  to  the  antenna,  in  watts. 
Hie  radiation  efficiency  (n)  is  defined  by: 


W 


in 


( n  is  a  miner ical  ratio,  0  <  n  <  1.) 
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From  Equations  2-1  and  2-5: 


0  g.(9,<J>) 
a 


4it  <&  ( 6  ,  4> ) 


g  (8,*) 


(2-6) 


Since  n  _<  1,  gp  _<  g^,  and  the  difference  between  gp  and  g^  can  be 
significant.  The  radiation  efficiency  can  also  be  expressed  as: 


R  .  +  R. 
rb  loss 


(2-7) 


where 


Rloss  =  losses  associated  with  the  antenna,  in  ohms. 


From  Equations  2-3,  2-6,  and  2-7: 


gp  (0,4>) 


r2  |  i(Q  ,4.)  |  2 

30  I  2  (R  .  +  R.  ) 
b  rb  loss 


(2-8) 


In  spherical  coordinates  r ,  9  ,  <J) 


E(8,<j>r 


|Ee*2  +  I*/ 


(2-9) 


The  power  gain  in  decibels  (Gp)  is  given  by: 


0^  (9,<t>)  =  10  log1Q  gp  ( 9  ,<)> ) 


(2-10) 
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FIELD -INTENSITY  CALCULATIONS 


Consider  a  linear  antenna  element  as  shown  in  Figure  1  with  the  XY-plane 
being  lossy  planar  earth.  The  electric- field  components  produced  by  the 
antenna  at  a  point  P  (r,8,$),  including  the  contributions  of  the  direct, 
ground-reflected,  and  surface  waves  are  given  in  spherical  coordinates  ( r , 0  , <}> ) 

|cos  a'  cos  (<j>  -  )  cos  0 


by : 


-jkr 


=  j3  0k 


-jskH  cos6 

l  jks  cost  s 

x  /  I(s)  e  (1  -  R  e  (2-11) 

0  v 


+  surface  wave  terms)  ds  -  sin  a'  sin0 

jks  cost 


-j2kH  cos0 
s 


x  /  I(s)  e 
0 


(1  +  R  e 
v 


+  surface-wave  terms)  ds I 
-jkr 

e  i  jks  cost 

=  -j30k  -  cos  a'  sin  (t  -  t")  /  I(s)  e 

r  0 


-j2kH  cos0 
s 

x  (1  +  R  e  +  surface-wave  terms)  ds 

h 


(2-12) 


where 

0-,  t-  component  of  the  electric  far  field,  in  volts/meter 
2  ir/X 

wavelength ,  in  meters 


15 


Arbitrarily  oriented  current  element  over  lossy  planar  earth. 
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r 


l 

I(s) 

s 

* 


*v 


«h 


0' 


radial  distance  from  the  origin  to  the  far-field  point  P(r,6,$),  in 
meters 

angle  between  the  antenna  element  and  its  projection  on  the  XY-plane, 
in  degrees 

angle  between  the  X-axis  and  the  projection  of  the  antenna  element  on 

the  XY-plane,  in  degrees 

length  of  the  antenna  element,  in  meters 

linear  current  density  for  the  antenna  element,  in  amperes/meters 
linear  coordinate  coinciding  with  the  antenna  element,  in  meters 
angle  between  the  antenna  element  and  the  line  from  the  origin  to  the 
far-field  point  P(r,9,$),  in  degrees 

Fresnel  reflection  coefficient  for  the  vertically  polarized 
component,  defined  below 

height  of  the  current  element  at  ds  above  the  XY-plane,  in  meters 
differential  element  of  length  along  the  antenna  element,  in  meters 
FTesnel  reflection  coefficient  for  the  horizontally  polarized 
component,  defined  below 

angle  between  the  Z-axis  and  the  antenna  element,  in  degrees. 


Ttie  angle  ij>  can  be  calculated  from: 


cosi|>  =  cos0  cos0'  +  sin  9  sinO'  cos  ( 4>  —0 ) 


(2-13) 


where  the  primed  coordinates  refer  to  locations  on  the  antenna  element  and 
unprimed  coordinates  refer  to  the  far-field  observation  point  P(r,9,$). 

The  Etesnel  reflection  coefficients  for  the  vertically  and  horizontally 
polarized  components  of  the  electric  field  (Ry  and  R^,  respectively)  are  given 
by : 


n2  cos  eT  ~  \ 

2  2 
n  -  sin  0 

r 

n2  cos  ®r  +  ^ 

~7  '  '  2 — 

n  -  sin  0 

r 

17 
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and 


V2  2 

n  -  sin  9 

R  =  - - -  (2-15) 

cos  9  +  V  n  -  sin  9 

r  ’  r 

where 

n  ■  refractive  index  of  the  medium  under  the  antenna,  defined  below 

9r  =  angle  between  the  line  from  the  image  of  the  current  element  at  ds  to 
the  far- field  point  P(r,9,iji)  and  the  Z-axis,  in  degrees. 

The  refractive  index  of  the  medium  under  the  antenna  ( n)  is  : 


n 


2 


e 

r 


18000a 
f  MHz 


(2-16) 


where 

er  =  relative  dielectric  constant  of  the  medium  (dimensionless) 
a  =  conductivity  of  the  medium,  in  mhos/meter 
f =  frequency,  in  MHz. 

The  input  resistance  of  an  antenna  (Rin) »  including  the  effects  of  lossy 
ground,  is  calculated  from: 

I  Rin  -  *11  (2-17) 

I 
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where 


2m 

C 


input  resistance,  in  ohms 

real  part  of  the  antenna  self- impedance,  in  ohms 

mutual  impedance  between  the  antenna  and  its  image  in  perfectly 

conducting  ground,  in  ohms 

factor  to  account  for  lossy  ground,  as  given  below. 


c 


cos  a"  +  j  R  sin  a') 
v 


(2-18) 


R^'  is  obtained  by  evaluating  Equation  2-15  at  9r  =  0°  to  give: 


«h 


1-n 

1+n 


(2-19) 


Rv'  is  obtained  by  evaluating  Equation  2-14  at  6  =  0°  to  give: 


R  " 


v 


n-1 
n+ 1 


(2-20 ) 


VARIATION  OF  ANTENNA  GAIN  AS  A  FUNCTION  OF  EAR- FIELD  DISTANCE 


The  definitions  of  directive  gain  and  power  gain  of  an  antenna  are  the 
result  of  considering  the  properties  of  an  antenna  located  in  free  space. 

Mien  the  fields  vary  as  1/r,  the  expressions  for  directive  gain  and  power  gain 
(see  Equations  2-3  and  2-8)  are  not  a  function  of  the  distance  from  the 
antenna  to  the  observer. 

However,  when  an  antenna  is  deployed  over  lossy  ground,  the  far  fields  no 
longer  vary  exactly  as  1/r.  Thus,  the  presence  of  lossy  ground  makes  the  gain 
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i 


a  function  of  the  far-field  distance  between  an  observer  and  the  antenna  under 
consideration.  Wfeeks  presents  the  following  discussion  of  this  point  (see 
Reference  6,  pp.  346-347)  . 


The  measurements  and  application  of  the  concepts  of  power  gain, 
directivity,  and  radiation  efficiency  are  particularly  difficult 
when  the  antennas  must  operate  in  lossy  environments.  The 
definitions  of  these  quantities  were  conceived  initially  to  provide 
simplicity  in  free-space  environments;  they  do  not  lend  simplicity 
elsewhere.  Honest  and  meaningful  evaluation  are  also  clouded  by  the 
divergent  motives  and  objectives  of  the  pure  antenna  engineer  and 
the  user.  R>r,  if  system  performance  is  degraded  by  an  environment 
over  which  he  has  no  control,  the  antenna  engineer  understandably 
would  like  to  describe  the  gain  and  efficiency  of  his  product  as 
they  wauld  be  in  an  ideal  environment,  since  after  all  "his"  part  of 
the  system  is  "working  fine." 

Oi  the  other  hand,  the  systems  engineer  and  user  are  concerned  with 
overall  performance  and,  understandably,  are  not  kindly  disposed 
toward  an  antenna  that  would  "work  fine"  in  an  ideal  environment  but 
fails  to  provide  conmunication  in  the  actual  environment.  So  to 
sell  his  product,  the  antenna  engineer  must  evaluate  his  product  as 
it  would  function  in  an  actual  environment. 

The  most  common  application  in  which  these  difficulties  are  manifest 
is  that  of  antennas  for  operation  on  or  close  to  the  surface  of  the 
earth,  at  frequencies  below,  say  30  MHz.  Here,  one  of  the  most 
basic  difficulties  is  that  the  definitions  of  power  gain  and 
directive  gain  have  in  them  inherently  the  assumption  that  the 
fields  vary  as  1/r  and  that  the  power  density  varies  as  1/r.  In  the 
actual  earth  environment,  this  is  not  the  case;  in  the  direction 
near  the  ground,  the  field  falls  off  faster  than  this,  perhaps  much 
faster.  Thus,  with  a  direct  application  of  the  usual  definition, 
the  gain  depends  on  distance  from  the  antenna.  The  radiation 
pattern  may  also  depend  on  the  distance  from  the  antenna  structure, 
even  though  the  distance  may  be  large  compared  with  the  free-space 
distant- field  criteria.  At  large  distances,  even  for  vertical 
antennas,  there  is  usually  essentially  a  null  at  the  horizon. 


Since  the  electric  far  field  due  to  the  surface  wave  attenuates  more 
rapidly  than  1/r  the  gain  obviously  is  a  function  of  far-field  distance  when 
the  surface-wave  contribution  is  included  in  the  total  field.  Less  obvious  is 
the  fact  that  gain  can  be  a  function  of  fai^field  distance  even  if  the 
surface-wave  contribution  is  not  included.  This  fact  will  be  addressed  below. 


ESD-TR-80-1 02 


Section  2 


ltie  objective  of  the  SKYVAVE  antenna  gain  routines  (see  Reference  1)  is 
to  calculate  vertical  gain  patterns  for  use  in  ionospheric  propagation 
predictions,  so  the  SKYfcAVE  gain  equations  do  not  include  distance  at  all. 
APACX  calculates  both  sky-wave  and  ground-wave  field  strengths  and  gains . 

A PACK  directive  and  power  gain  are  calculated  from  Equations  2-3  and  2-8, 
respectively,  from  the  electric  far  fields. 

As  long  as  the  field  strength  varies  as  1/r,  gain  is  independent  of 
distance.  The  reason  that  gain  can  vary  as  a  function  of  far- field  distance 
even  if  the  surface-wave  contribution  is  not  included  is  that  the  Fresnel 
reflection  coefficients  depend  on  the  angle  0r  which  changes  with  distance. 
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SECTION  3 

EXTENSION  OF  THE  FORMULATION  TO  SPHERICAL  EARTH 
WITH  THE  RADIO  HORIZON 


INTRODUCTORY  REMARKS 


The  formulation  of  the  far  fields  of  a  current  element  presented  in 
Section  2  assumes  that  the  element  is  located  above  planar  earth.  Wien  the 
feed  point  of  the  antenna  is  located  several  wavelengths  or  more  above  the 
surface  of  the  earth,  the  earth  can  no  longer  be  considered  planar. 

Rjr  rigorous  calculations  of  electric  field  intensities  over  a  spherical 
earth,  the  Bremmer  formulation  (see  Reference  17)  should  be  used.  However, 
since  the  Qremmer  series  requires  a  large  number  of  terms  for  convergence 
within  the  radio  horizon,  this  formulation  is  used  by  APACK  only  in  the 
diffraction  region  beyond  the  radio  horizon.  Within  the  radio  horizon,  the 
planar  earth  formulation  can  be  modified  to  account  for  the  curvature  of  the 
earth  without  resorting  to  the  rigorous  Bremmer  techniques  (see  Reference  17). 

Strictly  speaking,  a  spherical  reflection  coefficient  should  be  used  to 
include  the  effect  of  earth  curvature.  APACK  uses  the  Hresnel  reflection 
coefficient  within  the  radio  horizon  because  the  difference  between  the 
Fresnel  coefficient  and  the  spherical  coefficient  is  negligible  except  near 
the  horizon.  The  Fresnel  coefficient  must  be  appropriately  modified,  however, 
by  the  divergence  factor. 

CALCULATION  OF  THE  DIVERGENCE  FACTOR 

Wie  divergence  factor  (adiv)  ,  a  geometrical  quantity  independent  of 
frequency,  is  a  measure  of  the  extra  divergence  acquired  by  a  beam  of  rays 
after  reflection  from  a  spherical  surface  as  opposed  to  a  planar  surface.  The 
divergence  factor  is  defined  by : 
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a 


div 


a 

e 


b  b' 


(  C+C*  )  J 
sin6  ( C  b' 


sin  8  cosB 


cosa'  +  C^b  cosa) 


when 


(3-1  ) 


a  =  6370 

e 


0.005577N 


1  -0 . 04665e 


-1 


(3-2) 


-0.1057h 

N  =  N  e  * 

s  o 


(3-3) 


where 

a0  =  effective  earth  radius,  in  kilometers 

Ng  =  surface  atmospheric  refractivity ,  in  N-units 

Nq  =  surface  atmospheric  refractivity  reduced  to  sea  level,  in  N-units 
hg  =  elevation  of  the  surface  above  mean  sea  level ,  in  kilometers • 

(If  Nq  and  hg  are  not  given,  Ng  =  301  is  assumed.)  The  quantities  C,  C' ,  a, 
a’ ,  6,  and  0  are  as  defined  in  Figure  2.  The  quantities  b,  b’ ,  and  d  are 
defined  by: 


a  + 
e 

h2 

(3-4) 

a  + 
e 

h1 

(3-5) 

a  0 
e 

(3-6) 
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where 

h1  =  height  of  the  transmitting  antenna  feed  point  above  ground,  in 
kilometers 

h2  =  height  of  the  receiving  antenna  feed  point  above  ground,  in 
kilometers 

d  =  path  length  (measured  along  the  surface  of  the  spherical  earth) , 
kilometers 

The  distances  C  and  C'  and  the  angles  Y  and  y'(as  defined  in  Figure 
and  the  angles  a,  a'  and  8  can  be  calculated  from  the  antenna  height  (h1 
h2)  and  the  path  length  ( d)  by  using  the  nine- step  procedure  given  below 

Step  1 


Determine  the  ratio  of  the  heights  (u)  from: 


Step  2 

Solve  for  S  from : 


i 


3  2  S  1  +  u 

2  S  "  2  2 


(3-9) 


in 

2) 

and 
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Step  3 

Solve  for  0  from: 


h1 

i 

,\ 

2a 

e 

—  -  sv 
sv 

Step  4 

Solve  for  a  and  a'  from: 


Step  5 

Calculate  y  and  y'  from: 

t 

y  =  tt/2  -  (a  +  0) 

y'  =  Tt/2  -  (a'  +  0) 

Step  6 

Check  to  see  whether  or  not : 

0  =  tt  -  20  -  (a  +  a') 


(3-10) 


(3-11) 


(3-12) 


(3-13) 

(3-14) 


(3-15) 


with  8  given  by  Equation  3-6 
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If  the  equality  in  Equation  3-15  is  not  satisfied,  assvane  another  value  of  3 
and  repeat  Steps  4,  5,  and  6  until  the  equality  in  Equation  3-15  is  satisfied. 

Step  7 

Calculate  C  and  C'  from: 


C 


a 


e 


sin  y 
sin  a 


(3-16) 


a 


e 


sin  y' 
sina' 


(3-17) 


Step  8 

Solve  for  Rd  ( needed  for  the  calculation  of  the  field  strength  due  to  the 
direct  wave)  from: 


-V 


b  +  b  -  2  bb'  cosQ 


(3-18) 


with  b,  b'  ,  and  9  given  by  Equations  3-4,  3-5,  and  3-6,  respectively. 


Step  9 


Solve  for  from  Equation  3-1. 


SUMMARY 


Within  the  radio  horizon,  the  divergence  factor  is  used  by  APACK  to 


account  for  the  effects  of  earth  curvature  without  using  the  Btemmer 
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formulation.  The  contribution  of  the  direct  wave  at  the  observation  point  is 
calculated  using  the  antenna  heights  and  path  length.  The  contribution  of  the 
ground-reflected  wave  at  the  observation  point  is  calculated  from  the  antenna 
heights  and  path  length  by  using  the  Fresnel  reflection  coefficient  multiplied 
by  the  divergence  factor.  Since  the  surface-wave  contribution  is  negligible 
when  the  antennas  are  several  wavelengths  or  more  above  the  ground,  the 
surface  wave  does  not  have  to  be  included  in  this  region. 
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SECTION  4 

EXTENSION  OF  THE  FORMULATION  TO  THE  DIFFRACTION  REGION 
BEYOND  THE  RADIO  HORIZON 

INTRODUCTORY  REMARKS 

The  formulation  presented  in  Section  2  for  the  far  fields  of  a  current 
element  assumes  that  the  element  is  located  above  planar  earth.  R>r 
calculations  of  field  intensities  over  a  spherical  earth  in  the  diffraction 
region  beyond  the  radio  horizon,  APACK  employs  the  Bremmer  formulation  (see 
Reference  17).  This  formulation  gives  the  far  fields  of  a  linear  antenna  in 
terms  of  the  radiation  vector  which  accounts  for  the  geometry  of  the  antenna 
and  the  aremmer  secondary  factor  which  accounts  for  the  geometry  of  the  path. 

The  radiation  vector  is  useful  not  only  in  describing  the  far  fields  in 
the  diffraction  region  but  also  for  calculating  antenna  gain.  Therefore,  the 
radiation  vector  will  be  discussed  first,  followed  by  a  discussion  of  the 
diffraction  region  far  fields. 

THE  RADIATION  VECTOR 


The  radiation  vector 


—  19  20 

(N)  for  a  linear  antenna  is  defined  by:  ' 


=  frl( s)ejksCOS'1'  ds 
L 


(4-1  ) 


1  9 

Fbster,  D. ,  "Radiation  from  Rhombic  Antenna,"  Proc.  IRE,  Vol.  25,  No.  10, 
October  1937,  pp.  1327-1353. 

20Schelkunof f ,  S.A.,  "A  General  Radiation  formula,"  Proc.  IRE,  Vol.  27,  No.  10, 
October  1939,  pp.  660-666. 
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where  L  indicates  integration  over  the  current  elements  comprising  the  linear 
antenna  and  all  other  terms  have  been  defined  in  Section  2  following  Equations 
2-11  and  2-12  (see  also  Figure  1).  In  spherical  coordinates  (r,0,$),  the 
radiation  vector  can  be  written  as : 


n 


*r  +  ae  N0  +  a4 


N . 


(4-2) 


where  a  ,  a.,  and  a,  are  unit  vectors  in  the  r-,  0-,  and  <(> -directions, 
r  o  <p 

respectively . 


The  free- space  far- field  components  are  then  given  in  terms  of  the 
radiation  vector  by : 


Efl  =  j3  Ok 


e-jkr 
—  N0 


(4-3) 


•  jkr 


Ea  =  j30k 


(4-4) 


j30k  e _ 

120*  r 


-jkr 


(4-5) 


j3  0k  e 


-jkr 


120if  r 


(4-6) 


where 


N 


0 

k 


■  0-,  component  of  the  electric  far  field,  in  volts/meter 

*  9  — ,  <^—  component  of  the  magnetic  far  field,  in  amperes/meter 

>  9-,  component  of  the  radiation  vector,  in  volt-meters 

■  —  where  X  is  wavelength,  in  meters 

radial  distance  from  the  origin  to  the  observation  point,  in  meters. 
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The  time- average  Poynting  vector  (P^)  f°un^  from: 


P  =  ^  a  R  (ExH* ) 

r  2  r  e 


(4-7) 


where 


=  time-average  Poynting  vector,  in  watts/square  meter 
a^  =  unit  vector  in  the  r-direction  in  spherical  coordinates  (r,  9,  <|>) 

E  =  electric  far  field,  in  volts/meter 
H  =  magnetic  far  field,  in  amperes/ meter 

(In  Equation  4-7,  "Re"  denotes  the  real  part  of,  and  denotes  complex 

conjugate.)  R>r  the  components  of  E  and  H  given  by  Equations  4-3  through  4-6 


'r  ■  i‘,  R.  “s'  *  h  V 


(4-8) 


Substituting  Equations  4-3  through  4-6  into  Equation  4-7,  the  magnitude 
of  the  time  average  Poynting  vector  can  be  written  in  terms  of  the  radiation 
vector  as : 


i  S  K  •».'■] 


(4-9) 


33 


where 

g^(0,<}))  =  power  gain  as  a  function  of  the  spherical  angular  coordinates 

0  and  p  (numerical  ratio) 

W^n  =  power  input  to  the  antenna ,  in  watts . 

THE  BREMMER  FORMULATION 

The  divergence  factor,  presented  in  Section  3,  extends  the  basic 
formulation  for  an  elementary  current  element  over  planar  earth  to  include  the 
effects  of  spherical  earth  within  the  radio  horizon.  Beyond  the  radio 
horizon,  diffraction  phenomena  must  be  taken  into  account. 
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The  problem  of  electromagnetic  wave  propagation  over  a  lossy  homogenous 

21  22 

spherical  earth  was  solved  by  Qremmer  (see  Reference  17)  and  others  '  many 
years  ago.  In  this  solution,  the  transmitting  antenna  is  assumed  to  be  a 
Hertzian  dipole. 

The  assumption  of  a  Hertzian  dipole  antenna  does  not  fundamentally  limit 
the  solution  obtained  because  it  is  well  known  that  the  fields  of  an  antenna 
of  finite  length  can  be  obtained  by  integration  from  the  superposition  of 
Hertzian  dipole.  The  integration  is  not  straightforward,  however,  because  the 
fields  of  a  Hertzian  dipole  located  above  lossy  spherical  earth  are  given  by 
an  infinite  series. 

In  the  spherical- earth  theory,  the  vertical  component  of  the  electric  far 
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field  of  a  Hertzian  dipole  is  given  by: 
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where 

E  =  vertical  component  of  the  electric  far  field,  in  volts/meter 
IQ  =  Hertzian  dipole  current,  in  amperes 
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i  *  length  of  the  Hertzian  dipole,  in  meters 
k  =  wave  number  (~)  ,  in  meters-1 
eQ  =  permittivity  of  free  space,  in  farads/meter 
to  =  frequency,  in  radians/ second 

d  =  distance  along  the  surface  of  the  spherical  earth,  in  meters 
r  =  radial  distance  from  the  origin  to  the  observation  point,  in  meters 
=  Bremmer  secondary  factor. 

The  Er aiuner  secondary  factor  (  Fr)  is  used  to  describe  the  far  fields  of 
linear  antennas  in  the  diffraction  region  when  the  radiation  vector  for  the 
antennas  is  known .  Fr  is  given  by  : 


( 4-1 3a ) 


where 


fs<V  “ 

fg(h2)  - 


radius  of  the  spherical  earth,  in  meters 
effective  radius  of  the  spherical  earth,  in  meters 
parameter  associated  with  the  vertical  lapse  of  the  permittivity 
of  the  atmosphere  (dimensionless) 

height  of  the  transmitting  antenna  feed  point  above  the  surface  of 
the  spherical  earth,  in  meters 

height  of  the  receiving  antenna  feed  point  above  the  surface  of 

the  spherical  earth,  in  meters 

height  gain  factor  of  the  transmitting  antenna 

height  gain  factor  of  the  receiving  antenna 
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(for  a  vertical  element) 
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The  factor  xg  is  calculated  from  Riccati1  s  differential  equation: 


d6 

■—  -  26  2  T  +1=0  (4-13d) 

dx  e  s 

s 


Cbmputational  formulas  for  evaluating  fs  (h.,),  fg  (h2),  xg,  and  6  can  be 
found  in  Reference  23,  Appendix  X.  Although  the  formulas  presented  in 
Reference  23  are  not  amendable  to  manual  computation,  they  can  be  used  for 
automated  calculations • 

As  shown  in  Reference  20,  the  free- space  field  intensities  of  linear 
antennas  other  than  the  Hertzian  dipole  can  be  obtained  by  replacing  the 
dipole  moment  in  expressions  fof  the  field  intensities  of  a  Hertzian  dipole 
with  the  radiation  vector.  attending  this  principle  to  the  diffraction 
region,  the  equations  for  the'  electric  far-field  components  of  a  linear 
antenna  in  the  diffraction  region  are  given  by: 
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(A  similar  use  of  this  principle  was  made  by  Kuebler24  in  the  line-of-sight 
region  by  substituting  the  field  function  for  linear  antennas  for  the  dipole 
moment.  Hie  field  function  is  identical  to  the  radiation  vector  except  for  a 
constant .) 

Equations  4-14  and  4-15  were  used  to  formulate  the  far  fields  of  the  16 
types  of  linear  antennas  listed  in  TABLE  1  for  the  diffraction  region  beyond 
the  radio  horizon. 


4Kuebler,  W. ,  Ground-Wave  Electric  Field  Intensity  formulas  for  Linear 
Antennas ,  ECAC-TN-74-1 1,  Electromagnetic  Compatibility  Analysis  Oenter, 
Annapolis,  MD,  June  1974. 
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SECTION  5 

CALCULATIONS  OF  FIELD  STRENGTH  AND  GAIN 


The  formulation  presented  in  Section  2  for  the  far  fields  of  a  current 
element  above  planar  earth  was  used  to  derive  the  far-field  expressions  for 
the  16  types  of  antennas  considered.  The  extensions  to  this  formulation, 
presented  in  Sections  3  and  4,  were  utilized  to  derive  the  far  fields  over 
spherical  earth  within  the  radio  horizon  and  in  the  diffraction  region  beyond 
the  radio  horizon. 

This  section  presents  a  brief  introduction  to  each  of  the  16  types  of 
antennas  (including  a  figure  of  the  antenna  geometry  and  a  description  of  the 
geometrical  parameters) ,  TABLE  2  (which  refers  to  appropriate  equations  in 
APPENDIXES  A  and  D  for  calculating  the  components  of  the  electric  far  field, 
directive  gain,  and  radiation  efficiency),  and  supplemental  symbols  and 
formulas  that  are  frequently  used  in  the  equations  listed  in  TABLE  2.  Symbols 
occurring  in  the  equations  in  APPENDIX  D  are  described  with  corresponding 
equations  in  APPENDIX  A. 

TABLE  2  serves  as  a  reference  to  key  equations  for  calculating  the 
electric  far  fields,  directive  gains,  and  radiation  efficiencies  of  the  16 
types  of  antennas.  Summaries  of  the  derivations  of  these  equations  are 
presented  in  APPENDIX  A.  The  observation  point  P  (r,0,<J>)  indicated  in  all 
figures  and  referenced  to  in  TABLE  2  is  determined  by  the  standard  spherical 
coordinates  r,  9,  and  <j>. 

HORIZONTAL  DIPOLE 


The  horizontal  dipole  (or  doublet)  antenna  shown  in  Figure  3  is  a  basic 
antenna  type  commonly  used  over  a  wide  range  of  frequencies.  The  dipole  also 
serves  as  a  building  block  element  for  antenna  arrays.  The  fields  of  the 
dipole  are  a  function  of  its  length  and  height  above  ground.  R>r  this 
antenna,  the  feed  point  (i.e.,  the  point  at  which  the  antenna  is  excited  by 
transmission  line)  is  located  at  the  center  of  the  element. 


Figure  3.  Horizontal  dipole. 
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The  geometrical  parameters  of  interest  for  the  horizontal  dipole  are: 
l  =  half-length  of  the  dipole 

H  =  height  of  the  feed  point  (i.e.,  center  of  the  dipole)  above  ground. 

VERTICAL  MONOPOLE 


The  vertical  monopole,  also  known  as  a  "whip,"  is  identical  to  half  of  a 
dipole  antenna.  The  vertical  monopole  is  fed  at  its  base  and  produces  a  field 
that  is  not  a  function  of  the  spherical  angle  <j> .  Since  the  field  is  not  a 
function  of  <f>,  the  vertical  monopole  is  referred  to  as  being  omnidirectional, 
often  shortened  to  "omni."  The  only  geometrical  parameter  needed  for  the 
vertical  monopole  is  l,  the  length  of  the  monopole.  Figure  4  shows  the 
geometry  of  the  vertical  monopole. 

VERTICAL  MONOPOLE  WITH  RADIAL-WIRE  GROUND  SCREEN 


Most  vertical  monopoles  that  are  permanently  installed  include  a  ground 
screen  to  improve  the  radiation  resistance  and  to  increase  the  level  of  the 
fields  radiated  at  small  values  of  elevation  angle  (i.e.,  values  of  0  near 
90°).  One  form  of  ground  screen  commonly  used  consists  of  a  group  of  radial 
wires  placed  on  the  ground  that  are  centered  at  the  base  of  the  monopole  and 
spaced  equally  in  angle. 

The  geometrical  quantities  of  interest  for  the  vertical  monopole  with 
radial-wire  ground  screen  (see  Figure  5)  are: 

l  =  length  of  the  monopole 

C  =  radius  of  the  wires  comprising  the  ground  screen 

a  =  radius  of  ground  screen 

N  =  number  of  radial  wires . 
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ELEVATED  VERTICAL  DIPOLE 


The  elevated  vertical  dipole  antenna  is  simply  a  dipole  that  is  oriented 
with  its  axis  orthogonal  to  the  ground  plane  below  it.  The  fields  of  this 
antenna  are  not  a  function  of  the  spherical  angle  <t>,  so  the  elevated  vertical 
dipole  is  also  considered  "omni 

The  pertinent  geometrical  quantities  shown  in  Figure  6  are: 

£  =  half-length  of  the  dipole 

ZQ  =  height  of  the  feed  point  (i.e.,  center  of  the  dipole)  above  ground. 

INVERTED -L 


The  inverted-L  antenna,  shown  in  Figure  7,  consists  of  both  a  vertical 
and  a  horizontal  section.  This  antenna  is  fed  at  the  bottom  of  the  vertical 
section.  The  inverted-L  is  often  used  when  a  tall  vertical  monopole  is 
inconvenient  to  erect,  because  the  horizontal  section  acts  as  a  "top  load"  for 
the  vertical  section.  This  effectively  increases  the  length  of  the  vertical 
section. 

The  geometrical  parameters  for  the  inverted-L  are  H,  the  length  of  the 
vertical  section,  and  l,  the  length  of  the  horizontal  section. 

ARBITRARILY  TILTED  DIPOLE 


The  arbitrarily  tilted  dipole  is  a  dipole  that  is  inclined  at  an  angle 
with  respect  to  the  ground  under  it.  The  geometrical  parameters  as  shown  in 
Figure  8  are: 

l  =  half-length  of  the  dipole 

H  =  height  of  the  feed  point  (i.e.,  center  of  the  dipole)  above  ground 
a'  -  angle  between  the  axis  of  the  dipole  and  the  Y-axis  . 
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SLOPING  LONG-WIRE 


The  sloping  long-wire,  shown  in  Figure  9,  is  fed  at  the  base  and  produces 
both  vertically  and  horizontally  polarized  field  components  that  depend  on  the 
inclination  of  the  wire  with  respect  to  the  ground  plane.  The  conductivity  of 
the  ground  beneath  the  wire  can  have  substantial  effects  on  the  radiation 
characteristics.  Wien  the  ground  has  high  conductivity,  radiation  is 
reflected  off  the  ground  in  the  direction  of  the  high  end  of  the  wire.  Wien 
the  ground  has  low  conductivity,  radiation  directed  toward  the  ground  is 
absorbed,  and  radiation  directed  upward  and  in  the  direction  opposite  to  the 
high  end  predominates. 

The  geometrical  parameters  of  interest  for  the  sloping  long-wire  are : 

l  =  length  of  the  wire 

a'  =  angle  between  the  axis  of  the  wire  and  the  Y-axis. 

TERMINATED  SLOPING-V 


■Hie  terminted  sloping-V,  shown  in  Figure  10,  consists  of  two  wires  fed  at 
the  apex  and  terminated  with  appropriate  resistances  at  the  ends  away  from  the 
feed  point.  The  terminating  resistances  on  the  wires  cause  the  currents  in 
the  wires  to  result  in  traveling  waves.  Thus,  the  terminated  sloping-V  is  a 
traveling-wave  (i.e.,  nonresonant)  antenna  as  opposed  to  other  types  discussed 
previously  which  are  resonant  antennas.  Traveling-wave  antennas  have  the 
advantage  of  providing  operation  over  a  wide  range  of  frequencies  without  the 
need  for  matching  (coupling)  networks  between  the  feed  point  and  attached 
transmission  line. 

The  fields  of  the  terminated  sloping-V  depend  on  a  number  of  parameters 
including  the  wire  lengths,  angle  between  the  wires,  and  heights  of  the 
structure.  The  geometrical  parameters  are: 


48 


ESD-TR-80-1 02 


Section  5 


2 


B] 


ESD-TR-80- 1 02 


Section  5 


l  =  length  of  the  wire 
Y  =  half-angle  between  the  wires 

H  =  height  of  the  apex  (feed  point)  above  ground 

H'  =  height  of  the  terminated  end  of  the  wires  above  ground 

a'  =  angle  between  the  plane  containing  the  wires  and  the  X-axis 

6'  =  angle  between  the  projection  of  the  wires  in  the  XY-plane  and  the 

X-axis . 

TERMINATED  SLOPING  RHOMBIC 


The  terminated  sloping  rhombic,  shown  in  Figure  11,  can  be  considered  as 
being  made  up  of  two  sloping-V  antennas  placed  end-to-end.  Hie  apex  of  one  of 
the  sloping-V  antennas  is  used  as  the  feed  point,  and  the  apex  of  the  other 
sloping-V  antenna  is  terminated  in  an  appropriate  resistance.  Hie  terminated 
rhombic  is  also  a  traveling-wave  antenna  that  can  be  operated  over  a  wide 
range  of  frequencies  without  matching  (coupling)  networks. 

Hie  geometrical  parameters  are : 

1  =  length  of  each  of  the  four  wires  comprising  the  rhombus 

Y  =  half  angle  between  the  wires,  the  feed  point,  and  the  termination 
H  =  height  of  the  feed-point  apex  above  ground 

H'  =  height  of  the  terminated  apex  above  ground 

H”  =  height  of  the  center  of  the  rhombus  above  ground 
a'  =  tingle  between  the  plane  containing  the  rhombus  and  the  x-axis 

8'  »  angle  between  the  projection  of  the  feed- point  apex  in  the  XY-plane 

and  the  x-axis. 

TERMINATED  HORIZONTAL  RHOMBIC 


The  terminated  horizontal  rhombic  is  identical  to  the  terminated  sloping 
rhombic  ex-  pt  that  the  plane  of  the  rhombus  is  parallel  to  the  XY-plane.  The 
geometrical  parameters  for  the  horizontal  rhombic,  shown  in  Figure  12,  are 
identical  to  those  for  the  sloping  rhombic  except  that  H  is  the  height  of  the 
rhombus  above  ground. 
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SIDE-LOADED  VERTICAL  HALF-RHOMbIC 

The  side- loaded  vertical  hal c- rhombic ,  shown  in  Figure  13,  consists  of 
two  sectons  of  wire  fed  at  one  end  and  terminated  with  an  appropriate  resistor 
at  the  other  end.  The  side- loaded  vertical  half- rhombic  is  also  a  traveling- 
wave  antenna  and  radiates  a  combination  of  vertically  and  horizontally 
polarized  fields  depending  on  the  length  of  the  wires  and  their  inclination 
with  respect  to  the  ground  plane.  Note  that  if  the  ground  plane  were 
perfectly  conducting,  the  two  wire  sections  and  their  images  below  the  ground 
would  form  a  rhombic  antenna  in  the  YZ-plane. 

The  geometrical  parameters  of  interest  for  the  side- loaded  vertical  half- 
rhombic  are : 

i  =  length  of  each  of  the  wire  sections 
a'  =  angle  between  the  feed  point  or  termination  point  of  the  wires  and  the 
Y-axis. 

HORIZONTAL  YAGI-UDA  ARRAY 


Hie  horizontal  Yagi-Uda  array  is  a  coplanar  arrangement  of  dipole 
elements  of  different  lengths  with  variable  spacing  between  dipoles. 
Specifically,  a  Yagi-Uda  array  consists  of  a  single  driven  element  and  one  or 
more  parasitic  elements  that  can  function  either  as  reflectors  or  as  directors 
(  see  Figure  14)  . 

The  director  element  is  directly  coupled  to  the  transmission  line  at  the 
center  of  the  element.  The  parasitic  elements  are  not  coupled  directly  to  the 
transmission  line  but  are  rather  coupled  through  the  electromagnetic  fields 
emanating  from  the  driven  element.  The  parasitic  reflector  element  is  longer 
than  the  driven  element  and  is  placed  behind  the  driven  element  with  respect 
to  the  desired  direction  of  radiation.  The  parasitic  director  elements  are 
shorter  than  the  driven  element  and  are  placed  in  front  of  the  driven  element 
with  respect  to  the  desired  direction  of  radiation. 
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A  common  arrangement  is  Cor  the  reflector  to  be  slightly  longer  than  a 
half-wavelength,  the  driven  element  about  a  half-wavelength,  and  the  directors 
less  than  a  half-wavelength .  Typical  spacings  between  the  dipole  elements 
vary  from  about  one-tenth  to  four- tenths  of  a  wavelength. 

Hie  geometrical  parameters  for  a  Yagi-Uda  array  are  : 

=  half-length  of  the  i1"*1  element 

U. 

d^  =  spacing  between  the  itn  element  and  the  (i  +  I)1"  element 
H  =  height  of  the  array  above  the  ground. 

HORIZONTALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 

The  horizontally  polarized  log-periodic  dipole  array  is  also  a  coplanar 
arrangement  of  dipole  elements  with  varying  lengths  and  spacing  between 
elements.  The  log-periodic  dipole  array  differs  from  the  Yagi-Uda  array  in 
that  the  element  lengths,  spacings,  and  excitation  are  designed  to  provide 
operation  over  a  wide  range  of  frequencies  (typically  two  or  three  octaves) 
without  the  need  for  matching  (coupling)  networks  between  the  feed  point  and 
attached  transmission  line. 

The  geometrical  parameters  for  the  horizontally  polarized  log-periodic 
dipole  array  (see  Figure  15)  are  as  follows: 

=  half-length  of  the  ith  element 
d^  =  spacing  between  the  ith  element  and  the  (i  +  1)th  element 
H1  =  height  of  the  shortest  element  above  the  Y-axis 
HN  =  height  of  the  longest  element  above  the  Y-axis 
<x  =  angle  between  the  array  axis  and  the  line  connecting  the  tips  of  the 
elements 

9"  =  angle  between  the  Z-axis  and  the  array  axis. 
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VERTICALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 


The  vertically  polarized  log-periodic  dipole  array  is  the  same  basic 
structure  as  the  horizontally  polarized  log-periodic  dipole  array  except  that 
the  plane  containing  the  array  is  vertical.  The  geometrical  quantities  for 
this  antenna,  shown  in  figure  16,  are  as  follows: 

=  half-length  of  the  ifc^  element 

i,U  tU 

d^  =  spacing  between  the  i  n  element  and  the  (i  +  1  )'"n  element 

=  height  of  the  center  of  the  shortest  element  above  the  Y-axis 
HN  =  height  of  the  center  of  the  longest  element  above  the  Y-axis 
a.j  =  angle  between  the  array  axis  and  the  line  connecting  the  upper  tips  of 
the  elements 

(*2  =  angle  between  the  array  axis  and  the  line  connecting  the  lower  tips  of 

the  elements 

=  angle  between  the  line  connecting  the  lower  tips  of  the  elements  and 
the  Y-axis 

0"  =  angle  between  the  Z-axis  and  the  array  axis. 

CURTAIN  ARRAY 


The  curtain  array,  shown  in  Figure  17,  is  a  coplanar  arrangement  of 
dipole  elements  with  spacings  and  excitations  such  that  the  direction  of 
maximum  far- field  strength  can  be  varied  in  both  the  vertical  plane  (XZ-plane) 
and  the  horizontal  plane  (XY-plane) .  the  array  of  dipoles  is  installed  in 
front  of  a  wire  screen  to  provide  a  unidirectional  characteristic.  The 
dipole-to-screen  spacing  is  usually  less  than  a  half-wavelength. 

The  dipoles  are  arranged  in  bays  of  vertically  stacked  elements  with 
variable  spacing.  The  bays  of  elements  are  horizontally  separated  by  variable 
distances. 


ESD-TR-80-102 


Section  5 


Ihe  geometrical  quantities  of  interest  for  the  curtain  array  are: 

l  =  half-length  of  a  dipole  element 
X-j  =  spacing  between  the  bays  of  dipoles  and  the  wire  screen 

4"  Vi 

y^  ==  horizontal  position  of  the  center  line  of  the  (i  +  1)  bay  (where 

y0  A  0) 

Z^  =  height  of  the  ith  element  of  each  bay  above  ground. 

M  =  number  of  elements  in  each  bay 
N  =  number  of  bays . 

SLOPING  DOUBLE  RHOMBOID 


The  sloping  double  rhomboid,  shown  in  Figure  18,  consists  of  two  sloping 
rhombics  fed  at  a  common  apex  and  terminated  at  the  opposite  apexes  with 
appropriate  resistors.  The  sloping  double  rhomboid  is  also  a  traveling-wave 
antenna  that  can  be  used  over  a  wide  range  of  frequencies  without  the  need  for 
matching  (coupling)  networks  between  the  feed  point  at  the  common  apex  and  the 
attached  transmission  line.  Wien  the  geometry  for  this  antenna  is  carefully 
chosen,  the  result  is  a  pattern  with  higher  gain  in  the  mainbeam  direction  and 
lower  gain  in  the  sidelobe  regions  that  could  be  realized  with  a  single 
rhombic . 

The  geometrical  parameters  for  the  sloping  double  rhombiod  are : 

=  length  of  the  short  legs  of  each  rhombus 
&2  =  length  of  the  long  legs  of  each  rhombus 

a  =  angle  between  the  principal  antenna  axis  (Y-axis)  and  the  major 
rhomboid  axis 

8  =  angle  between  the  major  rhomboid  axis  and  the  long  leg 

Y  =  angle  between  the  major  rhomboid  axis  and  the  short  leg 

=  height  of  the  feed  point  at  the  common  apexes  above  the  Y-axis 

Hj  =  height  of  the  apexes  at  the  ends  of  the  short  legs  emanating  from  the 

feed  point  above  the  Y-axis 


Sloping  double 


rhomboid. 
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=  height  of  the  apexes  at  the  ends  of  the  long  legs  emanating  from  the 
feed  point  above  the  Y-axis 
=  height  of  the  terminated  apexes  above  the  Y-axis 
a'  *  angle  between  the  plane  containing  the  two  rhomboids  and  the  Y-axis. 

KEY  EQUATIONS 

To  facilitate  the  use  of  this  report  as  a  reference  for  key  equations  for 
each  of  the  16  types  of  antennas  considered,  TABLE  2  was  prepared.  TABLE  2 
indicates  the  appropriate  equations  in  APPENDIXES  A  and  D  that  give  the 
electric  far-f ield  components  in  each  of  the  three  regions  considered  : 

1.  Over  planar  earth 

2.  Over  spherical  earth  within  the  radio  horizon 

3.  In  the  diffraction  region  beyond  the  radio  horizon. 

Hie  criteria  used  to  determine  in  which  of  the  three  regions  the 
observation  point  is  located  are  given  in  Section  6  under  the  subsection, 
"Changeover  Criteria."  Hie  equations  for  the  electric  far-f ield  components 
indicated  in  TABLE  2  together  with  the  changeover  criteria  given  in  Section  6 
have  been  used  at  ECAC  to  write  a  computer  program  that  calculates  the 
electric  far  fields  of  the  16  types  of  antennas  over  lossy  ground  at  any  point 
in  the  far  field. 

TABLE  2  also  indicates  appropriate  equations  in  APPENDIX  A  for 
calculating  the  directive  gain  of  each  of  the  16  types  of  antennas  as  well  as 
appropriate  equations  or  values  for  radiation  efficiency.  Note  that  since  the 
equations  given  in  APPENDIX  A  for  directive  gain  are  expressed  as  numerical 
ratios,  the  directive  gain  in  dB  is  given  by: 


Gd(0,4>)  =  10  log1Q  gd  <0,4) 


(5-1  ) 
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where 

<  0 ,  ((> )  =  directive  gain  as  a  function  of  the  spherical  angular  coordinates 
9  and  <}>,  in  dB 

gd(  0  ,<)>)  =  directive  gain  as  a  function  of  the  spherical  angular  coordiantes 
9  and  <J>  as  a  numerical  ratio  ( as  given  by  the  equations  in 
APPENDIX  A) . 

Also,  the  power  gain  is  calculated  simple  from: 

G  (0,4.)  =  G,(0,4>)  +  n  (5-2) 

p  d  dB 

where 

Gp(9,4>)  -  power  gain  as  a  function  of  the  spherical  angular  coordinates  0 
and  41/  in  dB 

ndB  =  radiation  efficiency,  in  dB. 

The  formulas  for  directive  gain  and  radiation  efficiency  (presented  in 
TABLE  2)  and  the  formulas  for  transmission  loss  (presented  in  Section  6)  have 
also  been  programmed  at  ECAC.  Ihus,  APACK  provides  a  versatile  automated 
package  for  calculating  far- field  intensity,  directive  gain,  power  gain,  and 
transmission  loss  for  16  types  of  antennas  mounted  over  lossy  ground . 

The  equations  indicated  in  TABLE  2  are  the  primary  equations  used  in 
calculating  the  various  quantities.  Secondary  equations  defining  additional 
geometrical  quantities  and  terms  such  as  the  radiation  resistance  are  given  in 
APPENDIX  A  in  the  same  subsection  as  the  primary  equations.  Other  symbols  and 
formulas  that  appear  often  in  the  primary  equations  are  presented  below. 
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SUPPLEMENTAL  SYMBOLS  AND  FORMULAS 


The  following  supplemental  symbols  and  formulas  are  used  with  the 
equations  indicated  in  TABLE  2  to  calculate  field  intensities  and  gains. 


Symbols  frequently  used  are: 


k 


2tt 

X 


free-space  wave  number,  in  (meters) 


where 


X 

k2 

adiv 

Nr  ,  9  ,  <{) 


d 


R 


d 


9 


d 


wavelength,  in  meters 

wave  number  for  the  lossy  medium  under  the  antenna,  in  (meters)-1 
divergence  factor  used  for  calculations  over  spherical  earth  within 
the  radio  horizon,  given  by  Equation  3-1 

r-,9-,  or  ^-component  of  the  radiation  vector  used  for  calculations 
in  the  diffraction  region  beyond  the  radio  horizon 

Bremmar  secondary  factor  used  for  calculations  in  the  diffraction 
region  beyond  the  radio  horizon,  given  by  Equation  4-1 3a 
distance  along  the  surface  of  a  spherical  earth  between  the  antenna 
location  and  the  observation  point 

distance  from  the  feed  point  of  the  antenna  to  the  observation 
point 

angle  between  the  Z-axis  and  the  line  between  the  antenna  feed 
point  and  the  observation  point 

distance  from  the  feed  point  of  the  image  antenna  and  the 
observation  point 

angle  between  the  Z-axis  and  the  line  between  the  image  antenna 
feed  point  and  the  observation  point. 


Rarmulas  frequently  used  are: 


n 


_2 

k 


8000a 

fMHz 


(5-3) 
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where 


n 


e 


r 


a 

f  MHz 


complex  index  of  refraction  for  the  lossy  medium  under  the  antenna 
( dimensionless) 

relative  dielectric  constant  of  the  lossy  medium  under  the 
antenna  (dimensionless) 

conductivity  of  the  lossy  medium  under  the  antenna,  in  mhos/meter 
frequency,  in  MHz. 


-  jkR 


P 

e 


2  sin  0 


cos  9  + 

r 


yP~- 


2 

sin  6 


(5-4) 


where 


Pg  =  complex  numerical  distance  for  vertical  polarization  (dimensionless) 


and  all  other  terms  have  been  defined  previously. 


P 

m 


-jkR 

r 

2  sin20 

r 


(cos9  + 
r 


'  5 

sin  9 

r 


(5-5) 


where 


Pm  =  complex  numerical  distance  for  horizontal  polarization 
( dimensionless) 


and  all  other  terms  have  been  defined  previously. 
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F 

e 


-P 

-  e 

1-j  irP  e 
n  e 


erfc 


(5-6) 


where 


Fg  =  surface-wave  attenuation  function  for  vertical  polarization 
( dimensionless) 


and 


erfc  =  complementary  error  function  defined  by: 


2 

2  x  -3 

erfc(x)  =  1  -  —  /  e  dy  (5-7) 

■Jit  0 


-P 


where 

^  =  surface-wave  attenuation  function  for  horizontal  polarization 

( dimensionless) 


R 

v 


2 

n 


cos  9 

r 


2 

n 


cos  9 

r 


sin2  0 


•  2  o 
srn  0 


(5-9) 
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where 

Ry  =  ftesnel  plane-wave  reflection  coefficient  for  vertical 
polarization  (dimensionless) 


Rjj  =  Etesnel  plane-wave  reflection  coefficient  for  horizontal 
polarization  (dimensionless). 
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SECTION  6 

CALCULATIONS  OF  TRANSMISSION  LOSS 


The  basic  calculations  performed  by  APACK  are  those  for  the  electric  far 
fields  of  a  linear  antenna  over  planar  earth,  over  spherical  earth  within  the 
radio  horizon,  and  over  spherical  earth  in  the  diffraction  region  beyond  the 
radio  horizon.  Since  the  electric  far  fields  are  calculated,  APACK  can  be 
used  to  calculate  transmission  loss  as  well  as  directive  and  power  gains. 

This  section  presents  expressions  used  for  calculating  transmission  loss, 
a  discussion  of  field- strength  curves  presented  by  the  International  Radio 
Consultative  Cbmmittee  (CCIR),  and  a  discussion  of  the  changeover  criteria 
used  to  calculate  field  strengths  in  different  regions. 


TRANSMISSION-LOSS  EXPRESSIONS 


pc 

Transmission  loss  (L)  is  defined  by: 


L  =  PT  -  PR  (6-1) 

where 

L  =  transmission  loss  between  the  transmitting  antenna  and  the  receiving 
antenna,  in  dB 

PT  =  total  power  radiated  from  the  transmitting  antenna,  in  dan 
PR  =  resultant  signal  power  that  would  be  available  from  an  equivalent 
loss-free  antenna,  in  dBn. 


2  5 

Rice,  P.L. ,  Iongley,  A.G. ,  Norton,  K.A.,  and  Earsis,  A.P. ,  Transmission 
Loss  Predictions  for  Tropospheric  Cbmmunication  Circuits,  NBS  Technical 
Note  101,  Vols.  I  and  II,  National  Bureau  of  Standards,  Boulder,  CO,  1967. 
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Transmission  loss  can  also  be  expressed  as : 


L  =  V  +  V  +  \  +  \  (6_2) 

where 

P  '  =  power  input  to  the  feed  point  of  the  transmitting  antenna,  in  dBn 

P_"  =  resultant  signal  power  available  at  the  feed  point  of  the  receiving 

antenna,  in  dBn 

nT  =  radiation  efficiency  of  the  transmitting  antenna,  in  dB 
h  =  radiation  efficiency  of  the  receiving  antenna,  in  dB. 

R 

(Note  that  nT  and  nR  are  always  less  than  or  equal  to  0  dB  because  radiation 
efficiency  is  always  less  than  or  equal  to  100%.) 

Basic  transmission  loss  (1^)  which  is  the  loss  between  two  lossless 
isotropic  antennas  is  defined  in  Reference  19  as : 


L  +  G 

P 


(6-3) 


where 


=  basic  transmission  loss,  in  dB 
Gp  =  path  antenna  gain,  in  dBL. 

Wien  multipathing  and  antenna  efficiencies  are  taken  into  account,  the 
path  antenna  gain  is  not  the  sum  of  the  directive  gains.  Thus,  in  general: 
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GTCl  =  directive  gain  of  the  transmitting  antenna,  in  dBi 
GR(J  =  directive  gain  of  the  receiving  antenna,  in  dBi. 

The  basic  transmission  loss  can  also  be  expressed  as  : 

h  "  Sf  +  ADB 

where 


(6-4) 


I^f  =  basic  transmission  loss  in  free  space,  in  dB 

ADB  =  attenuation  relative  to  free-space  loss  (where  all  difficult  problems 
in  analyzing  propagation  phenomena  are  hidden) ,  in  dB. 

The  basic  free-space  transmissin  loss  (L^)  can  also  be  expressed  as  : 


Sf  =  20  logi( 


=  32.45  +  20  log. _  f  +  20  log„„  d 

^10  MHz  10  km 


(6-6) 


where 

d  =  path  length,  in  meters 
\  =»  wavelength,  in  meters 

f MHz  =  frequency,  in  MHz 
dkm  =  path  length,  in  kilometers. 
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5tom  the  above  equations : 


L  =  +  ADB  -  Gp 

=  32.45  +  20  log  f  +  20  log,„  d,  +  ADB  -  G  (6-7) 

10  MHz  10  Km  p 

where,  in  this  case 


G 

P 


M 


+  G 


Rd 


(6-8) 


The  resultant  signal  power  available  at  the  feed  point  of  the  receiving 
antenna  (Pr’)  is  then  given  in  terms  of  the  power  input  to  the  feed  point  of 
the  transmitting  antenna  (PT* ) ,  directive  gains  (G,^,  G^)  ,  and  radiation 
efficiencies  (nT,  nR)  by s 


PD  =  P_  -  32.45  -  20  log  f  -20  log„„  d,  -ADB 
R  T  ^10  MHz  ^10  km 


+  +  G„.  +  n  +  n 

Td  R3  T  R 


(6-9) 


where 


ADB 


20  log 

r10 


(6-10) 


with 
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_  | Actual  disturbed  field  at  observation  point | 

|  Rree-space  field  at  observation  point | 


(6-11  ) 


1E0|  =  magnitude  of  the  0-component  of  the  actual  disturbed  field  at  the 

observation  point,  in  volts/meter 

|Egf|  =  magnitude  of  the  0-component  of  the  fre^-space  field  at  the 
observation  point,  in  volts/meter 

I E^f |  =  magnitude  of  the  ^-component  of  the  free-space  field  at  the 
observation  point,  in  volts/meter 

I E^l  =  magnitude  of  the  ^-component  of  the  actual  disturbed  field  at  the 
observation  point,  in  volts/meter. 

Equations  6-9,  6-10,  and  6-12  are  the  expressions  used  by  APACK  to  calculate 

transmission  loss. 

CCXR  CURVES 

Various  propagation  models  have  been  used  at  ECAC  to  predict  ground  wave 

loss  over  a  spherical  earth.  These  models  include  curves  published  by  the 
26 

CCXR  and  other  models  referred  to  at  ECAC  as  IPS  and  NA  (see  References  3 


International  Radio  Consultative  (Committee  (CCIR),  Recommendations  and 
Reports  of  the  CCIR,  1978,  Propagation  in  Non- Ionized  Media,  Vol.  V,  XIVth 
Plenary  Assembly,  Kyoto,  Japan,  1978. 
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and  4,  respectively).  The  CCIR  curves  which  are  based  on  Bremmer's 
formulation  (see  Reference  17)  assume  that  both  the  transmitting  and  receiving 
antenras  are  located  on  the  ground  and  that  both  antennas  are  Hertzian 
dipoles,  lhe  CCIR  curves  based  on  these  assumptions  are  presented  in  terms  of 
field  strength  which  can  be  converted  to  basic  transmission  loss  (see 
Reference  26)  by: 


L(dB)  =  132.45  +  20  log,n  f  -  E  (6-13) 

1 0  MHz 

where 

L(dB)  =  basic  transmission  loss,  in  dB 
f Mjjz  =  frequency,  in  MHz 

E  =  electric  field  strength,  in  dB,  with  respect  to  1  microvolt/meter 
( dB  ( uv/m) )  . 

The  other  models  used  at  ECAC  (see  References  3  and  4)  predict  basic 
transmission  loss  without  calculating  field  strength  but  do  account  for  the 
heights  of  the  antenna  feed  points  above  ground.  APACK  predicts  basic 
transmission  loss  from  calculated  values  of  electric  field  strength.  the 
APACK  field  strengths  account  for  the  actual  antenna  structure  and  include  the 
effects  of  the  contributions  of  the  direct,  ground-reflected,  and  surface 
waves . 

CHANGEOVER  CRITERIA 


The  calculations  of  electric  far  fields  made  by  APACK  depend  on  the 
region  in  which  the  observation  point  is  located.  Three  regions  are 
considered : 

Region  1.  The  line-of- sight  region  when  the  antenna  is  near  the  ground 
so  that  the  earth  can  be  considered  as  planar. 
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Region  2.  ’Hie  line-of-sight  region  when  the  antenna  is  well  above  the 
ground  so  that  the  earth  must  be  considered  as  spherical. 

Region  3.  ltie  diffraction  region  beyond  the  radio  horizon. 

Hie  calculations  include  the  surface-wave  terms  in  Region  1  but  not  in 
Region  2,  because  the  surface  wave  is  negligible  when  the  antenna  is  well 
above  the  ground.  Hie  calculations  made  for  Region  3  (the  diffraction  region) 
can  be  used  for  all  three  regions.  However,  since  a  large  number  of  terns  are 
required  for  convergence  of  =he  Bremmer  series  within  the  radio  horizon,  APACK 
uses  the  Bremmer  formulation  only  in  the  diffraction  region  (Region  3). 

Hie  APACK  criteria  for  determining  the  appropriate  region  for  the 
calculations  have  been  adopted  from  Reference  22  and  are  described  below.  Hie 
region  in  which  the  observation  point  is  located  is  determined  by  the  path 
length,  the  effective  radius  of  the  earth  (in  terms  of  atmospheric 
refractivity) ,  the  heights  of  the  feed  points  of  the  transmitting  and 
receiving  antennas  above  ground,  and  frequency. 

Riur  quantities  (d,  dc,  dD,  and  M)  are  calculated,  and  the  appropriate 
region  is  then  selected  from  the  following  relationships  among  these 
quantities. 

1.  If  d  _<  dc  and  dc  >  MdD,  then  the  observation  point  is  in  Region 

1  (line-of-sight,  planar  earth). 

2.  If  d  £  Hip  and  dc  <_  Mi^,  then  the  observation  point  is  in  Region 

2  (line-of-sight,  spherical  earth). 

3.  If  ( d  >  dc  and  dc  >  MdD)  or  (d  >  MdD  and  dc  <_  M3D)  ,  then  the 
observation  point  is  in  Region  3  (diffraction  region). 

Hie  quantities  d,  dc,  dD,  and  M  are  defined  as  follows: 

d  =  path  length,  in  kilometers 


ESD-TR-80-1 02 


Section  6 


dc< 


1.0195 

148  *  ,0'4  im  •  °-,5i  <W  <  ’ 

f  MHz 


1.0195 

-4  ae 

=  148  x  10  - — r  ,  1  <  f  <  100 

0.5305  -  MHz 

fMHz 


1.0195 

-4  3e 

59.5x10  — — „  — —  ,  100  <  f  <  500 

0.3337  -  MHz  - 


MHz 


where 


(6-14) 


(6-15) 


(6-16) 


a 

e 


6370 


1-0 


.  04665e 


0.005577N  1 
s 


-1 


(6-17) 


-0.1057  h 

N  =  N  e  S  (6-18) 

3  O 


with 

ag  =  effective  earth  radius,  in  kilometers 

Ng  *  surface  atmospheric  refractivity ,  in  N-units 

Nq  =  surface  atmospheric  refractivity  reduced  to  sea  level,  in  N-units 
hg  =  elevation  of  the  surface  above  mean  sea  level,  in  kilometers. 
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(If  NQ  and  hg  are  not  given,  Ng  =  301  is  assumed.) 


d_  = 


2k  +  k  +  2  a  h.  xl  0 


V* 


+  2  a  xl  0 


-3 


(6-19) 


where 


k 


a 

e 


X  x10~3  3 

2ira 

e 


(6-20) 


h1  =  height  of  the  transmitting  antenna  feed  point  above  ground,  in  meters 
h2  =  height  of  the  receiving  antenna  feed  point  above  ground,  in  meters 
X  =  wavelength,  in  meters. 


I1*0'  fMHz  —  ^ 


(6-21  ) 

(6-22) 
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SECTION  7 

COMPARISONS  BETWEEN  GAINS  PREDICTED  BY  APACK 
AND  OTHER  DATA 


R>r  each  of  the  16  types  of  linear  antennas  considered  in  APACK, 
predictions  of  APACK  power  gains  as  a  function  of  elevation  (vertical)  and/or 
azimuth  (horizontal)  angle  were  made  for  antennas  mounted  over  soil  with 
various  permittivities  and  conductivities  and/or  over  sea  water.  The  APACK 
predictions  were  compared  with  other  available  data,  including  manufacturer's 
data,  predictions  from  the  SKYWAVE  computer  model  developed  by  ITS  (see 
Reference  1 ) ,  and  predictions  resulting  from  Ma1 s  three-term  current 
distribution  (see  Reference  14). 

APACK  gain  predictions  were  also  compared  with  gain  predictions  resulting 
from  modeling  the  antennas  with  the  Numerical  Electromagnetic  Cbde  (NEC) ,  a 
method-of-moments  computer  program  developed  at  the  Lawrence  Livermore 
Laboratory  under  the  sponsorship  of  the  Naval  Ocean  Systsns  Center  (NOSC)  and 
the  Air  R>rce  weapons  laboratory  (AFWL)  (see  Reference  2).  In  most  cases, 
there  is  reasonable  agreement  between  the  APACK  predictions  and  the  NEC 
predictions.  This  indicates  that  APACK,  based  on  sinusoidal  current 
distribution  for  resonant  antennas  and  exponential  current  distribution  for 
traveling-wave  antennas,  provides  reasonable  predictions  and  at  the  same  time 
large  savings  in  computer  time.  Rjr  example,  the  time  required  to  obtain 
gains  for  an  electrically  large  antenna  using  APACK,  such  as  a  sloping  double 
rhomboid,  can  be  approximately  1/1 000th  of  the  time  required  by  the  rigorous 
NEC  method-of-moments  program. 

The  comparisons  provided  in  this  section  are  not  exhaustive,  but  they  do 
provide  an  indication  of  the  versatility  and  reasonableness  of  the  APACK  model 
for  providing  rapid  predictions  of  gains.  Comparisons  performed  for  the 
various  types  of  antennas  are  described  below  and  shown  in  the  accompanying 
figures . 


Q7 
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HORIZONTAL  DIPOLE 


Figures  19  through  22  show  comparisons  between  APACK,  NEC,  and 
manufacturer's  data  for  the  Osllins  637T-1/2  half-wave  horizontal  dipole 
mounted  above  soil  for,  frequencies  of  5,  10,  20,  and  30  MHz,  respectively.  In 

all  cases,  the  gains  predicted  by  APAOC  differ  by  less  than  approximately  2  dB 
from  the  other  data. 

VERTICAL  MONOPOLE 


Figure  23  shows  comparisons  between  APACK  and  NEC  for  a  quarter-wave 
vertical  monopole  (without  a  ground  screen)  mounted  above  soil.  The  gains 
predicted  by  APACK  differ  from  those  predicted  by  NEC  by  less  than 
approximately  1  dB.  The  APACK  predictions  shown  in  Figure  23  are  identical  to 
those  of  the  ITS  SKYWAVE  program  (see  Reference  1). 

VERTICAL  MONOPOLE  WITH  RADIAL-WIRE  GROUND  SCREEN 


Figure  24  shows  comparisons  between  APACK,  NEC,  and  Ma's  three-term 
current- distribution  predictions  for  a  vertical  monopole  with  radial-wire 
ground  screen  mounted  above  soil  and  sea  water  at  an  operating  frequency  of 
10  MHz.  The  differences  between  APACK  predictions  and  Ma’s  predictions  are 
less  than  approximately  1  dB  for  both  soil  and  sea  water.  Hie  NEC  predictions 
are  identical  to  Ma's  for  the  monopole  mounted  over  sea  water.  Wien  the 
monopole  is  mounted  over  soil,  the  NEC  predictions  of  gain  are  considerably 
higher  (as  much  as  about  3  dB)  than  either  APACK  or  Ma.  The  reason  for  the 
higher  gains  predicted  by  NEC  has  not  been  determined. 

ELEVATED  VERTICAL  DIPOLE 


Figure  25  shows  comparisons  of  predictions  made  by  APACK,  NEC,  and  the 
ITS  SKYWVVE  program  for  a  vertical  half-wave  dipole  with  its  feed  point 
located  2.5  meters  above  soil  at  a  frequency  of  30  MHz.  The  differences 
between  the  APACK  predictions  and  the  NEC  predictions  are  less  than  5  dB  for 
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Figure  20.  Elevation  patterns  of  a  Collins  637T-1/2  half-wave 

horizontal  dipole  mounted  7.6  m  above  soil  (10  MHz). 
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Figure  21.  Elevation  patterns  of  a  Obllins  637T-1/2  half-wave 
horizontal  dipole  mounted  7.6  m  above  soil  (20  MHz) . 
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FREQUENCY :  AS  INDICATED 
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<r  =  0.01  mho/m 


23.  Elevation  patterns  of  a  quarter-wave  vertical 
monopole  mounted  above  soil. 
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FREQUENCY:  10  MHz 
SEA  WATER:  6r=80 
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Figure  24.  Elevation  patterns  of  a  vertical  monopole  with 

radial-wire  ground  screen  mounted  above  sea  water 
and  soil. 
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Figure  25.  Elevation  patterns  of  a  half-wave  vertical  dipole 
mounted  with  feed  point  2.5  m  above  soil. 
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low  elevation  angles  and  approach  7  dB  near  the  zenith.  The  gains  predicted 
by  the  ITS  SKYWAVE  program  are  considerably  less  than  those  predicted  by 
either  APACK  or  NEC.  The  low  values  of  SKYWAVE  gains  are  due  to  the  fact  that 
SKYWAVE  uses  the  free- space  value  of  radiation  resistance  in  computing  gain. 


INVERTED -L 


Figures  26,  27,  and  28  show  comparisons  between  APACK,  NEC,  and  ITS 
SKYWAVE  for  an  inverted-L  (mounted  above  soil)  operating  at  frequencies  of  10, 
20,  and  30  MHz,  respectively,  lhe  differences  between  the  predicted  gains  (as 
much  as  7  d  B)  are  believed  to  be  due  mainly  to  differences  in  the  values  used 
for  the  radiation  resistance,  lhe  reason  for  the  high  values  of  gain 
predicted  by  NEC  at  30  MHz,  as  compared  with  the  APACK  and  SKYWAVE  values,  has 
not  been  explained. 

ARBITRARILY  TILTED  DIPOLE 


The  horizontal  dipole  and  elevated  vertical  dipole  are  both  special  cases 
of  the  arbitrarily  oriented  dipole.  Since  comparisons  of  APACK  predictions 
with  other  data  for  both  horizontal  and  vertical  orientations  indicated 
reasonable  agreement,  comparisons  for  other  orientations  were  not  made. 


SLOPING  LONG-WIRE 


Figure  29  shows  comparisons  between  APACK,  NEC,  and  ITS  SKYWAVE 
predictions  for  a  sloping  long-wire  (mounted  over  soil)  operating  at  10  MHz. 
The  predictions  of  APACK  and  SKYWAVE  in  this  case  are  identical,  but  both 
APACK  and  SKYWAVE  predict  lower  gains  than  NEC.  lhe  difference  in  gains, 
varying  from  as  much  as  3  dB  at  low  elevation  angles  to  as  much  as  5  dB  near 
the  zenith  is  believed  to  be  due  to  the  approximations  made  by  APACK  SKYWAVE 
for  the  radiation  resistance. 
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Figure  26.  Elevation  patterns  of  am  inverted-L  mounted 
above  soil  (10  MHz). 


93 


GAIN  IN  dBi 


ESD-TR-80-102 


Section  7 


Figure  27.  Elevation  patterns  of  an  inverted-L  mounted 
above  soil  (20  MHz) . 
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Figure  29.  Elevation  patterns  of  a  sloping  long-wire 
mounted  above  soil. 
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Figure  30  shows  comparisons  between  APACK,  NEC,  and  Ma's  three-term 
current-distribution  model  for  a  terminated  sloping-V  (mounted  over  soil) 


operating  at  a  frequency  of  10  MHz.  Figure  30  gives  gains  as  a  function  of 
elevation  angle,  and  Figure  31  gives  gains  as  a  function  of  azimuth  angle  for 
the  same  antenna  and  operating  frequency.  The  gains  shown  in  Figure  31  are 
for  an  elevation  angle  of  14  degrees. 

Figure  30  indicates  excellent  agreement  in  the  elevation  pattern,  with 
gains  predicted  by  APACX  differing  from  other  data  by  about  1  dB  or  less.  The 
agreement  indicated  in  Figure  31  for  the  azimuth  pattern  is  also  very  good, 
with  gains  predicted  by  APACK  differing  from  other  data  by  3  dB  or  less. 

TERMINATED  SLOPING  RHOMBIC 


Figures  32  and  33  show  comparisons  between  APACX,  NEC,  and  Ma' s 
predictions  for  a  terminated  sloping  rhombic  over  sea  water  for  an  operating 
frequency  of  10  MHz.  Figure  32,  comparing  elevation  patterns,  indicates 
agreement  between  APACK  and  other  data  within  5  dB  at  most,  with  3  dB  being 
typical.  Figure  33,  comparing  azimuth  patterns  at  an  elevation  angle  of  20 
degrees,  indicates  agreement  between  APACK  and  other  data  within  approximately 
3  dB. 

TERMINATED  HORIZONTAL  RHOMBIC 


Figure  34  shows  comparisons  between  APACK,  NEC,  and  SKYVBVVE  for  a 
terminated  horizontal  rhombic  operating  at  10  MHz  mounted  over  soil.  The 
elevation  gains  shown  in  Figure  34  indicate  that  APACK  predictions  differ  from 
other  data  by  not  more  than  2  dB  at  low  elevation  angles  and  by  not  more  than 
5  dB  near  the  zenith. 
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Figure  30.  Elevation  patterns  of  a  terminated  sloping-v 
mounted  above  soil. 


98 


GAIN  IN  dBi 


ESD-TR-B0-102 


Section 


Figure  31.  Azimuth  patterns  of  a  terminated  sloping-v  (mounted  above  soil) 
at  an  elevation  angle  of  14°. 
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Figure  32.  Elevation  patterns  of  a  terminated  sloping  rhombic  mounted  above  sea 
water . 
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Figure  34.  Elevation  patterns  of  a  terminated  horizontal 
rhombic  mounted  above  soil. 
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SIDE-LOADED  VERTICAL  HALF-RHOMBIC 


Figure  35  gives  comparisons  between  APACK,  NEC,  and  SKYWAVE  predictions 
for  the  elevation  patterns  of  a  side-loaded  vertical  half-rhombic  mounted  over 
soil  and  operating  at  a  frequency  of  10  MHz.  m  this  case,  the  gains 
predicted  by  APACK  differ  from  other  data  by  3  dB  or  less. 

Figures  36  and  37  show  elevation  and  azimuth  patterns,  respectively,  of  a 
side- loaded  vertical  half-rhombic  operting  at  10  MHz  mounted  over  sea  water. 
The  comparisons  between  APACK,  NEC,  and  fte' s  model  for  elevation  patterns  (see 
Figure  36)  indicate  agreement  generally  within  3  dB.  Hie  comparisons  for 
azimuth  patterns  at  an  elevation  angle  of  2  degrees  (see  Figure  37)  indicate 
similar  agreement. 

HORIZONTAL  YAGI-UDA  ARRAY 


Figure  38  gives  comparisons  between  APACK,  NEC,  and  Ma' s  predictions  for 
gains  of  a  horizontal  Yagi-Uda  array  operating  at  10  MHz  over  soil  and  sea 
water.  Bccellent  agreement  is  indicated  with  APACK  gains  differing  from  other 
data  by  1  dB  or  less  for  all  elevation  angles  above  10  degrees. 


HORIZONTALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 


Figures  39  and  40  compare  APACK  predictions  with  Ma' s  predictions  for 
gains  of  a  horizontally  polarized  log-periodic  dipole  array  operating  at. 

12  MHz  mounted  over  both  soil  and  sea  water.  The  elevation-pattern 
comparisons  in  Eigure  39  indicate  excellent  agreement  for  soil  (within  1  dB) 
and  good  agreement  for  sea  water  (within  2  dB)  for  all  elevation  angles  above 
10  degrees.  The  azimuth-plane  comparisons  made  at  an  elevation  angle  of  36 
degrees  (see  Figure  40)  indicate  agreement  within  3  dB  for  all  azimuth  angles. 

Figures  41,  42,  and  43  show  comparisons  between  APACK,  NEC,  and 
manufacturer's  data  for  the  Cbllins  237B-3  horizontally  polarized  log-periodic 
dipole  array  (above  soil)  operating  at  frequencies  of  8,  12,  and  20  MHz, 
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Figure  35.  Elevation  patterns  of  a  side-loaded  vertical 
half- rhombic  mounted  above  soil. 
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ELEVATION  ANGLE 

Figure  36.  Elevation  patterns  of  a  side-loaded  vertical 
half- rhombic  mounted  above  sea  water. 
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FREQUENCY:  10  MHz 
SEAWATER:  6r  =  80 
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jure  40.  Azimuth  patterns  of  a  horizontally  polarized  array  (mounted  above  sea  water  and 
soil)  at  an  elevation  angle  of  36°. 


FREQUENCY:  12  MHz 
SOIL:  6r  =30 


Elevation  patterns  of  a  Collins  237B-3  horizontally  polarized  log-periodic  dipole  array 
mounted  above  soil  (12  MHz). 
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respectively.  The  predictions  of  APACK  agree  with  the  manufacturer's  data 
within  5  dB  or  less  in  general,  but  the  predicted  gains  of  the  SKYWAVE  program 
are  much  lower . 

Not  only  are  the  SKYWAVE  gain  values  much  lower  than  those  predicted  by 
APACK  and  the  manufacturer,  but  the  peaks  and  nulls  of  the  SKYVAVE  patterns  do 
not  coincide  with  the  other  predictions.  The  inaccuracy  of  the  SKYWAVE 
predictions  is  believed  to  be  due  not  to  the  theory  on  which  SKYVAVE  is  based, 
but  rather  due  to  coding  errors  in  the  horizontal  log-periodic  routine  that 
have  not  been  corrected  by  ITS. 

VERTICALLY  POLARIZED  LOG-PERIODIC 


Figure  44  compares  elevation- pattern  predictions  by  APACK  and  Ma' s  model 
for  a  vertically  polarized  log-periodic  dipole  array  mounted  over  sea  water  at 
frequencies  of  6  and  18  MHz.  The  comparisons  at  6  MHz  show  agreement  between 
the  two  models  to  less  than  2  dB,  while  the  comparisons  at  18  MHz  show 
agreement  to  within  5  dB. 

CURTAIN  ARRAY 


Figure  45  shows  the  elevation  pattern  predicted  by  APACK  for  a  two-bay, 
four-stack  curtain  array  operating  at  10  MHz  over  soil.  The  APACK  predictions 
shown  in  Figure  45  are  essentially  identical  to  those  obtained  by  Earghausen 
et  al .  (see  Reference  1,  p.  230). 

Figure  46  compares  elevation- plane  patterns  predicted  by  APACK  and  Ma's 
model  for  a  one-bay,  two- stack  curtain  array  operating  at  10  MHz  over  sea 
water.  The  predictions  agree  within  2  dB  for  all  elevation  angles  above 
5  degrees.  The  predictions  are  identical  for  elevation  angles  above  about 
60  degrees . 
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Figure  44.  Elevation  patterns  of  a  vertically  polarized 

log-periodic  dipole  array  mounted  above  sea  water. 
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Figure  45.  Elevation  pattern  predicted  by  APACK  for  a  two-bay, 
four-stack  curtain  array  mounted  above  soil. 
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Figure  46.  Elevation  patterns  of  a  one-bay,  two-stack 
curtain  array  mounted  above  sea  water. 
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47.  Azimuth  pattern  predicted  by  APACX  for  a  sloping  double 
rhomboid  (mounted  over  soil)  at  an  elevation  angle  of 
20°  . 
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SECTION  8 

COMPARISONS  BETWEEN  TRANSMISSION  LOSS  PREDICTED 
BY  APACK  AND  OTHER  DATA 

A  number  of  comparisons  were  made  between  transmission  loss  predicted  by 
APACK  and  transmission  loss  predicted  by  other  data.  In  all  cases, 
comparisons  are  made  with  respect  to  basic  transmission  loss,  i.e.,  the  path 
loss  that  would  result  if  the  two  actual  antennas  were  replaced  by  lossless 
isotropic  antennas . 

The  data  with  which  APACK  transmission- loss  predictions  were  compared  are 
field- strength  curves  given  by  the  CCIR  (see  Reference  26)  and  results 
obtained  from  two  other  ground-wave  propagation  models  used  at  ECAC,  referred 
to  as  the  IPS  model  (see  Reference  3)  and  the  NX  models  (see  Reference  4). 

The  CCIR  curves  were  converted  to  basic  transmission  loss  using  aquation 
6-12.  The  IPS  and  NX  models  predict  basic  transmission  loss  directly. 

The  CCIR  curves  which  are  based  on  Bremmer’ s  formulation  (see  Reference 
17)  assixne  that  both  the  transmitting  and  receiving  antennas  are  located  on 
the  ground  and  that  both  antennas  are  Hertzian  dipoles.  The  IPS  and  NX  models 
account  for  the  heights  of  the  antenna  feed  points  above  ground.  The  APACK 
predictions  of  transmission  loss  account  for  the  actual  antenna  structure. 

The  comparisons  predicted  below  are  made  to  indicate  the  accuracy  and 
versatility  of  APACK  for  predicting  transmission  loss  for  frequencies  between 
150  kHz  and  500  MHz.  The  comparisons  include  ground-wave  propagation  over 
both  average  soil  and  sea  water.  The  differences  indicated  between  the 
various  predictions  shown  in  the  figures  are  not  meant  to  represent  either  the 
"best"  or  "worst"  case.  Rather,  the  differences  should  be  considered  typical. 

Figures  48,  49,  and  50  show  the  basic  transmission  loss  predicted  by 
APACK,  CCIR  curves,  and  the  IPS  model  at  1,  3,  and  10  MHz,  respectively  over 
soil  for  vertical  polarization.  Typical  differences  between  losses  predicted 
by  APACK  and  the  other  models  are  5  dB  or  less . 


119 


FREQUENCY = 
SOIL  : 


Figure  48.  Comparisons  between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground- 
wave  propagation  over  soil  at  1  MHz  (vertical  polarization) . 
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Figure  49.  Comparisons  between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground 
wave  propagation  over  soil  at  3  MHz  (vertical  polarization) . 


FREQUENCY :  10MHz 
SOIL  :  Cr=  4 


Figure  50.  Comparisons  between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground 
wave  propagation  over  soil  at  10  MHz  (vertical  polarization ) . 
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Figures  51,  52,  and  53  show  the  basic  transmission  loss  predicted  by 
APACK,  CCIR  curves,  and  the  IPS  model  at  1,  3,  and  10  MHz,  respectively,  over 
sea  water  for  vertical  polarization.  Typical  differences  between  losses 
predicted  by  APACK  and  the  other  models  are  again  5  dB  or  less. 

Figures  54  and  55  compare  APACK  predictions  with  the  CCIR  curves  at 
150  kHz  and  10  MHz,  respectively,  for  soil  with  slightly  different  constants 
than  used  in  previous  comparisons.  Vertical  polarization  is  assumed.  The 
differences  at  150  kHz,  shown  in  Figure  54,  are  10  dB  or  less.  The 
differences  at  10  MHz,  shown  in  Figure  55,  are  5  dB  or  less. 

Figures  56  and  57  compare  APACK  predictions  with  those  of  the  NA  model  at 
frequencies  of  42.9  and  100  MHz,  respectively,  for  vertically  polarized 
propagation  over  soil.  Rsr  these  comparisons,  the  transmitting  antenna  is  an 
elevated  vertical  dipole  with  its  feed  point  located  10  meters  above  ground, 
and  the  height  of  the  feed  point  of  the  receiving  antenna  is  as  indicated.  In 
general,  the  APACK  predictions  are  within  2  dB  of  the  NX  predictions.  The 
reason  for  the  large  discrepancies  at  distances  exceeding  approximately  100 
kilometers  is  that  NX  includes  a  model  for  tropospheric-scatter  predictions, 
whereas  APACK  does  not. 

Figure  58  compares  APACK  and  NX  predictions  at  500  MHz  with  both  antennas 
located  on  the  ground  for  propagation  over  average  soil.  APACK  values  are 
generally  within  2-3  dB  of  the  NX  values  except  at  distances  above 
10  kilometers  where  troposcatter  propagation  is  dominant.  In  Figure  58,  the 
term  ”Yeh"  refers  to  the  tropospheric- scatter  model  included  in  NX.  Wien  NX 
predictions  are  made  without  accounting  for  troposcatter,  the  results  are 
essentially  identical  to  those  obtained  with  APACK  for  large  distances. 

Figures  59  and  60  compare  APACK  and  NX  predictions  for  vertically  and 
horizontally  polarized  propagation,  respectively,  over  sea  water  at  a 
frequency  of  2  MHz.  Bbr  Figure  59,  the  transmitting  antenna  is  an  elevated 
vertical  dipole  70  meters  long  with  its  feed  point  located  35  meters  above  the 
surface,  and  the  height  of  the  receiving  antenna  feed  point  is  varied.  Rir 


FREQUENCY :  I  MHz 
SEA  WATER  :  6r=  80 


between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground-wave 
over  sea  water  at  1  MHz  (vertical  polarization). 


FREQUENCY :  3  MHz 
WATER:  Cr=  80 


Comparisons  between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground -wave 
propagation  over  sea  water  at  3  MHz  (vertical  polarization). 


FREQUENCY  =  10  MHz 
SEAWATER  :  6r  =  8C 


53.  Comparisons  between  basic  transmission  loss  predicted  by  APACK,  CCIR,  and  IPS  for  ground-wave 
propagation  over  sea  water  at  10  MHz  (vertical  polarization). 


FREQUENCY  10  MHz 
SOIL  :  6r  =  4 


Comparisons  between  basic  transmission  loss  predicted  by  APACK  and  CCIR  for  ground-wave 
propagation  over  soil  at  10  MHz  (vertical  polarization) . 
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FREQUENCY  :  500  MHz 
SOIL  :  6r  =  4 


Comparisons  between  basic  transmission  loss  predicted  by  APACK  and  NX  for  ground-wave 
propagation  over  soil  at  500  MHz  (vertical  polarization) . 


35m, 300m  FREQUENCY  •'  2  MHz 

SEAWATER:  €rs80 


Comparisons  between  basic  transmission  loss  predicted  by  APACK  and  NX  for  ground-wave 
propagation  over  sea  water  at  2  MHz  (vertical  polarization) . 
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Figure  60,  the  transmitting  antenna  is  a  horizontal  dipole  70  meters  long 
located  35  meters  above  the  surface,  and  the  height  of  the  receiving  antenna 
feed  point  is  varied.  By  canparing  figures  59  and  60,  it  is  seen  that  the 
horizontally  polarized  wave  is  attenuated  more  rapidly  than  the  vertically 
polarized  wave. 

figure  61  is  included  to  show  direct  comparisons  of  field  strengths 
produced  by  different  transmitting  antennas  with  the  same  feed-point  height 
(10  meters)  and  polarization  at  a  fixed  frequency  and  receiving-antenna  feed- 
point  height  (200  meters) .  The  three  transmitting  antennas  compared  are  the 
Hertzian  dipole  used  by  Bremmer  (see  Reference  17),  a  short  vertical  dipole 
(1.6  meters  long),  and  a  half-wave  vertical  dipole. 

As  shown  in  figure  61,  the  field  strengths  produced  by  the  Hertzian 
dipole  are  lower  than  the  fields  produced  by  the  short  dipole.  The 
significant  difference  between  the  fields  of  the  short  dipole  and  the  fields 
of  the  half-wave  dipole  (both  predicted  by  APACK)  is  due  to  the  fact  that 
APACK  considers  the  actual  antenna  structure. 


RMS  FIELD  STRENGTH  (JAV/m) 
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Figure  61.  Comparisons  between  rms  field  strength  predicted  by  APACX 
and  Bremmer  for  ground-wave  propagation  over  soil  at 
42.9  MHz  (vertical  polarization.) 
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SECTION  9 
RESULTS 

Equations  have  been  developed  to  predict  electric  far- field  strengths, 
directive  gains,  power  gains,  and  transmission  losses  for  16  types  of  linear 
antennas  used  at  frequencies  between  150  kHz  and  500  MHz.  The  antennas  for 
which  equations  were  developed  are:  horizontal  dipole,  vertical  monopole, 
vertical  monopole  with  radial-wire  ground  screen,  elevated  vertical  dipole, 
inverted-L,  arbitrarily  tilted  dipole,  sloping  long-wire,  terminated 
sloping-V,  terminated  sloping  rhombic,  terminated  horizontal  rhombic,  side- 
loaded  vertical  half-rhombic,  horizontal  Yagi-Uda  array,  horizontally 
polarized  log-periodic  dipole  array,  vertically  polarized  log-periodic  dipole 
array,  curtain  array,  and  sloping  double  rhomboid  (see  Figures  3  through  18). 

The  equations  include  the  contributions  of  the  surface  wave  in  addition 
to  the  contributions  of  the  direct  and  ground-reflected  waves  so  that  gain 
predictions  for  ground-wave  as  well  a3  sky-wave  analyses  can  be  made 
directly.  The  inclusion  of  the  surface  wave  also  allows  predictions  of 
ground-wave  transmission  loss  accounting  for  the  actual  antenna  structures. 

The  basic  formulation  of  the  APACK  equations  is  for  current  elements 
located  over  lossy  planar  earth  ( see  Figure  1 ) .  This  formulation  was  extended 
to  consider  spherical  earth  within  the  radio  horizon  by  making  use  of  the 
divergence  factor  (see  Figure  2).  In  the  diffraction  region  beyond  the  radio 
horizon,  the  basic  formulation  was  modified  by  using  the  radiation  vector  and 
the  Bremmer  secondary  factor.  APACK  includes  criteria  to  determine  in  which 
region  the  far-field  observation  point  lies  so  that  appropriate  equations  are 
used. 

Cbmparisons  have  been  made  between  APACK  predictions  of  antenna  gain  and 
other  available  data  (see  Figures  19  through  47).  Reasonable  agreement  is 
indicated. 
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Comparisons  have  also  been  made  between  APACK  predictions  of  transmission 
loss  and  other  available  data  (see  Figures  48  through  61).  Reasonable 
agreement  is  again  indicated. 

The  APACK  equations  documented  in  this  report  have  been  automated  at  ECAC 
to  provide  a  rapidly  executing  model  for  predicting  gains  and  transmission 
loss  for  the  16  antennas  considered.  APACK  allows  the  antennas  to  be 
described  in  terms  of  simple  geometrical  parameters,  as  opposed  to  the 
detailed  structure  description  required  by  rigorous  method- of- moments 
programs,  to  simplify  the  use  of  the  model  by  EMC  engineers. 

The  analytic  expressions  for  the  direct  and  ground-reflected  waves  have 
been  derived  by  others.  However,  the  analytic  expressions  for  the  surface 
wave  represent  new  contributions  to  the  antenna  and  propagation  literature. 
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APPENDIX  A 


FORMULAS  FOR  CALCULATING  FIELD  STRENGTH  AND  GAIN  OVER  PLANAR 
EARTH  AND  IN  THE  DIFFRACTION  REGION  BEYOND  THE  RADIO  HORIZON 

HORIZONTAL  DIPLOLE 


The  field  intensities  for  a  horizontal  dipole  are  derived  from  those  for 

a  Hertzian  dipole  by  using  the  principle  of  superposition.  Therefore, 

consider  a  Hertzian  dipole  located  H  above  planar  earth  with  its  axis  in  the 

same  direction  as  the  X-axis.  The  electric  field  intensities  in  cylindrical 

27 

coordinates  (p,  <)> ,  Z)  are  given  as  follows: 
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27King,  R.  W.P.,  The  Theory  of  Linear  Antennas,  Harvard  University  Press, 
Cambridge,  MA,  1956. 
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Ep  and  Ez  can  be  converted  into  the  component  Eg  in  spherical  coordinates  (r, 
9,0)  as: 


E.  =  E  cos  9  -  E  sin  9 
9  p  z 


(A-4) 
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Assuming  9d  =  9f  =  9 : 
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where 


- -  =  30  k,  e,  =  e 
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In  the  far  field,  let: 
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Rsr  the  horizontal  dipole  : 
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and 


cos  =  sin  0  cos  $ 


(A- 10) 


Assume  that  the  current  distribution  on  the  antenna  is  given  byi 


I  = 


I  sin  k  (l  -  x) ,  x  >  0 
m 
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and  that  the  length  of  antenna  is  2i .  Then; 
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where  e^1**  cos  ^  is  the  phase  advance  of  the  current  element  dx  located  at  x. 
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Equations  A-15  and  A-16  are  identical  to  Equations  65  and  66  of  Reference  8  if 
the  surface-wave  terms  are  neglected.  Rsr  ground-wave  calculations,  these 
equations  may  be  written  as  follows  : 
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(  A- 1 7 ) 
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( A- 18) 


where 


Rr  = 


Jcr  +  (h  -  H)z 

/d2  +  (h  +  H)2 

horizontal  distance  from  the  transmitter  to  the  receiver 
height  of  the  receiving  antenna 
height  of  the  transmitting  antenna. 


Equations  A-17  and  A-18  are  identical  to  Equations  10,  11,  and  12  of  Reference 
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Hie  mutual  impedance,  referred  to  the  dipole  base,  can  be  calculated  from 
Reference  8  as : 
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The  self- impedance,  obtained  from  Equation  A-19  if  d 
radius  of  the  dipole ,  is  given  by : 


=  2a ,  where  a  is  the 
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(A-21) 
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Equation  A-19  is  not  valid  when  l  is  an  integral  multiple  of  a  half¬ 
wavelength.  In  these  cases,  the  impedances  should  be  calculated  based  on  the 
assumption  of  a  more  accurate  form  of  current  distribution  (see  Reference 
9 ) .  The  total  input  resistance  is : 


R. 

m 


Re  (Z,.  +  R.  '  Z  ) 

11  h  m 


R^  is  given  by  Equation  2-19. 


( A-22) 


The  directive  gain  can  be  found  using  Equations  A-1A,  A-16,  and  A-22  as: 
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The  radiation  vectors  are  given  by : 
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The  free-space  field  intensities  are  given  by : 
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Then  : 
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where  E0  and  E^  are  given  by  Equations  A-15  and  A-16.  The  attenuation 
relative  to  free  space  is : 
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The  electric  field  intensities  in  the  diffraction  region  are  given  by: 
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VERTICAL  MONOPOLE 


The  field  intensities  for  the  vertical  monopole  are  also  derived  from 
those  for  a  Hertzian  dipole.  Hie  vertical  and  radial  components  of  electric 
frield  intensity  due  to  a  vertical  Hertzian  dipole  ( see  Reference  27 )  are  given 
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Due  to  symmetry.  Ex  =  0. 


In  spherical  coordinates  ( r ,  9 ,  $ ) 
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Again,  the  following  approximations  are  made: 
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Now  assune  that  the  current  distribution  is  given  by : 


I  sin  k  (£- z) 
tn 


.-36) 


-37) 


149 


ESD-TR-80-1 02 


Appendix  A 


Then : 
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where 


,  rl  jkz  COS  0  ,  ,  , 

J  =  e  Sln  k  (*-2)  132 

( A-4 1 ) 

=  [cos  (kJt  cos  9)  -  cos  k&]  +  j  [sin  (k£  cos  9)  -  cos  0  sinki] 

k  sin2  0 


rl  -jkz  COS  0  ...  ,  . 

J  e  J  sin  k  {l- z)  dz 


=  [cos  (k l  cos  9)  -  cos  k l]  -  j  [sin  (kft  cos  9)  -  cos  9  sinkft] 

2 

k  sin  0 


( A-42 ) 


rl  -jkz  COS  0 
Jn  e 


-P 


sin  k  (A- z)  j  J rrPg  e  erfc  (j  ^P^ >  dz  (A-43) 


Equation  A-43  is  evaluated  in  APPENDIX  B  as: 
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Then  : 
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Simplifying  further  : 
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where 

Aj  =  cos  (kZ  cos  0)  -  cos  k Z  (A-47) 

B2  =  sin  (  kZ  cos  6)  -  cos  0  sin  k 2  (A-48) 

Hie  first  two  terms  of  Equation  A-46  (direct  and  ground- reflected  wave 
contributions)  are  identical  to  Equation  29  in  Reference  8.  The  last  two 
terms  represent  the  surface-wave  contribution  to  the  field. 

The  radiation  resistance  of  a  vertical  monopole  over  lossy  earth  is  very 
difficult  to  calculate.  In  fact,  no  accurate  representation  is  available  in 
the  literature . 

2  8  29 

Hie  method  of  Sommerfeld  and  Renner  '  cannot  be  used  here  because  it 

is  not  applicable  when  the  dipole  touches  the  ground.  Thus,  the  expression 
for  a  perfectly  conducting  ground  will  be  used  here.  Hie  input  resistance  for 
perfectly  conducting  ground  is : 


30 


in  .  2  . 

sin  kZ 


(HCin  2 


k 2)  Cin  (2  kZ)  -  —  cos  2<2  Gn  (4k2) 


(  A-49 ) 


+  -  sin  2  kZ  Si  (  4k  Z )  -  sin  2k2  Si  (2k2) 


2®Sommerfeld,  A.  and  Renner,  F.,  "Strahlungsenergie  und  Erdabsorption  bei 
Dipolantennen , "  Ann.  Physik,  Vol.  41,  pp.  1-36,  1942. 

2®Kuebler,  W.  and  Karwath,  A.,  Program  RENNER;  Normalized  Ifediation 
Resistance  of  a  Hertzian  Dipole  over  Arbitrary  Gkound,  ECAC-TN-75-024, 
Electromagnetic  Compatibility  Analysis  Center,  Annapolis,  MD,  January  1976. 
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where 

CLn( x)  =  0.577  +  In  x  -  Ci  (x) 

,  ,  ex  sin  u 

Si(  x)  =  Jn  -  du 

•'0  u 

.  rx  cos  u  , 

Ci(  x)  =  I  -  du 

J  at  u 

(0.577  is  Euler’s  constant) 

The  directive  gain  of  the  dipole  is  given  by: 


30  I2  sin2  (ki)  R , 
m  in 


( A-50 ) 


Since  gp  =  n  g^,  the  power  gain  of  the  dipole  is  given  by: 


n  r2|E9|2 

30  I2  sin2  (kA)  R, 
m  in 


( A-51 ) 


where 


The  radiation  efficiency  is  given  in  Reference  14  as: 


( A-52 ) 


ESD-TR-8  0-102 


n(dB) 


6416.702 
+  364.817  (yj 


+  6091.33 

-  25.646 


2179.89 


The  radiation  vector  is  obtained  by: 


„  .  jkz  cos 

=  -  I.  I  sin  k  ( l-z )  e 

'fl  m 


dz 


jm  {cos  (k&  cos  6)  -  cos  kfc}  +  j  (sin  (k&  cos  8)  - 


•  2  Q 

sin  6 


N.  =  N  sin  0  = 
o  z 


5.  A2  *  j  B2 

k  sin  0 


The  free-space  field  intensity  is  given  by  : 


j  30  I  e_^kr  A  +  j  B„ 
_ _ m _  2  J  2 

0f  r  sin  0 


Hie  attenuation  relative  to  free  space  is  : 
|E0| 

A  =  >V 

ADB  *  20  log  j^j- 


The  electric  field  intensity  in  the  diffraction  region  is  given  by 


Appendix  A 

(A-53) 


(A-54 ) 

cos  0  sin  k£. } 

(A-55) 


(A-56 ) 

(A-57) 

(A-58) 
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j30  I  e 
m 


■jkd 


sin  6 


2  F 


( A-59  ) 


VERTICAL  MONOPOLE  WITH  RADIAL-WIRE  GROUND  SCREEN 


Most  permanent  monopole  installations  include  a  ground  screen,  ihe 

30 

presence  of  a  ground  screen  has  been  studied  by  Wait  and  Pope,  Wait  and 
SUrtees,31  wait,32  and  teiley  and  King.33 

The  presence  of  a  ground  screen  will  modify  the  field  in  the  absence  of 
the  screen  by  a  correction  factor  A  Eg  given  by: 

Eg  =  Eg  +  A  Eg  (A-60) 

when  Eg  (the  field  without  the  screen)  is  given  by: 


30Wait,  J.R.  and  Pape,  W.A.,  "The  Characteristics  of  a  Vsrtical  Antenna  with 
a  Radial  Conductor  Ground  System,"  Appl .  Sci .  Res . ,  Section  B.4,  1954 
pp.  177-185. 

31Wait,  J.R.  and  Surtees,  W.J. ,  "impedance  of  a  Top- loaded  Antenna  of 

Arbitrary  Length  Over  a  Circular  Grounded  Screen,"  J.  Appl.  Physics,  Vol.  25, 
1954,  pp.  553-555.  " 

32Wait,  J.R.,  "Effect  of  the  Ground  Screen  on  the  Field  Radiated  from  a 
Monopole,"  IRE  Trans.  Antennas  and  Propagation,  Vol.  AP-4,  No.  2,  1956, 
pp.  179-181. 

33Maley,  S.W.  and  King,  R. J. ,  "Ihe  impedance  of  a  Monopole  Antenna  With  a 

Circular  Oonducting-Disk  Ground  System  on  the  Surface  of  a  Lossy  Half-space,” 
J.  Research  NBS,  65D  (Radio  Propagation)  No.  2,  1961,  pp.  183-188. 
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E-  = 


j30  I  e 
m 


"jkr 


A,  +  jB  A  -jB 

1  2  +  R  ■ 


sin  9 


v  sin  9 


5  nr 

+  H-Rv)  (sin  9  -  v— - 

n 


sin2  9 


cos  9 ) 


(A-61 ) 


A2  ^b2 
sin  9 


I -  ~Pe  2 

{ 1  —  j  v it  e  sin  9  erfc  (j  >/p^ ) } 

Jr\2  - 


-j2  F  (sin  9  - 
e 


— r-  — ~~2 — 

sin  9  „ .  vn  -  sin  0 

— -  cos  0  )  - - -  sin  kJl 

n  sin  0 


where 


A2  =  cos  (kJl  cos  0)  -  cos  k£ 


32  »  sin  (ki  cos  9)  -  cos  9  sin  k£ 


(A-62) 

( A-63 ) 


AE„ 


n  sin  9  /' 


ka 


*2*2 


-  e 


jx 


cos  k£]  J^xsinB)  dx 


120  tt  sin  k A  [cos  ( ki  cos  9)  -  cos  k£] 


(A-64) 


3  . 


n  sin  9  /^a  [e 


-j  Vx  + 


2  2 

k  ^  -e  cos  k£]  J1  (x  sin  0)  dx 


120ir  sin  k£  [cos  (k£  cos  9)  -  cos  k£] 


(A-65) 


where 

•  Bessel  function  of  the  first  kind 
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ju>Uo 

n  =  - : -  ,  the  characteristic  impedance  of  the  ground 

o  +  -joje  e 
or 

UQ  =  4xx10  ,  the  permeability  of  free  space 

—  1  2 

eQ  =  8.854x10  ,  the  permittivity  of  free  space 

er  =  relative  dielectric  constant  of  the  ground 
a  -  conductivity  of  the  ground 
a)  =  cuigular  frequency. 


(free- space  value) 

(perfectly  conducting  ground) 


Wiit  and  R>pe  (Reference  30)  have  shown  that  the  effect  on  the  radiation 
resistance  caused  by  the  ground  screen  can  be  expressed  as  : 


R '  *  R.  +  Re  (AZ,  +  AZ„)  ( A-66 ) 

in  m  12 


where 


R.  =  - - -  f  (1  +  cos  2  ki.)  Cin  (2  kJl )  -  —  cos  2  k£  dn  (4  kJl ) 

in  .  2  .  .  I  2 

sin  M  * 

-  sin  2  W  Si  (2  kA)  +  -jj  sin  2  kJl  Si  (4  kX )  ] 


(A-67) 
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and  A corresponds  to  the  self- impedance  of  the  monopole  over  a  perfectly 
conducting  discoid.  LZy  can  be  expressed  as: 


S  -  - - -  \  e  Ci  (2k  (r  +  £))  +  j  (-j  -  Si  (2k  (rQ  +  »))) 

4ir  sin  k£  1  ° 

+  e“j2kA  [d  (2k  (rQ-£))  +  j  (j  -  Si  (2k  (rQ-£)))j 


+  2  cos2  k£  CL  (2  ka)  +  j  (j  -  Si(2ka)) 


+  4  cos  k£ 


[a.  ik^ 


)  +  j  (-  -  Si  (kr 


i”  I 


(A-68) 


-  4  cos  k£  e"jkA  [a.  (k(  -  £))  +  j  (|  -  Si  (k( ?,-£))) 


4  cos  k£  eJ 


CL  (k(r1 


+  *)>  +  j  (j  -  Si(  k(  +  £)))j 


where 


/a2  +  £2 


r1  -  a  +  \/a2  +  t2 


( A-69 ) 


(A-70) 


AZj  accounts  for  the  finite  surface  impedance  of  the  radial- conductor  system 
and  is  given  by : 
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_  f  ^  n  Ce~jk  Vp2  »  -e~jkp  cos  k*]2 
0  +  11  2  n  p  sin2  kJL 


( A-71 ) 


where 


240  it  p 

i  NX  ln  NC 


C  =  radius  of  the  conducting  screen 


N  =  number  of  radial  conductors 


(A-72) 


The  directive  gain  for  the  monopole  over  a  radial-wire  ground  screen  becomes : 


2  .  '  i  2 

r  ’V 

30  l2  sin2(k£)  R. 
m  m 


(The  radiation  efficiency  is  still  given  by  Equation  A-53.) 


The  free-space  field  intensity  is  given  by: 


j30  1  e_:)kr  A  +  j  B„ 
J  m  2  2 


The  attenuation  relative  to  free  space  is 


(A-73) 


( A-74 ) 


(A-75) 
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ADB 


20  1o’  W 


(  A -76  ) 


The  electric  field  intensity  in  the  diffraction  region  is  given  by: 


j3°  I  e 

HI 


-  j  kd 


A2  +  j  B2 

sin  6 


2  F 


(  A -77  ) 


ELEVATED  VERTICAL  DIPOLE 


The  current  distribution  assumed  for  the  elevated  vertical  dipole  is: 


34 


I 


I  sin  k  (A  +  z  -  z)  , 
m  o 


z  <  z  <  z  +  A 

o  —  -  o 


(  A-78) 


I  sin  k  (A  -  z  +  z)  , 
in  o 


z  -  A  <  z  <  z 
o  —  —  o 


(A-79) 


The  coordinate  system  is  identical  to  that  considered  for  the  vertical 
monopole.  From  Equation  A-38: 


34 

Jordan,  E.C.  and  Balmain,  K.G. , 
Prentice  Hall,  Englewood  Cliffs 


Electromagnetic  Waves  and  Radiating  Systems, 
,  NJ,  1968. 
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-jkr 


E  =  j30  kid*  — — 

o  tn  r 


sin  6 


jkz  cos9 


2  2 

.  _  -jkz  cos0  .  .  _  _  ,  _  -jkz  cos0  ,  .  2a  V n  -sin  0  „\ 

+  R  e  J  +  (1-R  )  F  e  J  (sm  0  -  - - -  cos  0  I 

v  v  e  2  / 


WLth  the  current  distribution  as  given  by  Equations  A-78  and  A-79 


“jkr  Z 

E_  =  j30  kl  -  sin  0  /  .  sin  k(£-Z  +  z) 

0  m  r  J  Z  -l  p 

o 


jkz  cos0 


12  2 

_  -jkz  cos0  2  „  Vn  -  sin  0 

+  R  e  +  ( 1  -R  )  ( sm  8  -  - - -  cos 


0)  Fe-3^cos0 


Z  +*. 


/  °  sin  k  (4  +  Zq-z) 
o 

l~T—T, ; 

vn-f 


jkz  cos  0  ,  _  -jkz  cos 

>  +■  R  Q 


.......  2n  Vn  -sm  0  _  -Tkz  cos  0 

+  ( 1  -R  )  ( sm  0  -  - - — -  cos  0  )  F  e  J 

v  2  e 

n 


dz 


( A-80 ) 


Equation  A-81  is  rewritten  as : 


E  =  j30  kl 
a  m 


-jkr 

— —  sin  0  +  Rv  (J2-W2) 


+  ( 1 -R  )  (sin20  - 


\/n2  -  sin20 


cos  0)  (J3  +  J3)] 


(A-81) 
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where 


o  jkz  cos0  ,  ,  „  .  , 

z  e  sin  k(A-z  +  z)  dz 

o 


/ 

Jz  - 

3 

jk  (Z  -A)  cos0  jkZ  cos  0 


-e 


(cos  kA  -j  cos  0  sin  kA) 


k  sin  0 


(  A-82 ) 


J-  = 


Z  +J1 

C  o  jkz  cos©  ,  , 

e  sin  k  (A+z  -z)  dz 

Zo  ° 

jk(  Z  +A)  cos©  jkz  cos0 
o  o 

_  e _ 2® _ ( cos  kA  +  jcosQ  sin  kA ) 

.  2  „ 
k  sin  0 


(  A-83 ) 


J1  +  J1 


=  2  e 


jkz  cos0 


cos  (kA  cos©)  -  cos  kA 
k  sin^  0 


( A-84 ) 


I  ° 

JZ  -A 
o 


-jkz  cos© 

e  J  8  sin  k  (A-Z  +  z)  dz 

o 

-jkZ  COS0  , _ ^  Q  (  A-85 ) 


- jk  ( Z  -A)  cos© 

o  -e 


(cos  kA  +  jcos  0  sin  kA) 


k  sin  0 


Z  +A  -jkz  cos@ 
r  o  o 

e  sin  k(  A+Z  -z )  dz 

Zo  ° 

-jk(Z  +A)  cos0  -jkZ  cos0  ,  „ 

o  J  o  ( cos  c  A 

e _  -e 


(  A-86 ) 


j  cos©  sin  kA) 


k  sin  0 
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+  Jo 
2  2 


=  2  e 


-jkz  cosQ  ,,  ,  „ 

o  cos  (k l  cosQ )  -  cos  k& 


(A-87) 


k  sin  8 


/z°_£  e  cos®  sin  k  (4-Zq+z)  [  1  —  j  e  6  erfctj^P^)]  ds 

o 


( A-88 ) 


/  °  e  jkzcosB  gin  ^^£+2  _z)  [1- j  ^/ttP  e  Pe  erfc  (j  /p~)]dz 

^  OS  »  6 

°  (A-89) 


J3  +  J3 


_  -ikz  cos0  cos  (k£  cos0)  -  cos  kZ 
2  e  o  - - - 

k  sin  0 


2e 


-jkz  cos@ 
o 


1  -  R 


\l~ 2  " 

V n  ~  SI 


•  2n 
_sxn  0 

2  ' 
n  sin0 


x  [cos  ( k£  cos0)  -  cos  k£]  e  e  erfc  ( j  \/^^— 


( A-90 ) 


Substituting  Equations  A-84,  A-8  7,  and  A-90  into  Equation  A-81  : 
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_  e  cos  (k£  cos0 )  -  cos  k l 

Ee  =  360Im“ - sine - 


tjkZ  cos0 
° 


-jkZ  cos  0  r 2 

+  R  e  °  +  (1-R  )  <sin20  -  - - — 

v  v  2 

n 


(  A-91  ) 


cos  0 ) 


x  e 


-jkZ  cos0  „  -P 

O  ,  .  ^  _  e  .  ..  /r)} 


f  / -  2  6 

{ 1  —  j  J  irPg  sin  0  e  erfc  (j 


=  j60  I 


-jkr  jkZ  cose 
e  e  o 

m  r 


cos  (k i  cos0 )  -  cos  k l 
sin  0 


( A  -92 ) 


Equation  A-91  has  been  derived  previously  in  the  literature  for  the  case  ZQ  - 
0.  The  second  and  third  terms  of  Equation  A-91  represent  the  ground- 

reflected  wave  and  surface  wave,  respectively. 

If: 

R,  =  r  -  Z  cos  0 

d  °  (  A-93 ) 

R  =  r  +  Z  cos  9 
r  o 


Equation  A-91  may  be  written  as: 


^^Ramo,  S. ,  Whinnery,  J.  R. ,  and  Van  Duzer,  T. ,  Fields  and  Vfaves  in 
Communication  Electronics,  John  Wiley  and  Sons,  New  York,  NY,  1965. 
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E0  =  *0Im 


cos  ( k&  cos9)  -  cos  k& 
sin6 


-jkR 


-jkR 


+  R 


V  R 


+  ( 1  -R  )  ( sin  9 

v 


f 2  ~T7 

yn  -  sin  9 


(  A -94  ) 


cos  0  ) 


•jkR. 


2  -pe 

{ 1  -  j  sin  9  e  e  erfc  (j  ^P^)} 


For  calculation  of  the  ground-wave  field  intensity.  Equation  A-94  may  be  used 

for  better  accuracy.  It  can  be  obtained  frcrn  Equations  1  and  2  in  Reference  20 

2 

if  9r  =  0^  =  0  and  sin  0=1  inside  the  braces. 

An  analysis  of  the  impedance  of  am  antenna  of  length  2 l  and  a  small  radius  a 
at  a  height  zQ  >  t  over  a  lossy  ground  is  an  intricate  problem  and  has  not 
been  solved  completely  (Reference  27)  .  An  approximate  solution  is  given  by 
Sommerfield  amd  Renner  (Reference  28).  Et>r  a  Hertizan  dipole  or  for  a  very 
short  end- loaded  amtenna  of  length  21  with  uniform  current  distribution, 
Sommerfeld  amd  Renner  obtain: 


=  120  (  kJl ) 


+  2'i 


„  sin  2  kZ  -  2  kZ  cos  2  kZ 

|  ♦  2  i - 2 - _a - a , 

( 2kz  ) 

o 

cos  (2kZ  -  y)  +  2kZ  sin  ( ZkZ  -  y) 
_ o _ o _ o _ 

( 2k Z  )2 

o 


(  A-95 ) 


I  [' 


-  cos  Y  Ci  (ZkZ  )  -  sin  y  { Si  (2kZ  )  -  ~  } ] 

O  1  O  2  1 
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where 


k_  k„ 

—  =  |-2  |  e"3Y 
k  'k  1 


( A-96 ) 


2  2 

Bguation  A-95  is  multiplied  by  the  factor  (— )  for  sinusoidal  current 
distribution  as  discussed  in  Reference  28  to  give  : 


R  =  (^)  2  RV 

r  tt  r 


*lhe  directive  gain  is  given  by : 


( A-97 ) 


2  •  .2 

r  E9 

J  j  n  ^ 

30  I  sin  (  k£)  R 
m  r 


( A-98 ) 


the  radiation  vector  is  given  by : 


-2  I  cos  ( k£  cos  8)  -  cos  k£  jkZ  cosS 
m  o 

N  =  - : —  - - -  e 

z  k  ■  2  a 

sin  0 


( A-99 ) 


2  1  , .  o  n  \  jkZ  cos6 

m  cos  (k£  cos8 )  -  cos  ki  o 

Nq  _  _  - -  e 


(  A- 1 00 ) 


Hie  free-space  field  intensity  is  given  by  : 


j60  I  e 
m 


cos  (kft  cos8 )  -  cos  k£ 


( A- 1 0 1  ) 
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the  attenuation  relative  to  free-space  is : 


A  = 


( A- 102) 


and 


ADB 


20 


( A-103) 


the  electric  field  intensity  in  the  diffraction  region  is  given  by: 


•fin  t 

_  3  me  cos  (k&  cosO?  -  cos  k£  . 

E0  d  sin  0  r 


(A- 104 ) 


INVERTED -L 


The  inverted-L  is  analyzed  by  considering  the  vertical  and  horizontal 
sections  separately  and  then  combining  the  results . 

R>r  the  vertical  portion : 
a'  =  tt/2 
r'  =  0 

Hs  =  z 

3  =  Z 

COS  v  =  cos  ® 


I  (S)  =  I  sin  k  (H+fc-z)  z  >  0  (A-105) 

v  m 
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cos  4*. 
n 

^  {s) 


sin  3  sin  <j> 

I  sin  k  (Z-y) ,  y  >  0 


( A-106  ) 

(A-107) 


Following  the  approach  used  in  the  previous  subsections: 


-jkr 


E0,v  "  j3° 


m  r 


A4  +  jB4 


+  (1  -  R  )  (sin  0  - 

v 


+  R 


A4  “  jB4 


sin  0  v  sin  0 


f~2  ~2 

Vn _ -  sm  0 


cos  0 ) 


A4  -  jB4 
sin  0 


-P 


X  {l  -  j  y/nPe  e  e  Sin2  0  erfc  ( j  ^P^)} 


2 

-  j  (sin  0  - 


Jn2  -  sin2  0 


cos  0  ) 


JJT  ■ 2-  2  F 


sin  0 


sin0 


( A- 108 ) 


{sin  k  (H  +1)  -  sin  kZ  cos  (kH  cos  0)  +  j  sin  kl  sin  (kH  cos  0)} 


Rjr  the  horizontal  portion  : 


J8,h 


.  ‘3oi  ainij)  cos8 

m  .  2  . 
sin  , 
h 


■jkr 


A  2  +  B  2  ib"  ri-R  e“j2khcos6 
A5  B5  ejb  [  v  (A- 109) 


rr 

vn  -s 


...  sin20  ,  .  2„  •/n“-sin‘0  „ ,  -  j2kHcos0 

+  ( 1  -R  )  F  x  — - -  (sin  0  -  - r -  cos0  )  e 

v  e  2  .  2 

n  cos0  n 


>,h 


1301  ^  ,/a  2  l  B  2  .jb” 

m  r  sm2*;  V  5  5 


,  .  _  -j2kH  COS0  , .  _  ,  _  -j2kHcos0 

x  1  +  \  6  +  (1“V  Fm  e 


( A- 1 10 ) 
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where 


A,  =  cos  (kH  cos9)  coskJi  -  sin  (kH  cos0)  cos0  sin  k l  -  cos  (k  (H +£)) 

4 


B,  =  cos(kHcos0)  cos0sink£  +  sin(kHcos0)  cosk£  -  cos0  sin  k(H+£)  (A- 

4 


A,  =  cos  ( kJl  cos<|;.  )  -  cos  k£ 
5  n 


(A- 


B,.  =  sin  (kA  cosi);,  )  -  cos  ij>.  sin  kJl 

5  .  h  h 


(  A- 


b"  *  tan 


-1  B5 


(A- 


E0  E0,v  +  E0,h 


(A- 


E<l>  Vh* 


(  A- 


The  radiation  resistance  is  given  by: 

Rin  =  Rinv  +  Rinh 


where  H^nv  is  given  by  Equation  A-49  except  that  l  is  replaced  by  (H  + 
Rinh 


( A-111) 


112) 


113) 


1 14) 


115) 


1 16) 


117) 


is  given  by  : 


ESD-TR-80-102 


R.  . 
inh 


=  60 


('• 


415  +  In  (i|)  -  Cl  (2k£)  +  2£-U 


-Hie  directive  gain  is  given  by  : 


r2  (  I EQ  | 2  +  1  | 2 ) 

30  I  2  sin2  [k  (H+i)]  R. 
m  in 


The  radiation  efficiency  n  is  given  for 


0.2 


(Reference  14)  by 


n(dB)  =  20  log1Q  [6.335  (|)  +  67.95  <|)2  -  693  (|)3  +  1600  (|)4 


The  power  gain  is  given  in  dB  by : 


G  (dB)  =  ri(dB)  +  10  log  g 
p  lu  a 


The  radiation  vectors  are : 


0  ,v 


sin6  I1 


I  sin  k  (H 
m 


+ 1  - 


.  jkz  cos0  . 
z)  e  dz 


Bn  A4  *  jg4 
k  sin  0 


N 


<|»,v 
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(A-118) 


( A- 1 19) 


(A-120  ) 


( A- 1 2 1 ) 

4 


(A-122) 


(A- 123) 
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.  =  -cos0  sin$  /  I  sin  (k(£-y))  e 

0  ,h  o  m 

I  A_  +  jB_ 

*  .  ,  m  5  5 

=  -coso  sm<j»  —  - ^-7 — 

sin  i|>, 
n 


jky  cos')/ 


(A-124) 


,  jt 

.  =  cos 6  J  I  sin  (k(4-y))  e 

<p  ,h  o  m 


jky  cos'll 


=  cos4>  — 


I  A_  +  j  B. 

m  5  5 


•  2 
sin 


(A-125) 


where  A4,  b4,  Ag,  and  Bg  are  given  by  Equation  A-111,  A-112,  A-113,  and  A-114, 
respectively. 


nq  +  nq  . 

0  ,v  0  ,h 


1  U,  x  j®, 
m  4+4 


A5  +  jB5 

■  2  S 

sin  s|>h 


cosG  sinij) 


( A-126) 


*■  A5  +  ^B5 


■  2  i 

sm  ip 


—  cos 


( A-127) 


The  free-space  field  intensities  are : 


I  e-jkr  [aa  + 
m _  4 

r  sir 


A5  *  jB5 
sin  'll. 


cos6  sin<|> 


(A-128) 
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j30  I  e“:’kr  A  +  j  B 
m  5  J  5 

V  =  - ; - T7C^ 

sxn  lk 


(A- 129 ) 


Hie  attenuation  relative  to  free  space  is  given  by : 


A  = 


'i^9i2  - 
'V^  'V'2 


(A-130) 


where  Eg  and  E^  are  given  by  Equations  A-116  and  A-117,  respectively. 


ADB  =  20  log 


|A| 


(A-131) 


The  electric  field  intensities  in  the  diffraction  region  are  given  by: 


j30I  e 
m 


-  jkr 


»4  *  ^4 
sin6 


*5  »  jB5 

•  2  / 
sin  ♦. 


cos0  sin<)> 


2  F 


(A-132) 


j30  I  e 
m 


■jkr 


A5  +  j  Bg 


— - - -  cosiji  2  F 

.  2  '  r 


(A-133) 


sin 


1 
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ARBITRARILY  TILTED  DIPOLE 


R>r  this  antenna: 


0  =  ir/2  -  a  ,  =  tt/2 ,  H  =  2  =  H  +  s  sin  a' 

s 


cosip  =  cos0  sin  a'  +  sin0  cos  a'  simp 


I 

m 

sin 

k 

a  + 

s) , 

-A  <  s  <  0 

I 

m 

sin 

k 

(A  - 

s)  , 

0  <  s  <  A 

Using  the  same  approach  as  in  the  previous  subsections, 
are  given  by  : 


( A-134) 


( A-135) 


(A-136) 


the  field  intensities 


e  3  ]kHoos0 

E.  =  -]60I  -  e 

0  m  r 


A, 


•  2, 

sin  ip 


(cos  a'  simp  cos6  -  sina'  sin0 ) 


(A-137) 


-R  - — — —  (cos  a'  simp  cos0  +  sin  a'  sin0)  e  32kHcos0 

V  sin2  ip 


7  -j2kHcos0  ,  _ 

+  - - — 7  e  J  (1-R  )  F 

2  v  e 


sin  ip 


/  2 
V  n  -s 


2  /2  2~ 

.  ,  2„  V  n"-sin  0  „ .  ,  .  ,  ,  vn  -  sin  0  ,  ..  . 

x  (sin  0  -  - - -  cos0 )  (cos  a  simp  v - - - sin  a  sin0) 


E .  =* 


_  e  jkHcos0 

j60I  -  e  cos  a  cos 0 

m  r 


sin2ip 


+  *h 


-jZkHcosQ  _ 7_ 


sin2  ip 


+  <1_RJ  F 

n  m 


sin2  ip 


7 _  j2kHcos0 


( A- 1 38 ) 
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where 


Ag  =  cos  (k£  cos<|<)  -  cos  k £  (A-139) 

A7  =  cos  (k£  cos  i|>' )  -  cos  k£  (A- 140) 

costy'  =  cosi|>  -  2  sin  a'  cos0  =  -cos0  sin  o"  +  sin0  cos  a'  sin<|>  (A-141) 

The  field  equations  for  a  horizontal  dipole  are  obtained  if  a'  =  0  is 
substituted  in  Equations  A- 137  and  A-138.  The  substitution  a'  =  ir/2  gives 


the  field  equations  for  the  vertical  dipole,  except  for  a  slightly  different 
surface-wave  term  (because  Fg  was  assumed  constant  in  evaluating  the 
integrals) . 


The  directive  gain  is  obtained  from: 


120  (  |  F0  |  2  +  |  F  |2) 

gd  =  - R - 2 -  (A'142) 

in 


where 


F0  = 


L  sin  ip 


(cos  a'  sin$  cos0  -  sin  a'  sin0 ) 


( A-143) 


n_ 

-  R  - - — 7  (cos  a"  sin<£  cos0  +  sin  a"  sin0)  e  cos 

V  sin  <|> 


sin 


7  - j2kHcos0  , ,  _  .  _  ,  .  2a 

7  ®  ( 1  -R  )  F  ( sin  0  - 

a.,  V  0 


\/n2-sin20 


COS0  ) 


x  (cosa'  sin$ 


JT- 


sin  9  -  sina*  sin0 ) 
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A  7 

6  -  -j2kHcos0 

.  2.  +  .  2  '  ®  _ 
sm  t|>  sin  i p 

.  .  _  .  A7  -  j2kHcos0 

+  (1-R,  )  F  - - — v  e  cos  a  cosiji 

n  m  2  . 

sm  ij> 


(A-144) 


R^^  =  R1 1  +  Re  [Z21  (cos  a'  -  j  sin  a')  (R^  cos  m'  +  j  sina')]  (A-145) 


where  R^'  and  Ry'  are  given  by  Equations  2-19  and  2-20.  Hie  self  impedance 
R1 1  can  be  calculated  from  R21  by  replacing  the  element  separation  by  2a, 
where  a  is  the  radius  of  the  element,  to  give: 


=  -30  /  1 


[—  (sin  2trr1 


S  +  Z  +  — 
z  o  2 


( A-146) 


Sz  -  Zq - 2  Sz  +  Zo 

+  sin  2it  - - 2  cos  tt£.  sin  2tr  - 

2  r2  1  r 

x  (S  2  +  Y  S  +  S  2)] 
x  o  y  y 


sin  2wr  cos  ir£  sin  2irr  sin  2irr  |sin  [2v  (— —  -  |s|)] 

+  [2  - - - - - -  ]  S  - - -  ds 

r  z  s 


-30  /  1 


—  (cos  2rrr1 
P 


s  +  z  +  — 

z  o  2 


(A-147) 
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the  free-space  field  intensities  are  given  by  : 


'0f 


60  I  e 
m 


-jkr 


cos  (k&  cosip)  -  cos  k& 

.  2, 

sin  ip 


(-cosa'  simp  cos6  h 


V 


60  1  e 
m 


-jkr 


cos  (k£  cosip)  -  cos  k£ 

■  2. 

sm  ip 


cosa'  cosip 


The  attenuation  relative  to  free  space  is : 


A  = 


sin2ip  \  / 1  Fq  I  2  +  IF^I2 


[cos  (k£  cosip)-coskil]  /(-cosa'sinipcosQ  +  sina'sin0)  n 


ADB 


20  log 


The  electric  field  intensities  in  the  diffraction  region  are: 


60  I  e 
m 


-jkd 


cos  (kicosip)  -  cos  k& 

d  4  2, 

sin  ip 

x  (-cosa'  simp  cos0  +  sina'  sin0 )  (2  F^ 


60  1  e 
m 


■jkd 


cos  (k&cosip)  -  cos  k& 


sin2  ip 


Appendix  A 


sina'  sin0) 
(A-151 ) 


(A- 152  ) 


(A- 153) 


2 

(cosa'cosip) 


( A- 1 54  ) 


( A-155  ) 


(A- 156) 


x  (cosa'  cos< b)( 2  F  ^ 
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SLOPING  LONG -WIRE 

R>r  this  antenna : 

0'  =  it/2  -  a",  4>  =  it/2,  H  =  z  =  s  sin  a'  (A-157) 

s 

cosij;  =  cos0  sin  a'  +  sin0  cos  o'  sin(j)  (A-158) 

Ihe  current  distribution  is  given  by : 

I  =  I  sin  k(l-s) ,  s  >  0  (A-159) 

s  m 

Ihe  final  expressions  for  the  field  strength  are: 


e-j*r  Ag  +  jBg 

-j3  0  — - —  - - -  (cosa'  siniji  cos0  -  sina'  sin0 ) 

sin  t 

m 

^7  +  ^ 

-  R  - - — —  (cosa'  sin<j>  cos0  +  sina'  sin0 ) 

V  sin  iji 


(A-160 ) 


A?  +  jB7  „  n  ,  r  ,  .  2  7n2  -  sin20 

+  - - — —  (1-R  )  F  (sin  0  -  - - -  cos0 ) 

.  2  .  v  e  2 

sin  n 


/  2_  .2 

x  (cos  a'  sin<}>  — — —  -  sina'  sin©) 

n 


The  first  two  terms  of  Equation  A-160  are  identical  to  Equation  54  of 
Reference  8,  and  the  third  term  represents  the  surface  wave. 


180 


f 


AD-A102  622 


I IT  RESEARCH  INST  ANNAPOLIS  MD 

APACK#  A  COMBINED  ANTENNA  AND  PROPAGATION  MODEL. (U) 

jul  ai  s  Chang »  h  c  haddocks  Fi962B-eo-< 


F/G  20/14 

:-0042 

NL 


UNCLASSIFIED 


ESD-TR-80-102 


ESD-TR-80-102 


Appendix  A 


-3  tar 


=  j30I 


cosa'  cos$ 


+  (1-R  ) 
v 


m 


A7  *  3B7 
sin2  i|/ 


•  2, 

sxn  ty 


+  R, 


A7  *  ijl 
h  sin2  <|>' 


The  first  two  terms  of  Equation  A-161  are  identical  to  aquation 
Reference  8,  and  the  third  term  represents  the  surface  wave.  In 
160  and  A-161 : 


Ag  =  cos  (ki  cosip)  -  cos  ki 


B.  =  sin  (ki  cosi|>)  -  sin  ki  cos  ip 
o 


A7  =  cos  (ki  cosij/  )  -  cos  ki 


•  I 

B?  =  sin  (ki  cost|»  )  -  sin  ki  costy 


cosij/  =*  cosip  -  2  sina'  cosS  =  -cos0  sino'  +  sin0  cosa'  sirup 


The  directive  gain  is  : 

r2  (  |  EQ  | 2  +  E  |2) 

g  *  ■  ■  -  ■  ■ 

3  30  I  2  sin2  (ki)  R 

m  in 


m 


(A-161  ) 


55  of 

Equations  A- 


( A-162) 

( A- 163) 

( A- 164) 

(A-165) 

(A-166 ) 


(A- 167) 
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where 


R. 


in 


R,  +  R  ( CZ  ) 
11  e  m 


*11 


60  [in  ki 


Ci  (2  ki)  + 


sin  2ki 
2  ki 


0.423] 


Z  =  120  (in  - -  0.60  +  in  sina')  -  jl 70 

m  27ra 


with  a  =  wire  radius • 


The  radiation  vectors  are  given  by  : 


(  A-168) 


( A-169) 


(A-170  ) 


Xm  A6  +  jB6 


Nq  *  -(cosa'  sinip  cosQ  -  sina'  sin8)  —  2 


sin  ip 


I  A,  +  jB^ 

m  6  J  6 

N.  =  cosa  cos<j>  —  - - - 

<P  Y  k  .  2 


sin  ip 

The  free-space  field  intensities  are  given  by  : 


30  Im  e’jkr  A  +  jB 
m  6  6 


'0f 


■  2  , 
sm  ip 


30  1™e_jkr  ft,  +  jB, 
tn  6  6 


v 


sin^  ip 


cosa'  cosip 


The  attenuation  relative  to  free  space  is 


iv2  * 


.  2  ,  ,2 
SaJ  +  I  1,1 


( A- 1 7 1  ) 


( A- 1 72 ) 


(cosa'  sin<p  cos8  -  sina'  sin6 )  (A- 173) 


(A- 174) 


(A-175) 
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where  Eg  and  E±  are  given  by  Equations  A-160  and  A-161,  respectively. 


AOB 


20  log 


(  A-176) 


The  electric  field  intensities  in  the  diffraction  region  are: 


30  I  e“jkd  A  +  jB, 
m  6  J  6 


— - —  (cosa'  sin$  cos9  -  sina'  sinQ)(2  f)(A-177) 
sin  t|»  V  r/ 


30  1  e-jkd  A  +  j  B 
_ m _  6  J  6 

d  .  2  . 

sxn  if; 


(cosa'  cos$)^2  F^) 


(A-178) 


TERMINATED  SLOPING-V 


Let  the  angle  between  the  antenna  wires  be  2y ,  and  the  angle  between  the 
antenna  wires  and  the  ground  plane  be  a’ .  The  projection  of  the  two  wires  on 
the  ground  surface  makes  an  angle  0’  with  the  X-axis.  The  length  of  each  wire 
is  l  meters. 

Ibr  wire  #  1 : 

9'  -  tt/2  -<*,$•-  8' 


cos  i|>  -  cos9  sina'  +  sin9  cosa'  cos  ($  -S') 


(A- 179) 


Ibr  wire  #2: 
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9'  =  ir/2  -  a'  ,  <J> •  =  -0 ' 


cos  ^  =  cos0  sina"  +  sinQ  cosa'  cos(<{)  +  6') 


(A- 180 ) 


where 


(ir/2  -  y)  <  ot'  <  ir/2  -  Y 


( A-181 ) 


.  -1 
sin 


(sin  Y 
cos  a 


(A- 182  ) 


The  height  of  the  current  element  ds  above  ground  is : 


H  =  z  =  H  +  s  sina' 
s 


( A-183 ) 


The  current  distribution  is : 


I^s) 


I2(s) 


I 


m 


e 


-jks 


( A- 184) 


WLth  the  above  current  distribution,  the  components  of  the  electric  field 
become : 
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Eq .  -  30  I  - - 

01  m  r 


cosa^  cosQ  cos  ($-0^ ) 


■jWu, 


1-e 


v  U  '  v 


(A- 185 ) 
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+ 


E 


4>2 


The  directive  gain  is  given  by: 


9d  = 


r2  (  1  Eg  | 2  +  I  | 2 ) 

30  I  2  R. 
m  in 


R.  =  120  (An  — —  +  An  sin  Y  -  0.60) 

in  2ira 


(This  is  a  free-space  approximation  where  a  is  the  wire  radius.) 
In  practice,  Rin  =600  0  is  used. 

The  radiation  efficiency  of  the  terminated  sloping-V  is : 


n  =  -2.7  dB  or  (54%) 


The  power  gain  in  decibels  is  given  by: 


G  (dB)  =  -2.7  +  10  log.,  g. 

p  lu  a 


{cosa'  cos0  cos  ($  -  6')  -  sina'  sin0)} 


{-cosa'  cos0  cos($  +  0')  +  sina'  sin©} 


N  s  — 


m 

j* 


-  jkAU. 


1-e 


■jkAh 


1-e 


(A-193 ) 


( A- 1 94 ) 


( A-195 ) 


(A-196  ) 


( A- 197  ) 


(  A-198) 
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N,  =*  — 


m 

jk 


-  jkiU 


1-e 


sin  (<|>  -  S' ) 


“1 

-jWO, 


1-e 


sin  (<j>  +  6' ) 


cosa 


(A-199) 


Ttie  free-space  field  intensities  are  given  by: 


-jkr 


"9f 


jk  30 


(A-200  ) 


-jkr 

jk  30  - -  N 

r  4> 


(A-201 ) 


The  attenuation  relative  to  free  space  is  : 


1  -V 


ADB 


20  log1Q 


( A-203 ) 


The  electric  field  intensities  in  the  diffraction  region  are: 


-  jkd 

Eq  „  jk  30  (  Nq)  {2  p_) 


( A-204  ) 


jk  30 


-jkd 

~ —  (N,  )  (2  F  ) 


(A-205) 
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TERMINATED  SLOPING  RHOMBIC37'  38 

Equations  A-185,  A-186,  A-187,  and  A-t88  still  apply  for  wires  #1  and  #2,  with 
the  understanding  that: 


-1 

sin 


.  -1 
sin 


H'  -  H 
21 


(  A-206 ) 


where 

H'  =  height  of  the  termination 

H  =  height  of  the  feed 

H"  =  height  of  the  other  two  vertexes. 

For  wires  #3  and  #4: 


H 

s 


I 


I 


H"  +  s  sina 


-jc  i 

I  e  e 

m 


_  -jki  ■ 
-I  e  e 

m 


jks 

-jks 


(A-207) 

(A-208) 

(A-209) 


Since  one  end  of  wires  #3  and  #4  are  not  located  at  the  origin,  the  radial 
distance  for  these  wires  ( r'  =  r  -  l  cos  must  be  used  in  the  exponent. 
Assume  r'  3  r  in  the  denominator.  These  approximations  yield  the  following 
expressions  for  the  fields  from  wires  #3  and  #4; 


*5  *7 

Bruce,  E.,  "Developments  in  Siort-Wave  Directional  Antennas,"  Proc. 
IRE,  Vol.  19,  No.  8,  pp.  1406-1433,  August  1931. 

38Bruce,  E. ,  Beck,  A.C.  and  Dowry,  L. R. ,  "Horizontal  Riombic  Antennas," 
Proc.  IRE,  Vol.  23,  No.  1,  pp.  24-46,  January  1935. 
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e-jkr  -jk£U2  -,_e 

E._  =  30  I  -  e  {-cosa'  cos0  cos  ($  +  0')  - ~ 

03  m  r  u 

-jk£U 

1  -e  -j2kH"  cosQ 

-  - - -  R  e 

u3 


2  .  2 
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-jk£U1 


( A-210 ) 


■T2kH*  cos6  . .  _  .  _  vn  -  sm  0  ,  .  2  A 

J  (1-R  )  F  - - -  (sm  0 

v  e  2  a 

n  cosQ 


/  2  .  2  A 
vn  -  sin  0 


cos0 )  +  sina^  sinQ 


-jkJfcU  -  jk 

1-e _  1-e 

.  °1  U3 


R  e 
v 


-j2kH"  cos8 


J  1  2  2 
1-e  . .  .  _  -i2kH"  cos0  .  .  2„  Vn  -  sin  0 

+  — 5 -  (1-R  )  F  e  J  (sm  0  -  - - -  cos0 ) 

U3  v  e  n2 


g-jkr  -jk£U1  !_e  ' 

30  I  -  e  cosa'  cos©  cos  ( 4>  +  0')  - — 

m  r  U. 


1-e  -  j2kH*'cos0 

- - -  R  e 

°4 


-jk£U2 


(A-211) 


-jk£n 


12  2  /~~2  2 
■i2kH"cos0  , .  _  .  _  vn  -sin  0  ,  .  2„  Vn  -sin  0 
( 1  -R  )  F  — -  ( sm  0  -  - - - 


v  e  2  Q 
n  cos0 


COS0  ) 


-jk£U  -jk£U 

„  .  „  1-e  ,  1-e  „  - j2kH"cos0 

sma  sin0  - - -  +  - - -  R  e 

2  4  I 


j  /~2 - 2 — 

1-e  _  - j2k H"cos©  .  .  2.  Vn  -  sin  0 

+  -  ( 1  -R  )  F  e  ( sm  0  -  - - -  cos0  ) 

U^,  v  e  2 

4  n  J 
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'♦3 


'<t>4 


e-jkr  -jk£U2 

-30  I  -  e  cosa'  sin  (d>-6'  ) 

m  r 


1  -e 


•jkiu1 
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1 


(A-212 ) 


-jkiU  -jk£U 

1-e  --j2kH"cos0  1-e  ,  „  _  ,  _  -j2kH"cos9 

-  5^  e  J  +  - - -  (1-RU)  F_  e 


U. 


U, 


e-jkr  -jMO 

30  I  -  e  cosa'  sin  (<j>  +  8'  ) 

m  r 


h  m 
-jk£n 


1-e 


u_ 


-jk£U, 


(A-213 ) 


1-e 


U. 


Rh  6 


•  j2kH"cos0 


-jkAU 


1-e 


u. 


(1-R  )  F  e 
n  m 


-  j2kH"  cosQ 


The  total  field  components  are: 


'01 


+  E92  +  E03  +  E04 


=  E, 


+  E 


4,2 


+  E  +  E 
<t>  3  <p4 


( A-214) 


(A-215) 


The  expression  for  directive  gain  is 


r2  (1E0|  ±  |E»|2) 

30  1  2  R. 

m  in 


(A-216) 


If  the  rhombic  antenna  is  properly  terminated,  may  be  taken  as  the 

termination  resistance.  In  practice,  R.  =  600  ohms  is  used.  The  radiation 

in 

efficiency  is  given  by: 
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n  =  -2.3  dB  (59%) 


and  the  power  gain  is  given  by : 


G  (dB)  =  -2.3  +  10  log 

p  *10  *d 


The  total  radiation  vectors  are : 


I  -jk*U,  -jkiU. 

M  I®  /  j  I  .  ,  .  £ 

Nq  =  —  1-e  )  (1-e 


cosa '  cos6  cos  (4»+B' )-sina' 

U„ 


cosa^  cosQ  cos  -  sing'  sin8 

U. 


I  -jk£U  -jk£.U 

N<(>  =  3^  (1-e  )  (1-e  )  cosa' 


sin  (<ji-g')  _  sin  (i fi+B' 


The  free-space  field  intensities  are  given  by: 


E 


9f 


-jkr 

j3  0k  - -  N. 

r  9 


•jkr 


'♦  f 


j30  k 
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(A-217  ) 

(  A-218) 

sin6 

(A-219) 

]  (A-220 ) 


(A-221 ) 

(A-222) 


The  attenuation  relative  to  free-space  is : 
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'V2  *  'V2 
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( A-223 ) 


( A-224) 


ADB  =  20  log 


The  electric  field  intensities  in  the  diffraction  region  are: 


Eq  =  j30k  — (N0)  (2  Fr) 


E  =  j3 Ok  — - —  (N  )  (2  F  ) 
9  d  r 


TERMINATED  HORIZONTAL  RHOMBIC 


The  following  equations  can  be  obtained  from  those  for  the  terminated 
sloping  rhombic  by  letting  a'  =0  and  6'  =  Y  • 


(A-225 ) 


(A-226) 


The  6  -  components  are  given  by : 


-jk£u, 

-]kr  ,  J  1 

0  J  —  0 

30  j  -  - - -  cos6  cos  ($-Y) 

m  r  U1 


-  j2kHcos0 


( A-227) 


n - 2  /  2  2 

_  - j2kHcos0  Vn  -  sin  0  ,  ,  2„  <Jn  -sin  0 

+  ( 1  -R  )  F  e  - - -  (sin  0  -  - - -  cos0  ) 

v  e  2  n  2 

n  cos0  n 

_  -jk£U 

_  e  -1  r  1-e  „  _  _  -j2kHcos0 

-30  I  -  - — -  cos0  cos  (A+y)  1-R  e 

m  r  U  v 


+  (1-R  )  F  e 


r 2  2  - 7“ 

-  jZkHcos©  Vn  -  sin  0  2  /n  -  sin  0 

(sin  0  -  '* _ 


(A-228) 


2 

n  cos0 


cos  0 ) 
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E_  „  =  -30  I  -  e 

93  m  r 


rjkr  -jk£U2 


-jk£U1 


cos9  cos(^-Y) 


1  -R  € 
V 

p 


■  j2kHcos0 


(A-229 ) 


+  (1-R  )  F  e 
v  e 


- j2kHcos0  Vn2  -  sin29  (sin2@  _  7n  -  sin  9 


n2cos0 


COS0  ) 


.  .  -  j 

-jkr  -3kiU  1-e  2  - j2kHcos9 

E  =  30  I  - -  e  - - -  cos6  cos ( <f>+Y )  1  -R  e 

04  m  r  u  v 

-i2kHcos0  \/n2  -  sin20  /_,_2a  'In’  -  sin^O 
+  (1-R  )  F  e  J  - : -  (sin  0  -  r  coso 

v  e  n2  cos6  nZ 


( A-230 ) 


The  ^-components  are  given  by : 


- jkiU, 
-T kr  ,  1 
a  1-e 


m  r  U, 


sin  (<}>-y)  1  +  R^e 


•  j2kHcos0 


+  (1-R.)  F  e 
n  n\ 


-  j2kHcos0 


(A-231 ) 


-jk*U_ 

«  *  -j2kHcos0 

-30  I  — -  -  sin  (<t>+Y)  1  +  R.® 

m  r  U_  n 


+  (1-R.)  F  e 

n  ttl 


-  j2kHcos0 


( A-232) 


=  -30  I 


— jk&U.  j Tt  r 

1-e  1  r  -j2kHcos0 

5 -  — -  e  sin  ($-y)  1  +  J 


+  n-RJ  p  e 

n  m 


-  j2kHcos© 


( A-233) 


-jkr  1_e’jkAa2  -3**0, 

=  30  1 - - -  e  sin  (<j>+Y) 

m  r  U 


+  (1-RJ  F  e 
n  in 


-  j2kHcos0 


1  +  Rh 


-j2kHcos0 


( A-234) 
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where 
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directive  gain  is  given  by : 

r2  (lE0|2  +  ( | 2 ) 

30  I  2  R. 
m  xn 


where  Rin  =  600  in  practice. 
The  radiation  efficiency  is: 


n  =  -2.3  dB  (59%) 


The  power  gain  in  decibels  is  given  by : 


Gp  =  ~2*3  +  10  lo*i0  gd 


The  radiation  vectors  are : 


The 


(A-239) 


( A-240 ) 


( A-241  ) 


-iktu.  ktu  Mu, 

8  1  2  2  sin  — —  sin  — — 

me  e  2  2  „  .  , 

N0  - - 7, — 7 -  cos6  sin<(>  smy 


N  . 


jk 

U1  °2 

8  X 

2 

"jkZU2  k*U1  kJtU„ 

2  sin  — —  sin  — — 

m 

e 

e  2  2 

jk 

U1  °2 

(A-242) 


siny  (sin0  cosy  -  cos«j>) 


k 
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The  free-space  field  intensities  are 


£  ,  i30*  .-Jte 

E0f  r  N0 


“♦f 


j  30  k  e 


-jkr 


The  attenuation  relative  to  free  space  is 


{ A-244 ) 


( A-245 ) 


|E/ 


'V2*  'V2 


ADB  =*  20  log 


Hie  electric  field  intensities  in  the  diffraction  region  are: 


E„  *  319*.  —  (N  )  (2  F  ) 


0 


’0' 


*♦  »  <V  ^  V 


SIDE-LOADED  VERTICAL  HALF-RHOMBIC 


Etor  wire  #1 


0  *  ir/2  -  a',  4)'  ■  0,  H  ,  «*  z  *  s  sina' 

s  l 
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( A-246 ) 


( A-247) 


(A-248 ) 


( A-249 ) 


(A-250) 


r 

1 

i 

f 

l 

I  ESD-TR-8  0-102 


cos4<1  =  cosQ  sina'  +  sin6  cosa'  cos$ 

Itor  wire  #2 : 
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(A-251  ) 


0'  =  rr/2  +  a,  <j>'  =  0,  H  _  =  (i  -  s)  sina' 

s2 


(A-252) 


cosi|<2  =  -cos0  sina'  +  sin0  cosa'  cos$ 


(A-253) 


The  current  distribution  on  the  two  wires  is  given  by : 


I1(S)  =  Xme 


■jks 


( A-254 ) 


I_(s)  =  I  e 

z  m 


-jk l  e 


-jks 


( A-255 ) 


Following  the  same  approach  used  in  the  previous  subsections,  the  final 
expressions  for  the  components  of  electric  field  strength  are: 


-jkr 


'01 


30  I 


m  r 


+  (1-R  )  F  F 
v  e  2 


-cosot'  COS0  COSp 


F,  -  R  F„ 
1  v  2 


( A-256 ) 


y/n2  -  sin^0  ,_J_2„  v4^~sin^0 


2  e 

n  cos 


(sin  0  - 


COS0  ) 


+  sina'  sin0 


n. - t~ 

F1  +  Rv^2  +  FeP2  ^sin2®  “  “ - 23in  -  cos8) 

n 
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-02 


-jkr  -jk IV 

-30  I  5 -  e 

m  r 


cosa'  cosQ  cosf 


F  -  R  F  e 
2  v  1 


-jk2k£  sina' 


where 


+  (1-R  )  F  F 
v  el 


l 


2  .  2a 
sm  9 


n2cos0 


/  2  2 

.  .  2a  Jn  sin  9  -j2ki  sina'  cos0 

x  (sm  9  -  ’ - - -  cos9 )  e 


+  sina'  sin9 


F2  +  Rv  P,  e 


■j2ki  sina'  cos9 


(A 


/2  2~ 

.  2n  vn  -  sin  0  -j2k£  sina'  cos0 

+  (1-R  )  F  F„  (sm  0  -  v - - -  cos0)  e 

v  e  1  2 


>1 


-jkr 


-30  I 


m  r 


cosa'  sin$ 


"-V  F.  P2 


.-3kr 


%2  =  -3°  1 


m  r 

-jkiU 


cosa'  sin<t> 


-jk  *U 


F2  6 


+  RhF'e 


1 


.  -jWO, 

+  "“V  Fm  F1  e 


(A 


(A 


U1  *  l-cos^  =*  1-cos0  sina'  -  sin9  cosa'  cos<j> 


(A 


U2  =  1-cos4'1  +  2  cos9  sina'  =  1-cost2 


(A 


=  1  +  cos9  sina'  -  sin9  cosa'  cos<j> 


cos9 


-257) 


-258) 


-259) 


-260) 


261 ) 
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-jkiU 


1  -  e 


U. 

1 


,  i  =  1,  2 


(A-262 ) 


The  total  field  components  become: 


I  -  j  kr 

E  =  E  +  E  =  _  p  ' 

E0  E01  E02  r  F0 


(  A- 2 63  ) 


E  =  E  +  E 
<J>  $1  <p2 


1  -jkr 
30  m  e 


(A-264 ) 


where 


Fq'  =  -cosa'  cos0  cos<j> 


-jkAU, 


-jkZUn 


F1  +  F2  @ 


"  Rv  (F2  +  F!  6 


(A-265) 


+  (1-R  )  F 
v  e 


y/n2  -  sin29 


+  sina"  sin0 


n  cos0 


„  2a  -jkin 

.  .  2„  vn  -  sxn  0  2, 

x  (sin  0  -  - - -  cos0)  (  F2+Fie  ) 


■jkAU, 


-jk«, 


F  -  F  e 
1  2 


+  Rv  (F2  "  F1  6 


+  (1-R  )  F  (sin20  - 
v  e 


y 2 

Vn  -  si 


2 

sin  0 


-j kiU, 


cos0)  ( F2  -  F1  e 


F,"  »  cosa'  sin$ 

<P 


-jWO  -jkAO 

F,  +  F2  e  ^  \  (F2  +  F1  6  > 


(A-266) 


-jwo. 


+  (1"V  Fm  (F2  +  F1  « 


> 


The  directive  gain  of  the  sid&- terminated  vertical  half-rhombic  antenna  is 
then : 


30  (  |  FQ'|2  +  I  F  '|2) 


in 


( A-267) 
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where 


60  (In  — —  +  In  sin  a' 
2na 


=  wire  radius 


-  0.6) 


Equation  A-268  is  exact  when  the  ground  is  perfectly  conducting. 


The  radiation  efficiency  is : 


(A-268) 


n  =  -2.3  dB  (59%) 


(A-269 ) 


The  total  radiation  vectors  become: 


N0  “ 


I 

—  -F.  ( cosa' 

3k  1 


cosQ  cos4>  -  sina'  sin0 ) 


-  F2  e 


-jk£U1 


(cosa'  cosQ  cos$  +  sina'  sin0) 


( A-270) 


-3WU, 


•jWD, 


=  —  -cosa"  cos9  *cos$  ( F1  +  F2  e  )  +  sina'  sin0  (F^  -  e  ) 


Im  -3k£U1 

N  «  -  —  [F,  +  F_  e  ]  sin$  cosa" 

4>  jk  1  2 


( A-271 ) 


The  free-space  electric  field  intensities  are  given  by: 
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!0f  =  j3°  k  — T  Ne 


( A-272 ) 


V  "  330  k  — 


(A-273) 


The  attenuation  relative  to  free  space  is  given  by: 


<v2  +  'y2 

lE9f|2  +  'V2 


( A-274 ) 


»  20  log 


( A-275 ) 


The  electric  field  intensities  in  the  diffraction  region  are: 


E*  =  330k  — —  (H0)  (2  Fr) 


(A-276) 


E,  -  j  30  k  — - —  (N  )  (2  F  ) 
9  d  9  r 


(A-277) 


HORIZONTAL  YAGI-UDA  ARRAY 


Ebr  this  antenna : 
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Q'  =  TT/2,  4>  =  0,  a'  =  0,  H  =  H,  s  =  X ,  cos<|>  =  sin©  cosifi 

s 


( A-278  ) 


I  (s)  =  I  .  sin  k  (l.  -  |x|),  i  =  1,2,...,  N 

l  mi  l 


^'<ri  cos  (k i  cosip )  -  cos  k£  . 

e  nkHcoso  i 

E- .  =  -36O  I.  -  e  - ; - 

01  l  r.  .2. 

l  sin  ip 

a  ,  .  _  -  j2kH  COS0  , 

x  cos0  ( 1  — R  e  ) 

v 


( A-279 ) 


,  Jn  -  sin  0  --)2kHcos0  2  .  \Jn  -sm  8  „ . 

+  (1-r  )  * - - -  e  J  F  (sm  8  -  - - -  cos0  ) 

v  2  e  2 

n  n 


^ri  .  cos  (kJt  cosij;)  -  cos  k£ . 

e  ikH  cos©  l 

j60  j  __  e  - — - 

l  sin  p 


-j2kHcos0 

X  1  +  \  6  J  + 


•  2, 

sin  p 

-■iTVBfncfl 


( A-280 ) 


f*  Vi 

is  the  distance  between  the  base  of  the  i  element  and  the  far-field 
point,  which  is  related  to  r1  by: 

=  r1  -  y^  sin0  sin<}>  (A-281) 
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where 


y  =  0,  y  =  d  ,  y  =  d  +  d  ,  .  .  .  y  =  £  d  (A-282) 

1  2  13  12  n  i=1  i 


The  total  field  components  can  be  obtained  by  summing  the  contributions  from 
all  the  elements  to  give: 


-jkr 


1 


e  jkHcos0  cosij) 

E  «  - j60  -  e  - 

i=1  0i  r  2 


1 


sin  ^  L 


- j2kHcos0 

cos0  ( 1-R  e  )  +  ( 1-R  ) 

v  v 


JT. 


sin  0  -j2kHcos0  2 

-  e  F  (sin  0 


rr 

V  n  - 


2 

sin  0 


COS0) 


(A  +  jB  ) 
11  11 
( A-283) 


■jkr 


1 


e  jkHcos©  sin<}> 

E  =*  j60  -  e  - 

i=1  <J>i  r  2 


1 


sin  \|; 


1  +  R  e 
h 


-  j2kHcos0 


( A-284) 


+  (1-R  )  F  e 
h  m 


-j2kHcos0 


(A  +  jB  ) 
11  11 


where 


L  +  jB  =  J 
11  11  1=1 


lei 
mi 


jk  y  sin0  siniji 


[cos  (k£  -  cos^)  -  cos  k£  ]  (a-285) 

i  i 


The  input  impedence  looking  into  the  feed  point  of  the  active  element  can  be 
determined  by  considering  the  circuit  relations: 
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211 

+ 

V 

ziT 

Z  1  2  + 

V 

V  • 

.  .  .Z,  + 

IN 

V 

Z1N ' 

Tb1 

Z21 

+ 

V 

z2l' 

Z22  + 

V 

Z22^  * 

.  .  .Z„  + 

2N 

V 

Z2N 

Ib2 

ZN1 

+ 

V 

zn/ 

ZN2  + 

V 

ZN2^  * 

#  •  •  z  + 

NN 

V 

ZNN 

• 

• 

XbN 

where  R^,  the  reflection  coefficient  for  horizontal  polarization  with  9=0°, 
given  by: 


V 


and 


1  -  n 
1  +  n 


k  -  k. 

iT+T 


( A-287) 


2ii 


open-circuit  self- impedance  of  the  ith  real  element,  calculated  by 
Equation  A-289  with  l  - 


Zii 


open-circuit  mutual  impedance  between  the  ith  and  jth  real  elements, 
calculated  by  Equation  A-290  with  d^j  =  | y^  -  y^ | 


Zii' 


■  th 


open- circuit  mutual  impedance  between  the  i  real  element  and  its 
own  image,  calculated  by  Equation  A-238  with  l  =  and  d  =  2H 


V 


=  open-circuit  mutual  impedance  between  the  ith  real  element  and  the 
image  of  the  real  element,  calculated  by  Equation  A-290  with  d  ^  ^ 

|/h2  +  (y.  -  y.)2 


2  =  z - -  -■  rr  e"j2kZ  [K(U  )  -  2  K  (U  )]  +  ej2kA  [K(V  )  -  2K(V  )] 

m  1  -  cos  2  k*  I  o  1  o  1 


+  2  [K(U  ')  -  K(U, )  -  K  (V  )  ]  +  2  K(U  ')  (1  +  cos  2  kJ i) 
oil  o 


(A-298 ) 


is 
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where 


K  ( x)  =  Ci  (x)  -  j  Si(x)  =  /*  c°-— dy 


/•x  sm  y 

3  J  - -  dy 

1  o  y 


( A-288a) 


k  [  \/d2  +  4 l2  -  21} 


(A-288b) 


Vq  =  k  [  \j&2  +  4A2  +  21] 


( A-288c) 


U  '  =  kd 


( A-288d) 


k  [  \/d2  +  £2  -  i] 


( A-288e) 


V4  =  k  [ 


\/a2  +  i2  + 


4] 


(A-288f ) 


d  =  2H 


(A~288g) 


Equation  A-288  is  also  used  to  calculate  the  self  impedance  if  d  =  a  2,  where 
a  is  the  radius  of  the  dipole,  i.e.: 


Z 


1 1 


( A-289) 


Since  Equation  A-288  is  not  valid  for  k£  =  irP,  P  =  1,  2,  . the  more  exact 
formulas  in  Reference  9  should  be  used.  lb  save  computer  time.  Equation  A-288 
can  be  used  satisfactorily  in  most  cases. 
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The  mutual  impedance  between  two  parallel  dipoles  of  different  lengths, 

•5Q 

referred  to  the  antenna  base,  is  found  from: 


60 


m  cos  W_  -  cos  W 

2  1 


-jW, 


e  [K(U  )  -  K(U  )  -  K( U  )  ] 

o  1  2 


jw,  -jw 

+  e  [K(V  )  -  K  (V  )  -  K  (V  )  ]  +  e  [K(U  ')  -  K(U  )  -  K  (V 

0  12  o  1  2 

(A 


+  e  [K(Vq')  -  -  K(U2 )]  +  2  KCfcM  (cos  W1  +  cos  w2> 


where 


K( x)  =  Ci(x)  -  j  Si(x) 


(A 


U  "  k  t  \  k.  2  +  U.  +  l.  )2  -  (1.  +  l  .)] 
o  V  i]  13  13 


(A 


V  =  k 
o 


c^a..2  *  a.  *  «.,2  +  at  ♦  J . 


)] 


(A 


King,  R.W.P.  and  Mi,  T.T.,  "Currents,  Charges,  and  Near  Fields  of 
Cylindrical  Antennas",  Radio  Science,  Vol.  69D,  No.  3,  pp.  429-446, 
March  1965. 


>] 

-290  ) 


-290a) 

-290b) 

-290c) 
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k  [\fd iT^i-v5- 


v  '  -  k  (\/d..2  +  a.  -  i.)2  +  a.  -  £.n 

o  \/  13  13  13 


U1  -  k  [ 


\l  d.  ,2  +  £  . 2  -  £.] 

V  13  1  1 


V1  -  k  f\/dij2 +  £i2  +  v 


°2  *  k  [  \/di32  +  S'  -  S1 


v2  =  k  [ 


,/d..2  +  £/ 

V  13  3 


+  V 


W,  =  k  (£.+£. ) 
1  1  3 


W  =  k  ( £ .  -  £  . ) 

2  13 


W 

o 


k  d.  . 
13 


dij 


distance  between  the  two  elements  considered. 


Appendix  A 

( A-290d) 

( A-290e) 

( A-290f ) 

(A-290g) 

( A-290h) 

( A-290i ) 

( A-290  j ) 

(A-290k) 

( A-290 1) 


It  is  easy  to  show  that  Equation  A-288  can  be  obtained  from  Equation  A-290  by 
letting  £^  =  £^  =*  £. 


Solving  for  Jbi  from  Equation  A-286,  the  current  maxima  can  be  calculated 
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from : 


I  . 
mi 


bi 


sin  k£ . 


i  -  I,  2, 


N 


(A-291  ) 


The  input  resistance  is  given  by  : 


Again,  the  directive  gain  is  found  by: 


ri2  (  'V2  *  |E<j,|2) 


30  I,  _  R. 
b2  in 


{ A-293) 


Since  1^,  Eg,  and  E^  are  directly  proportional  to  V2,  V2  cancels  in  Equation 
A-293.  Eg  and  E^  are  given  by  Equations  A-283  and  A-284  and  Rin  by  Equation 
A-292.  I^2  is  found  in  Equation  A-286.  The  radiation  efficiency  of  a  Yagi- 

Uda  array  is  given  by : 


n  =  -1.8  dB  (65%) 


( A-294 ) 


The  power  gain  in  decibels  is  given  by: 


G  =  -1.8  +  10  log.  _  g .  (A-295  ) 

p  310  d 
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The  radiation  vectors  are  obtained  by  summing  the  contributions  from  all  the 
elements  to  give: 


N  =  -cos<p  cosQ 
9 


2  I  cos  (ki  cos^)  -  cos  ki  jky  sin9  sinij) 
mi  i  i  i 


i-1  k 


2 

sin 


(A-296) 


N  =  sincf) 


2  I  cos  (ki  cos4>)  -  cos  ki  jky  sin0  sin<)> 
mi  i  i  i 

-  -  e  ( A-297) 


i=1  k 


2 

sin  ip 


The  free-space  field  intensities  are  given  by : 


-  jkr 


9f 


j30k 


N 


r  9 


(A-298) 


-jkr  jky  sin0  sin<]> 

j60  e  cos<j>  cos0  N  i 

- E  i  e  [cos  (ki  cosi|>)  -  coski  ] 

r  2  mi  i  i 


-jkr 

j60e  sin$ 


4>f 


r  2 

sin  i=1 


sin  f  i=1 

I 


jky  sin0  sin<{> 
i 


1  e 
mi 


(A-299) 


x  [cos  (ki  cos  ip)  -  cos  ki  J 
i  i 


The  attenuation  relative  to  free-space  is 
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'V2  *  i*V2 


"W2  ♦  'V1' 


ADB  =  20  log 


the  electric  field  intensities  in  the  diffraction  region  are: 


-  j  kd 

-j60  e  cosg  cos0 


2 

sin  g 


N  jky  sin6  sing 

I  I  e  i 


i=1 


mi 


x  [cos  (ki  cosg)  -  cos  ki  ]  >  (2  F  ) 
i  i  \  r 


-  j  kd 

j6  0e  sing 


d  2 

sin  g  I  i=1 


jky  sin0  sing 
N  i 

E  I  e 
mi 


x  [cos  (k£  cosg)  -  cos  ki  ]  >  (2  F  ) 

i  i  r 


HORIZONTALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 


40 


For  this  antenna. 


0,  0 '  *  ir/2. 


0,  s  *  x,  H  ”  H, ,  cosg  «  sin0  cosg 
s  i 


I  sin  k  U-s)  -  I  .  sin  k  (A  -  |x|) 

(Ti  mi  i 


(A-300 ) 


( A-301) 


( A-302 ) 


( A-303) 


( A-304 ) 


g"  =  angle  between  the  far- field  point  (0,  g)  and  the  array  axis. 


40 

Ma ,  M. T.  and  Walters,  L.  C. ,  Computed  Radiation  Ibtterns  of  Log-Periodic 
Antennas  Over  Lossy  Plane  (around,  ESSA  Technical  Report  IER  54-ITSA  52, 
Boulder,  CO,  1967. 
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cost"  =  cos0  cos0”  +  sin0  sin0"  sint 


The  lengths  of  the  dipole  elements  are  related  by: 


l 

-  Ti  n  *  1,  2,  .  .  .  N-1;  0  <x  <1 

n 

The  spacing  between  dipole  elements  is  related  to  the  element 


d 

n 

l  +1 
n 


(1 -T )  cota 


i=1 

r  *  r  -  (  l  d  )  cost"  *  r  -  y  csc0"  cost" 
i  1  n=*1  n  1  i 


i  *  1 ,  2,  .  .  .  . ,  N 


y  =  0 

1 


1r1 

H  =»  H  +  (  l  d  )  cos0"  =  H  +  y  cot  9" 
i  1  n*1  n  1  i 


The  field  components  contributed  by  the 


,  th 


element  are: 


(A-305 ) 

(A-306) 

lengths  by : 

(A-307) 

( A-308) 

(A-309) 
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E9i  " 


-jkri 

j60I  .e  cos  (k i.  cosip)  -  cos  k i. 

mi  1  i 


r .  .2, 

l  sin  ip 

.  -j2kH.  cos9 

x  I  cos0  ( 1  -  R  e  1  ) 


COS!{> 


*  <1-Rv> 


rr~. 

vn  -  si 


sin20  _  .  .  2n 

2 -  Fe  (sm  9  - 


/n2  ^ 


sin  9 


cos9 ) 


-jkr. 

j60  I  .e  cos  (k£.  costp)  -coskJl . 

_ _ mi _  _ l _ l 

r .  .2. 

l  sin  ip 


sinip 


1  +  V 


+  ( 1  -R  )  F  e 
h  m 


-j2kH.  cos9 


Summing  the  contributions  from  all  the  dipole  elements  gives  : 


_  i60 


-jkr 


cosip  cos9 

.  2, 

sin  ip 


j60 


e 


iini.  s 

.  2  ,  S 
sin  ip 


(A-310 ) 

-  j2kH^cos0| 

e  J 

j2kHL  cos6 

(A-311) 


(A-312  ) 

(A-313) 


where 
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jky  csc0"  cost" 
N  i 

Die 


9  i*1  mi 


^cos  (k£ 


cost)  ~  cos  k£ 


J 


-j2kH  cos0 


i 

x  1  -  R  e 
v 


-  j2kH 


+  (1-R  )  e 
v 


COS0  f 
i  Vn 


2  2 
sin  0 


_  n.  — 2~ 

,  .  2„  vn  -  sm  0 
x  (sin  0  -  - - -  cos© 

n 

jky  csc0"cost" 

N  i 

Die 
i=1  mi 


'] 


e  2 

n  cos0 


-j2kH  cos0 
i 

x|1  +  R  e  +  (1-R  )  F  e 

h  h  m 


(■ 


£cos  (k i  cost)  “  cos  k&i^ 
COS0  I 

1  J 


■j2kH 


I  . 
mi 


Ibi 


sin  k i . 

i 


(A-314  ) 


( A-3 1 5 ) 


(A-316) 


The  directive  gain  of  the  array  is  given  by : 


120 

R. 

in 


2  2  2  2 
cos  t  cos  9  I Sfl |  +  sin  t  I s 


£ 


■  4, 

sin  t 


( A-3 17  ) 


Calculation  of  the  input  impedance  over  lossy  ground  is  extremely  difficult. 
The  main  problem  is  the  determination  of  the  dipole  feed-point  currents 
accounting  for  the  effects  of  lossy  ground.  Rjllowing  Reference  8,  two 
different  types  of  ground  will  be  considered  separately. 


a.  fbr  ground  with  relatively  low  conductivity,  such  as  sea  ice,  polar 
ice  cap,  or  poor  ground,  the  presence  of  the  ground  will  be  ionored  for  the 
purpose  of  determining  the  dipole  feed-point  currents.  The  poor  ground  is  not 
likely  to  have  a  substantial  effect  on  the  current  distribution.  Wien 
calculating  the  input  impedance,  the  imperfect  ground  effect  is  taken  into 
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account  by  adding  necessary  terms  pertaining  to  the  horizontally  polarized 
reflection  coefficient. 

b.  R>r  ground  with  relatively  high  conductivity,  such  as  sea  water  or 
fresh  water,  the  feed-point  currents  are  determined  as  if  the  antenna  were 
above  a  perfectly  conducting  ground  using  R^  to  account  for  the  imperfect 
ground . 

Case  1  (Ground  with  Relatively  low  Conductivity) 

The  matrix  of  dipole  feed-point  currents  is  given  by: 


[IJ  =  {[0]  +  [Y  J  [Z  ]}_1  [1]  (A-318) 

3  L  cl 

where 

[U]  *  N  X  N  unit  maxtrix 

[I]  =  [10  0...  0]T 

(The  current  source  at  the  feed  point  of  the  shortest  dipole  is  normalized  to 
unity.) 

[ZaJ  =  open-circuit  antenna- impedance  matrix  given  by: 


z. 

Z 

Z 

1  la 

12a 

INa 

Z21a 

• 

• 

• 

• 

2  22a  * 

.  .  .  Z 

2Na 

. 

* 

• 

% 

z  „ 

Z 

i  •  •  z 

Nla 

N2a 

NNa 

(A-319) 
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•ftie  main-diagonal  terms  represent  the  self- impedance  of  the  dipoles,  which  can 
be  calculated  by  Equation  A-289.  The  off-diagonal  elements  represent  the 
mutual  impedances  between  dipoles.  These  mutual  impedances  can  be  calculated 
by  Equation  A-290  with : 


j?1 

d  =»  l  d  . 
ij  n=i  n 


The  response  voltages  appearing  at  the  dipole  feed  points  are  found  by: 


[va]  = 


rzj  [ij 


( A-320  ) 


The  input  impedance  Zin  of  the  entire  array  is  numerically  equal  to  the 
voltage  V1  across  the  feed  point  since  the  input  current  has  been  assumed  to 
be  unity.  Including  the  effects  of  lossy  ground: 

Z'=V=Z  I  +  Z  I  +....+Z  I 
in  1  11a  la  12a  2a  iNa  Na 

( A-321 ) 

+  R  "  I  Z  I  =  R  "*  +  j  X 

h  i=*1  1  (N+i)a  ia  in  in 

where 


V  '  ,Rh  'lth  9  * 01 


( A-322) 


216 


ESD-TR-8  0 - 1 0  2 


Appendix  A 


Z1 (N  +  i)  a  =  open-circuit  mutual  impedance  between  the  first  (shortest) 
dipole  and  the  image  of  the  itk  dipole. 

Hie  mutual  impedance  can  be  calculated  from  Equation  A-290  with  : 


dij 


2  3"1 

4  H  +  (  Z 
n=1 


2 

d  )  +  4H 

n 


j-1 

(  Z  d  )  cos9 " , 
n=1  n 


d 

o 


0 


(A-323) 


[Y,]  = 


where 


short-circuit  admittance  matrix  (N  X  N)  associated  with  the  trans¬ 
mission  line  feeding  the  array 


[Y< 


Y11f  Y12f  °* 


Y  Y  Y  0 
1 2  If  22f  23f 


0  Y 


3Sf 

i 


Y11f  =  -jY0  oot  kd1 


Y22f  *  -^Yo  (GOt  kdl  +  COt  kd2) 


0 

0 


(N-1  )Nf 


Y(N-1 )Nf  YNNf 


Y(N-1  )  (H-1)f  =  -j  Yo  (COt  k  Vi  +  COt  k  dN-2) 


( A-324) 


( A-325) 


( A-326 ) 


( A-327 ) 
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NNf 


V'  -  j  V  cot  kd 
T  o  N-1 


Y  ' 
T 


=  Y 


cos  kd..  +  jY  sin  kd.. 

_ N  J  o  T _ N_ 

o  Y  Z_  cos  kd.,  +  j  sin  kd., 
o  T  N  N 


Y12f  =  Y21f  =  -j  Yo  CSC  kdt 


Y  =  Y  =  -j  Y  esc  kd. 
23f  32f  J  o  2 


Y(N-1 )Nf  YN(N-1)f  ~  “ j  Yo  CSC  kdN-1 


(A-328) 


( A-329  ) 


( A-330) 


(A-331  ) 


( A-332  ) 


Yq  =  characteristic  admittance  of  the  transmission  line 

ZT  =  termination  impedance  connected  to  the  last  (longest)  dipole  at  a 
distance  dR  =  &n/2  • 

case  2  (Ground  with  Relatively  High  Conductivity) 

The  dipole  feed-point  current  in  this  case  is  given  by: 


[Ia'J  »  {[0]+  [Yf]  (ZaJ  +  [Yf]  CZa']}_1  [II 

where 


Z1(N+1 )a 

Z1(N+2)a 

*  *  *  *  Z1(2N)a 

tZa'] 

- 

Z2(N+1 )a 

Z2(N+2)a 

*  *  *  *  Z2(2N)a 

ZN(N+1 )a 

ZN(N+2)a 

•  *  '  *  ZN(2N)a 

( A-334) 
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[Z  ']  =  open-circuit  mutual- impedance  matrix  between  the  real  dipoles 

and  their  images 

Zi(N  +  j)a  =  open-circuit  mutual- impedance  between  the  ith  real  element 
and  the  image  of  the  j  element. 


All  the  mutual  impedances  can  be  calculated  by  Equation  A-290  with 


d 

in 


With  the  dipole 
by: 


3?1 

+  ( 

I 

n=  i 

jpl 

l 

d 

n=  1 

n 

( A-335 ) 


'■] 


1/2 


feed-point  currents  determined,  the  input  impedance  is  given 


Z  " 
in 


V 

Z  I  '  +  Z 

X  +  .  •  . 

.  +  Z 

1 

11a  la  12a 

2a 

INa 

+  R  ' 

|z  1  '  + 

Z  I 

+  •  .  . 

h 

1  (N+l )  a  la 

1 (N+2 )a  2a 

N 

=  l 

(Z  +  R  "  Z 

)  I'  = 

R  "  + 

i=1  lia  h  1  (N+i)a  ia  in 


I  ' 
Na 


+  Z  I 

s  ( 2N )  a  Na 

Ta-336) 

j  X  -- 
in 


The  radiation  vectors  are  given  by : 


N  =  -COS0  COS<{>  T 
0  i=1  k 


2  1  cos  (k£  cosi())  -  cos  k£  jky  csc0"  cos^" 
mi  i  i  i 


2 

sin 


( A-337 ) 
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=  s 


in*  .? 


2  I  cos  (k l  oosi|i)  -  cos  kJt  jky  cscQ"  cos’t" 
mi  i  i  i 


i=  1  k 


2 

sin  ip 


-  jkr 

®0f  =  330k  —  N 


r  6 


e'jkr 

V  "  330k  — 


The  attenuation  relative  to  free  space  is  : 


A 


/  2  2  ,  .  2  .  .2  ,  ,2 

/cos  Q  cos  ^  |Sel  +  sin  ^  |S^| 


7 


2  2  2 
cos  „  cos  .  +  sin 

0  ffl 


where 


S 

f 


jky  cscQ"  cosip” 
i 

I  e 
mi 


[cos  (k Z  cost)  -  cos  k£ 


4 


The  electric  field  intensities  in  the  diffraction  region  axe: 


_  2 


60  e 


■  jkd 


cos9  cost 

•  2, 

sxn  * 


( Sf )  (2  Fr) 


=  m. 


-jkd 


(A-338) 


( A-339) 


(A-340) 


(A-341 ) 


( A-342) 


( A-343  ) 


sin 


(Sc)  (2  F  ) 


(A-344) 
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VERTICALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 


The  lengths  of  the  dipole  elements  are  related  by: 


l 

=  t;  H  *  1,  2 . 8-1j  0  <  t  <  1  (A-345) 

Vi 


The  spacing  between  the  dipole  elements  is  related  to  the  element  lengths  by: 


d 

n 


l 


n+1 


_ 1  -  T _ 

sin  (<*2  +  a^)  -  tan  a3  cos  (a2  +  a3  ) 


The  angles 


a2,  and  a3  are  related  by: 


sin  a1 

cos  (a1  +  a2  +  a3> 


sin  a2 
cos  a3 


(A-346) 


( A-347  ) 


The  height  of  the  center  of  the 


ith 


element  above  ground  is  determined  by: 


H. 

l 


tan  (a  +  a  ) 

l  - - - - - 

i  tan  (a  +  a  +  a  )  -  tan  (a  +  a  ) 


(A-348 ) 


Since  the  dipole  elements  are  parallel  to  the  Z-axis,  only  the  9-component  of 
the  field  exists.  The  final  expression  for  the  electric  field  strength  is: 


j60  e 


-jkr. 


*12  1  ^12 

sin0 


and 


( A-349 ) 
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E,  =0 


( A-350 ) 


where 


I  jkH  cos0 

N  bi  i 

A  +  jB  *  I  -  e 

12  12  i=1  sin  k£ 

i 


cos0 

Icos  (k£ 


cos0 )  -  cos  k£ 


( A-351 ) 


- j2kH  cos0  .  /  2  2 

I  i  2  l/n  -sin  9  -j2kH  cos0 

x  I  1+R  e  +  ( 1  -R  )  (sin  9  -  * -  cos 9)  e  i 

1_  v  v  2 

n 


—  ~P  2  _  ]  jky  sec  (a  +  a  )  cosip' 

x  [1"j  y»P  e  sin  9  erfc  ( j  )]J  e 


cos4<'  =  cos9  sin  (a2  +  a3>  +  sin6  cos  (a2  +  ct3)  sin<*>  (A<-352) 


Yi  =  y-coordinate  of  the  feed-point  of  the  ith  dipole. 


ip'  =  angle  between  the  far- field  point  (0,  <p )  and  the  array 


(9"*  "  I'  “2  -  V  3f>  * 


The  short- circuit  admittance  matrix  for  the  transmission  line  is  given  by 
Equation  A-324. 

The  elements  of  the  open-circuit  impedance  matrix  for  the  dipoles  are  given 
by  : 
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[Z  J  =  . 

cL  • 


Z 1  la  Z12a . Z1Na 

Z21a  Z22a . Z2Na 


( A-353 ) 


ZN1a  ZN2a . ZNNa 


The  mutual  impedances,  when  two  dipoles  are  either  arranged  in  echelon  or 
collinear  as  shown  in  Figure  A-2  are  as  follows.41  The  diagonal  elements  are 
found  by  Equation  A-20. 

For  echelon  arrangement : 


ija  cos  k(£.  -  Si.)  -  cos  k(Jl.  +  £,.] 
13  13 


..  +  l.) 


-jk  (£ ,-h . .) 


1  13 


[K(U  )  -  K(U  )  ] 
o  1 


(  A-354 ) 


jk(  i  ,-h  .  . ) 


[K  ( Vq  )  -K  (V  ^ )  ]  +  e 


-jkU.+h.  .) 


[K(U  ')-(K(U_ )] 

O  2 


-jkU.-2Jl  -h.  .  ) 


+  e  3>  [K  ( V  " )  -K  (V  )  ]  +  e  1  3  13  [-K  (U  )  +K  (U_ )  ] 


o  2 


V  3 


-jkU.-2£.-h.  .) 


-jk(Z.+2l  +h . . ) 


+  e  1  3  13  [-k(V,)+K(V,)]  +  e  1  3  13  [-K(uo)+K(U/l )  ] 

13  2  4 


jk  (l  +2£.+h  )  - jkh  .  . 

+  e  3  3  [-K(V2)+K(V4)]  +  2  cosk£i  e  13  [KCW^-KO^ )] 

jkh  -jk(2£  ,+h  ) 

+  2  coskA  e  3  [K(Y  )-K(Y  }]  +  2  cosk*.  e  3  13  [K(W  )-K(w  )] 

*  «  1  1  4  J 


jk(2i  +  h  ) 

+  2  cosk£i  e  3  [K(Y2)  -  K(Y3)]} 


King,  H.E.,  "Mutual  Impedance  of  Uhequal  Length  Antennas  in  Echelon," 

IRE  Trams .  Antennas  and  Propagation ,  Vol.  AP-5,  No.  3,  pp.  306-313,  1957. 
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a)  COLLINEAR 
ARRANGEMENT 


b)  ECHELON 

ARRANGEMENT 


Figure  A-2.  Two  parallel  dipoles  of  arbitrary  lengths 
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R>r  collinear  arrangement : 


Z .  . 
x:a 


where 


30 


cos  k  (i  . -i .)  -  cos  k  (i.+i.) 

3-3  il 


-jk  (i . -h .  . ) 

e  1  13  [K(U  )-K(U,)] 

O  1 


jk(i ,-h . . ) 

+  e  1  13  An 


h.  .  t  i  .  +  i  . 
_ i _ I 


h. .  -  i . 
3-3  i 


jk  (i.+h.  . ) 


-jk(i.+h.  . ) 


+  e 


[K(V  "* ) -K (V  )  ]  +  e 
o  i 


in 


h. .  +  i .  +■  i  , 
3-3  3.  3 

h.  .  +  i  . 

I  30  l  J 


-jk(i.-2i  -h. .) 

+  e  1  3  13  [-K(U1)+K(U3)]  +  e  1  3  lj  In 


jk(i  -21  -h. .)  [h. .-i.+i.  I 

3-3  i  -I  1 


jk(i.+2i.+h. .)  -jk  (i  .  +2i  ,+h  .  . ) 

+  e  1  3  3  [-k(v2)«(v4)]  +  e  1  3  13 

L  3-j  i  j. 


jkh .  . 


+  2  coski.  e  13  [-K(Y  )+K(Y  )]  +  2  coski.  e  LJ  in 


“jkh 


jk( 2i ,+h . .) 

+  2  coski t  e  3  13  [K(Y2)-K(Y3)] 


-jk(2i ,+h . . ) 

+  2  coski .  e  3  13  in 

1 


h . . +2i  . 

U _ I 


h. .  +  i . 
,  3-j  3 


( A-355 ) 


uo  -  k  [  vvT(hij  -  v2  +  (hij  -  v] 


vq  =  k  [  +  (h  -  V2  -  (h  -  i,)] 


[  ^  -  (h±j  +  i/  “  Ch^  ti.)] 
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V  ' 
o 


U1 

V1 


U2 

V2 

U3 

V3 

°4 

V4 

W1 

Y1 

W2 

Y2 

*3 

Y3 


k  [ 
k  [ 
k  [ 

k  [ 
k  [ 
k  t 
k  [ 
k  [ 
k  [ 
k  [ 
k  [ 
k  [ 
k  [ 
k  [ 
k  [ 


13 


(hij  .£.)  +  (hij  +tl)] 


v4,  2  + 


ij 


(h. .  -  ZJ  +  Z .)  +  (h. .  -  Z. 
13  i  3  13  x 


2  +  (h.  .  -  Z.  +  Z .  )2  -  (h..  -  Z. 


13 


13 


y/d.  2  +  (h.  .  +  Z.  +  Z  J2  -  (h.  .  +  i. 
13  13  1  j  13  1 

v4.  2  +  (h. .  +  Z.  +  Z. )2  +  (h..  +  Z. 
13  13  1  3  13  1 


v4.  ,2  +  (h. .  -  Z+  2Z.)2  +  (h  -  Z . 

13  i3  i  j  ij  1 

y/d,  2  +  (h.  .  -  Z.  +2 Z.)2  -  (h.  .  -  Z. 

13  i3  i  j  ij  1 

y/d.  *  +  ( h .  .  +  A .  +  2Z  . ) 2  -  (h.  .  +  r 

13  13  1  3  13  i 


’V ♦  <hi3  ♦  \  *  2v  +  ihi3  *  *i 
V  *  hi32  -  V 

y/d.  2  +  h.  i  +  h.  .] 

1.3  13  13 

A*  *  lhij  *  v2  -  lhi3  *  v 


•V^Vv  *  ,hi3  “j*1 


y/d±2  +  (hij  +  2  Aj)2  '  (hij  +  2  £j)] 


V&tj2  +  (hij  +  2  £j)2  +  (hij  +  2  +j)3 


+  *.)] 

+  ij] 

+  ij)] 

+  ij)] 

+  2i j  )  ] 

+  2i  )] 

+  2i  )] 

+  2  v] 


d^j  is  the  perpendicular  distance  between  the  two  dipoles ,  and  h^  j  is  the 
distance  measured  from  the  feed  point  of  the  i1"*1  dipole  to  the  nearest  end  of 
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The  mutual  impedances  between  the  itk  and  jth  real  elements  is  calculated 

ija 

from  Equation  A-354  with: 


jrtl 

d  =  (  l  d  )  cos  (a  +  a  ) 

ij  rt=i  n  2  3 


( A-357 ) 


j?1 

h  =  (  l  d  )  sin  (a  +  a  )  -  l 

ij  n=i  n  2  3  j 


( A-358) 


R>r  the  typical 
when  i  *  j  with 
i  *  j,  Equation 


element  Z^N  +  jja  in  Equation  A-334,  Equation  A-354  is  used 
the  same  d^  as  Equation  A-357,  and  h—  =  Hj  +  •  Mien 

A-355  is  used  with  h^j  =  2H^  - 


The  feed-point  currents  are  determined  by  Equation  A-318  if  the  free-space 
approximation  is  applied. 


The  input  impedance  for  the  free-space  approximation  is  found  by: 


Z  =^[Z  +  R  '  Z  ]  I  =  R  +  j  X  (A-359) 

in  i=1  11a  v  1  (N+i)a  ia  in  in 


where 


R  " 
v 


k2  +  k 


(R  with  9 
v 


0) 


(  A-360) 


Fbr  the  perfect- ground  approximation.  Equation  A-333  is  used. 

The  input  impedance  for  the  per feet- ground  approximation  is  found  by  : 


Z"»  ?  (Z  +R'Z  ]  I  '  «  *  ”  +  j  X  "  (A-361 ) 

in  i*1  11a  v  1 (N+i)a  ia  in  in 
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The  expression  for  the  directive  gain  for  the  vertically  polarized  log- 
periodic  dipole  array  becomes: 


gd  = 


a 

120  12 


R. 

in 


+  B 


12 


sin2  9 


(  A-362 ) 


where  R^  is  given  by  Equation  A-359  for  the  free- space  approximation,  and  by 
Equation  A-361  for  the  perfect-ground  approximation.  A12  +  3Bt2  is  9iven  by 
Equation  A-351. 


The  radiation  vector  is  given  by : 


N 

9 


2  I 

bi 


i=1  sin  ki 

i 


cos  (ki  cosQ )  -  cos  ki  jky  sec(a 
i  i  i  2 

-  e 

sin  6 


+  a  )  cosiji^ 
3 

( A-363 ) 


The  free-space  field  intensity  is  given  by: 


E 

9f 


-  jkr 

1 

j60  e 


r 

1 


I  cos  (ki  cos6 )  -  cos  ki 

bi  i  i 


sin  ki  sin9 

i 


jky  sec  (a  +  a  )  costfi' 
i  2  3 

x  e 


The  attenuation  relative  to  free-space  is  : 
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ADB 


20  log 


I  A  | 


The  electric  field  intensity  in  the  diffraction  region  is  found  by: 


-  jkd 

Eq  =,  j30k  •_  (Mq)  (2  Fr) 


0 


CURTAIN  ARRAY 


Btor  a  curtain  array: 


Z 

Z 


1 

i 


zi"zi-1 


Yi 

VYi-1 


N 


M 


X1 

J l 


=  height  of  the  lowest  element  above  ground 

=  height  of  the  i  element  on  each  bay 

=  vertical  spacing  of  elements 
=  horizontal  position  of  the  (i  +  n)th  bay,  YQ 
=  horizontal  spacing  of  bays  (bay  separation) 

=  number  of  bays 

=  number  of  elements  in  each  bay 

=  reflector  spacing 
=  half-length  of  element  - 


0 


Ihe  field  components  produced  by  the  first  bay  are: 


'81 


■  j60  I 


- 

e  _  cos  (k&  cosfr)  -  cos  k& 

m  r  .2. 

sin  ifi 


siraj)  (A  +  j  B1 3  ) 


( A-366  ) 


(A-367  ) 


(A-368) 


(A-369) 
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•m  n 

_  .  _  e  J  cos  (k l  cosii/)  -  cos  k£ 

E  .  -  ^60  I  -  - i - 

4>  1  m  r  .2. 

sin  ip 


cos<t>  (A..  +  j B  )  (A-370  ) 

14  14 


where  r  is  the  distance  between  (Y ,  Z)  =  (0,  0)  and  the  far- field  point. 


jkz  (cos0  -  cos6  ) 

M  i  Or 

A  + j B  =  E  C  e  cos9  (1-R  e 

13  13  i=1  i  l  v 


j2kZ  cos0 


( A-371 ) 


A  +  j  B 
14  14 


/ 2  2  /2  2  - j2kZ  cos0 

vn  -  sin  9  2  vn  sin  9  i 

+  (1-R  )  F  -  (sin  9  -  -  cos0)  e 

v  e  2  2 

n  n 

jkZ  (cos9-cos0  }  -j2kZ  cos0 

M  i  Of  i 

E  C  e  1+R  e 

i=1  i  l  h 


(A-372) 


+  (1-R  )  F  e 
h  m 


■j2kZ  cos9 
1  1 


cosi|>  =  sinQ  sin$ 


represents  the  relative  amplitude  excitation  of  the  ith  dipole  with  C, 
and  0Q  is  the  desired  position  of  the  beam  maximum. 


The  array  factor  for  N  bays  is : 


jky  sinQ  (sin$  -  sin9  ) 
N  n-1  o 

S  =  E  e 
y  n-1 


( A-373) 


where  determines  the  position  of  the  beam  maximum  in  an  azimuthal  surface 
assuming  that  the  amplitude  excitations  for  corresponding  elements  in  each  bay 
are  the  same. 
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The  array  factor  resulting  from  the  arrays  of  the  real  elements  and  their 
images  with  respect  to  the  screen  at  x  =  -x1  is: 

-j2kX  cos<()  -jkX  cosiJj 

S  =  1-e  X  =  2j  e  X  sin  (kX,  cosiJ>  )  (A-374) 

X  l  x 

where  cos,px  =  sin  9  cos  <j>  . 

The  total  field  components  become : 


E 


9,t 


E 


91 


S 

x 


S 

y 


(  A-375) 


>1 


s  s 


( A-376) 


The  input  resistance  of  individual  elements  can  be  expressed  as  : 


Re  (Z±  ),  i  =  1,2,...,  M;  n  =  1,2, ...,N 


( A-377  ) 


(here  Z^n  representing  the  input  impedance  looking  into  the  ith  element 
of  the  nth  bay  can  be  calculated  by  considering  the  following  four  groups  : 

first  group  -  open-circuit  self  or  mutual  impedances  among  the  real  elements 

nN 

themselves.  R>r  example,  z.  is  the  mutual  impedance  between 

1 M 

the  i*"*1  element  of  the  nfcl1  bay  and  the  Mfcl1  element  of  the  Nth 
bay,  and  Z^”  is  the  self- impedance  of  the  itfl  element  of  the 
nth  bay. 
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Second  group  - 


Third  group  - 


fourth  group  - 


open-circuit  mutual  impedances  between  the  real  elements  and 

nN 

their  imperfect  images  under  the  ground.  for  example  Z.  is 
-  1  Mp 

the  mutual  impedance  between  the  ith  element  of  the  nth  bay 
and  the  imperfect  image  of  the  M'"n  element  of  the  Ntn  bay. 

open-circuit  mutual  impedances  between  the  real  elements  and 
their  perfect  images  with  respect  to  the  vertical  conducting 
screen . 

Open-circuit  mutual  impedances  between  the  real  elements  and 
the  imperfect  images  of  the  perfect  images  behind  the  vertical 
conducting  screen. 
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LThird  Group 
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nN  nN 
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n1  +  znl  + 

ilt  i2t 


nN 
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(A-378) 


The  total  power  supplied  to  the  entire  system  when  all  the  elements  are  of 
equal  length  21  is : 
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W, 

in 


+  c 


R 

m 


R!2  + 


+ 


+ 


( A-379  ) 


N  2  N  2  N  2 

+  R  +  C  R  +  +  C  R}  sin  kit 

12  2  mm 


2  2  H  M  2  n 

=  |  I  |  sin  kJ l  l  I  C  R  ( C  =  1 ) 

m  n=1  i=1  i  i  1 


* 

C.  =  1  unless  otherwise  indicated. 

The  directive  gain  for  the  curtain  array  is  : 


480  |S  |2sin2(kX  cos^  )  [sin2<j>(A  2  +  B  2  )+cos2<MA  2+  B  2)] 
y  ^  _  x  i  j  i  j  1 4  i  q 

sin2  ki  ?  ?  C.2  R.n 

n=1  i=1  i  i 


(A-380) 


[cos  (kX.  cos<(>)  -  cos  kX]' 

.  4, 
sm  ij; 


the  fre^-space  field  intensities  are : 


j60  I  e 
m 


•jkr 


cos  (k X  cos^)  -  cos  kX 


9f 


sin<j>  cos0 


jin  >t< 


(A-381 ) 


X  ? 


i=1 


jkZ  (cos0  -  cos9  ) 
i  o 

C  e  ]  S  S 

i  x  y 
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j60  I  e 
m 

E  *  - 

tf  r 


cos  (k£  cost)  -  cos  k i 
2 

sin  t 


cos  <(> 


.  ? 


i=  1 


jkZ  (cos6  -  cos9  ) 
i  o 

C  e  ]  S  S 

i  x  y 


The  attenuation  relative  to  free  space  is  : 


A  = 


'iy2  - '  v2 
>y 2  * 1  v2 


Hie  electric  field  intensities  in  the  diffraction  region  are: 
-jkd 

j60  1  e 

m  cos  (k£  cost)  -  cos  k£ 

E  =  -  -  sint  cos0 

0  d  2 

sin  t 


jkZ  (cos0 
M  1 

x  I  C  e 
i 

i=1 


COS0  ) 
o 


]  S  S  2F 
x  y  r 


-jkd 

j60  I  e 

m  cos  <k£  cos^)  ”  cos  kl 

d  2 

sin  t 


cost 


jkZ  (cos0  -  cos0  ) 

Mi  o 

x  l  C  e  ]  S  S  2F 

i  x  y  r 

i=1 
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( A-382  ) 


(A-383  ) 


( A-384  ) 


(A-385) 
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SLOPING  DOUBLE  RHOMBOID 

R>r  the  sloping  double  rhomboid: 

a  =  angle  between  the  principal  antenna  axis  and  the  major  rhomboid 
axis ,  in  degrees 

6  =  angle  between  the  major  rhomboid  axis  and  the  long  leg,  in  degrees 
Y  =  angle  between  the  major  rhomboid  axis  and  the  short  leg,  in  degrees 
*  short  leg  length 
J&2  =  long  leg  length 

H1  -  feed  height 
H4  =  termination  height 

a'  ~  angle  subtended  between  the  wires  and  the  horizontal  plane,  in 
degrees 


L  =  l2  cos  (0  +  a)  +  Jlj  cos  (a  +  t) 


cota^ 


cos  (a  +  y)  +  l 2  cos  (a  +  S) 


( A-386 ) 


(A-387) 


R>r  wire  1 : 

9'  =  9'  -  a ",  <)>'  =  (ff/2)  +  a  +  y 

cos^  =*  cosG  cos  ((ir/2)  -  a' )  +  sin9  sin  ((tt/2)  -  a')  cos  ( 4>  —  (tt  /2  )-a-Y ) 
=*  cosQ  sina'  +  sin0  cosa'  sin  ($-a~Y) 

R>r  wire  2 : 

9'  -  (ir/2)  -  a',  =  (tt/2)  +  a  +  0 

cos0  sina'  +  sin9  cosa'  sin  (<fc-a-0) 


cosi|>. 


(A-389) 
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Ebr  wire  3 : 


0'  - 
cosij^  = 

(tt/2)  -  a.' ,  <t>'  =  (it/2)  - 

cos0  sina'  +  sin0  cosa' 

(a+B) 

sin  (<ji+a+B) 

( A-390 ) 

Ebr  wire  4 

: 

0'  = 

cos<(/,  = 

4 

(IT/2)  -  o'  ,  <fi'  =  (TT/2)  - 

cos0  sina'  +  sin0  cosa' 

(a+Y) 

sin  ( 4>  +  a  +  Y ) 

( A-391 ) 

The  height 

of  the  current  element 

above  ground  for  wires  1,  2,  3, 

and  4  is : 

H 

s 

H1  +  s  sin  a" 

(A-392) 

Ebr  wires 

6  and  7,  it  is : 

H 

s 

H2  +  s  sin  a' 

( A-393 ) 

Ebr  wires 

5  and  8,  it  is  : 

H 

s 

+  s  sin 

(A-394) 

Hie  current  distributions  on 

1 

the  wires  are  given  by: 

I^s)  * 

-I  e~jks 
m 

( A-395 ) 

I3(s)  = 

I  e"jks 
m 

(A-396) 

I8(S)  * 

,  -j*,  -jM2 

I  e  e 

m 

► 

First  Rhomboid 

(A-397) 

1,(3)  - 

-jki, 

_  -;jks  1 

-I  e  e 

( A-398 ) 
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X2(S) 

= 

-X  e'jkS 
m 

1.(8) 

4 

= 

I  e'jkS 
m 

I?(S) 

= 

X  e‘jkS  e"jkA1 
m 

► 

I5(s) 

= 

-I  e"jkS  e"jk*2 
zn 

where 

* 

H2 

= 

H .  +  l  cos  (a+Y) 

tan  a' 

H3 

= 

H ^  cos  (<*+T) 

tan  a' 
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The  field  components  are: 


-jkr 


E.  =■  30  I  - 

01  m  r 


cosa'  cos0  sin  (4>-a-Y) 


■  j2kH 


a  “  j  ki.O/ 

COS0  „  1  1 


1-e 


R  e 
v 


1- 


U. 


( A-399 ) 


( A-400 ) 

Second  Rhomboid  (A-401) 


(A-402) 


(A-403) 


(A-404) 


(A-405) 


- j2kH  cos0  /2  .  2.  _  ri  ~ 27* 

.  , .  _  ,  „  1  vn  -sin  0  ,  ,  2„  vn  -sin  0  1-e 

+  ( 1  -R  )  F  e  — -  (sin  0  -  v - - -  cos0)  - — 

v  e  2  .  2  U„ 

n  cos0  n  1 


-Jkt, 


sinct'  sin0 


1-e 


a. 


+  R  e 
v 


-j2kH  cos0  ^ki,1U1 

1  1-e 
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+  (1-R  )  F 

v  e 


-  j2kH.(  cos0 


-jksur 

1-6  . -  (sin20  - 


^/n2,  -  sin20 


COS0  ) 


.11 
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-  jkr 


Ea,  =  -30  I 


'03 


m  r 


cosa'  cos0  sin  ($+a+0) 


1-e 


•)M2o3 


u. 


R  e 
v 


- j2kH  co£Q  ”^ki2U3 

i  I  -e 


(A-406 ) 


U. 


+  (1-R  )  F  e 

v  e 


-j2kH  cos0  f~l  '  2. 
1  Vn  -sin  0 
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n  cos0 


( sin  0 


n — ~~t 

✓n  -sin  0  1-e 


■jk*2vi 


COS0  ) 
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+  (1-R  )  F  e 
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-j2kHlCoS0  1_e"j,a2V 


V 


(sin20  - 


y/^-sin"© 


j  _2, 


COS0  ) 


'08 


-jkr  -jk£  U 
30  I  - -  e 


m  r 


-cosa'cos0sin(4>-a-Y ) 


-jk£  U  ' 

1-e  1  1  -j2kH  cos0 

— — =-s -  Re  J 

U-  v 


1-e 


( A-407 ) 


1 


-j2kH,cos0  .  -**lV  J~2 - —a 

*  (1-R  )  F  e  3  ^ - 5 -  vn_-sin_0  2 

v  e  u_  2  v  s 


'1 


n  cos0 


sin  0 


COS0  ) 


+  sina'  sin© 


r  “jk£,U,  -  jk£ „  U  ' 

1  1  •j.g  1  1  -j2kH^cos0 
tT!  +  u7  Rve 


1-e 


1 


1 


- j2kH  cos0  -jkA1Ur 

+  <1-R,)  F  e  1-6  .  - -  { sin20  - 


sin  0 


cos© ) 
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e~jkr  -jkZ1U1  J 

30  I  — -  e  cosa'  cosQ  sin  (A+a+0)  — - — 

m  r  T  U 

'jk£2U3  -j2kH  cos0 

1-e  J  2 

- 57 -  Rv  •  - 


i*V3 


(A-408) 


- j2kH  cos@  “jk£2U3  f~2  ~27 

x  Fg  e  - p -  v— -  (sin  0 

3  n  cos0 


/  2  2  J 
/n  -sin  0  „  1-e 

- - -  cos6)  -  sma  sm0  - 

n  J  3 


-Jk<2°3 


•jk£2U3  - j2kH_cos0  - j2kH  cos0  ,  _:ik£2U3 

+  -  Re  2  +  ( 1  -R  )  F  e  2  — - : - 

U3  V  V  e  U3 

,  j  2„  \m  -  sin20 
x  (sin  0  -  - - -  cos0 ) 

n 


=  -30  I  -  I  cosa'  cos0  sin  (<j>+a+y) 


-  j2kH.)cos0 

x  R  e  +  (1-R  )  F 

v  v  e 


-jk£1°4 


-jk£iV 


(A-409) 


- j2kH  cos0  ”^k£1U4  /l  .  2.  _ 

x  e  1  -  Vn_Z|in_0  (gin2 

4  n  cos0 


/t— 27  r.  -jk£i°4  -^k£iV 

- - cos0)  +  sina'  sin0  — -  +  — - - 

n  U4  °4 


-j2kH.cos0  -jZkH.cosO  1  ^k£1U4 

x  Rve  +  (1-Ry)  Fe  e  — *  „  7 -  (sin20 

4 


/  2  "  2  n 

fn _ -  sin  0 


CO80  ) 
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'87 


-Jkr  -JM  D 

30  I  -  e 

m  r 


-cosa'cosSsinf^-a-B) 


1-e 


-jM2o2 


1-e 


-*W2V 


(A-410) 


U„ 


-j2kH  cos6 

Re  +  (1-R  )  F 

v  v  e 


x  e 


•  j2kH_ cos6  ,  ~jia2U2  f~2 
2  i  -e 


[2  .  2. 

Vn  -  sm  9  ,  .  2„ 

— - -  {  sin  0 


n  cosQ 


[2  727 

vn  -sxn  9 


COS0) 


+  sina"  sin0 


-JM2«2  -jm2d2- 

1-e  _  1-e 


U_ 


U. 


-j2kH  cos0  -j2kH  cos0  ^k^2U2 

xRe  +  ( 1  -R  )  F  e  — -  (  3in  0 

v  v  e  U_ 


v42  ■  2 


sin  0 


COS0  ) 


'05 


-jkr  -jkZ  U 

30  I  - -  e 

m  r 


cosa'  cos0  sin  (<j>+a+Y) 


1-e 


"^kS°4 


0. 


1-e 


-^kAiV 


( A-411 ) 


U. 


-j2kH  cos0 

Re  +  (1-R  )  F 

v  v  e 


- j2kH,cos0  .  ~iki^4  /*2  2 

3  1-e _  Vn  -sin 


x  e 


-  sina"  sin0 


9  (sin2  0 


y42-sin20 


COS0  ) 


4  n  cos0  n 

-jk£.U.  -jki.U  '  „ 

1-e  1  4  1-e  1  4  -j2kH3cos0 

-  +  - — -  x  Rv  e 


- j2kH  cos0  "jkA2U4 

+  <1-RV)  e  — — -  (sin  9  - 


/I  7 ~a 

Vn  -sin  9 


cos© ) 
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!02  =  30  rm  r  |cosc^  cosQ  sin  (<t>-a-B) 


-)M2o2 


-  j2kH.)cos0 

Re  +  (1-R  )  F 

v  v  e 


( A-412 ) 


-  j2kH  cos0  12  ~2 ~  !~ 2 - 5  "jk 

■  '  'tS-SHlS  (sin2e  .  4j3J»!i  cose) 

n  cosO  n2  U2 


-jfc«2o2- 


-  sina'  sin0 


-Jkt2U2  , 

T-e  1 -e 


x  R  e 
v 


- jZkHj cos0 


-j2kH  cos0  ~jia2V 
i  1-e 


+  ^-RJ  *  *  -  (sin20 


/n _ -  sin  0 


cos0 ) 


0  j 

“30  Im  — ~ —  cosa'  cos  (4>-a-y )  — - 


-3Ml0, 


-j2kH1cos0 

x  R  e  +  (1-r  )  F  e 

n  h  m 


-j2kH,cos0  ~jk£1Ui' 
i  1-e 


(A-413) 


cosa'  cos  ($+a+B) 


"jkA2Ui 

1-e  2  3  .  i-e  2  3 


- j2kH  cos0  - j2kH  cos0  „  “jia2V 

x  *h  e  +  < 1  -Rh)  V  - 

3 


(A-414) 


-jkr  -jMD 

30  i  2 -  e  2  3 

m  r 


'  -jk^O,  ~jkA1u1' 
cosa'  cos  (4-a-Y)  -^2® - -  +  ll® _ _ 

U1  V 


-  j2kH  cos0  -  j2kH  cos0  ,  -jkAiUi‘ 

x  Ry,  •  +  (1-R.  )  F  e  3  — - 

“  h  m  u  ' 

1 


(A-415 ) 
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>6 


e-jkr  -jk^1U1 

-30  I_  — - —  e  cosa'  cos  (<|>+a+B) 


m  r 


1-e 


•JM2u3 


(A-416) 


1-e 


u„ 


-j2kH  cos0  -j2kH  cos0  “jkA2U3 

e  +  (1 -Ru )  F_  e  ‘ 


h  m 


U. 


-jkr 


E  =  -30  I  -  cosa'  cos  ( 4>— ot— 6  ) 


r 

1-e 


♦2 


m  r 


jk^2U2  -j2kH  cos0 
1-e  1 

+ - U7 - V 


u„ 


(A-417) 


+  (1-R.  )  F  e 
n  m 


-j2kH  cos0  ■’k^2U2 
1  1-e 


U 


-jkr 


E  =  30  I  -  cosa"  (<|>+a+Y) 


>4 


m  r 


1-e 


^k*1U4 


j  lei  1 U4  -j2kH  cos0 
♦  -  Rh  e  1 


(A-418 ) 


+  (1-R.)  F  e 
n  m 


-j2kH  cos0  ~3k£1°4 
J  1  1-e 


U. 


>7 


e-jkr  -jk£1U4 

-  30  l_  — - —  e  cosa'  cos(i}(-a-0) 


m  r 


1-e 


•jkV2 


U„ 


~jk*2U2  -j2kH  cos0 

^  1-e _  _  J  2 

U'  Rh  ® 


+  (1-R.)  F  e 


■j2kH2cos0  1 


_e-jkV2'] 

rt  I 


( A-419 ) 
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>5 


e-jkr  -jk*2U2 

-  30  I  -  e  cosa'  cos  (di+a+y) 

in  r 


1-e 


1-e 


-jkSU4 


U. 


V 


Rh  6 


-j2kH3cos6 


(A-420 ) 


+  (1-R.)  F 
n  m 

- j2kH^cos0 

1-e 

G 

V 

COSl|<1  ■* 

=  -cos0  sina' 

+  sin0  cosa' 

sin  ( 4>  — at— Y ) 

(A-421 ) 

cos  4<2 

-  -cosQ  sina' 

+  sinQ  cosa" 

sin 

(A-422 ) 

COSt3' 

-  -cos0  sina' 

+  sin0  cosa' 

sin  (iji+a+0) 

( A-423 ) 

cosi|»4' 

-  -cos©  sina' 

+  sin0  cosa" 

sin  (<j>+a+y) 

(A-424 ) 

U±  -  1-cosij^,  i  -  1,2, 3,4 


(A-425 ) 


U±  =  1-cosi^',  i-1,  2,  3,  4 


(A-426 ) 


The  total  fields  are 


8 

Z  E 
i-1 


0i 


(A-427) 


8 

Z  E 
i-1 


♦  1 


(A-428) 
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The  expression  for  directive  gain  is : 


gd  = 


r2  (|  Eg  |2  +  lEj2) 


120  I  R. 
m  in 


Since  the  two  rhombics  are  connected  in  parallel,  the  total  input 
2  3^  e~^ks,  and  Rin  may  be  taken  as  =  300&  . 

Hie  radiation  vectors  are  given  by : 


N 


01 


02 


03 


04 


N 


05 


N 


06 


07 


I 

m 

-j*V, 

1-e 

jk 

I 

m 

-jkt2<j2 

1-e 

jk 

°2 

I 

m 

-Jkt2o3 

1-e 

jk 

°3 

X 

m 

1-e 

jk 

°4 

X 

m 

-jkt,u4 

1-e 

jk 

°4 

I 

m 

1-e 

jk 

°3 

X 

m 

-i«  U 

1-e 

jk 

U2 

[cosa'  cos©  sin  (<f>-a-Y)  -  sina'  sin0] 


[cosa'  cos©  sin  (4>-a-0)  -  sina'  sin0] 


[-cosa'  cos©  sin  (<£+a+0)  -  sina'  sin0] 


[-cosa'  cos0  sin  ($+<x+y)  -  sina'  sin0] 


i  k£  u 

2  2  [cosa'  cos©  sin  (ij>+a+Y) 
-  sina'  sin0] 


■^1U1  [cosa'  cos©  sin  (<{>+a+Y) 
-  sina'  sin0] 


^k*1U4  [-cosa'  cos0  sin  ($-a-0) 
+  sina'  sin0] 


( A-429 ) 


current  is 


(A-430 ) 


(A-431 ) 


(A-432) 


( A-433 ) 


(A-434 ) 


(A-435) 


( A-436 ) 
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The  total  radiation  vectors  are : 


8 

N  =  Z  N  ( A-446 ) 

0  i=  1  01 


8 

N  =  E  N  (A-447 ) 

<j>  i=  1  (j>  1 


Hie  free-space  field  intensities  are  given  by: 


'0f 


j30ke 


-jkr 


N- 


(A-448) 


V 


j30ke 


-jkr 


(A-449) 


The  attenuation  relative  to  free  space  is : 


A 


,/'Ee|2  -  'V2 
v  iv2  +  >e 


2 


( A-450  ) 


and 


ADB  -  20  log.,,  1 

10  TTT 

The  electric  field  intensities  in  the  diffraction  region  are: 


(A-451 ) 


'0 


$ 


j30ke 


-jkd 


j30ke 


■jkd 


(N„)  (2  F  ) 
o  r 

(N.  )  (2  F  ) 
4>  r 


( A-452  ) 


( A-453 ) 
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APPENDIX  B 

INTEGRALS  ENCOUNTERED  IN  VERTICAL-MONOPOLE  CALCULATIONS 


In  this  appendix,  the  following  integral,  encountered  in  vertical- 
monopole  calculations,  is  evaluated: 


l  -jkZ  cos9 

J  =  /  e  sin  k(l-Z)  j  y/rp~  exp  (-P  )  erfc  ( j  y/P~ )  dz  ( B-1 ) 

3  0  e  e  e 


From  King  (see  Reference  27,  pp.  757-758): 


— — -  Jtt P  exp  ( -P  )  erfc  ( j  J~P~)  -  k  'Ax 

z  *  e  e  *  e 


2  -  sin2e  /Q“  e  -jkg  dZ'  ( B-2 ) 


where  z'  is  a  variable  introduced  by  Van  der  Pol  (see  Reference  17) 


g  *  (n2  cos0  +  V^n2-sin29  )  Z"  +  ^ —  sin20 

r  r  2r  r 


(B-3) 


From  Equations  B-1  and  B-2: 

[■ 2  T 

£  -jkZcos8  2kV/n  -sin  0  00  -jkg 

J  *  /  sin  k(Z-Z)  e  j - - -  /  e  dZ'  dZ  (  B-4) 

3  0  1  -R  0 
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where 


h  h  +  z 

—  and  cosd  =  - 

r,  r  r 


2k  >/n2-sin29  rl 

1-R  ^0  '0  ® 


-J"  (l  *  4) 


sin  k(Ji-Z) 


x  exp 


“jk  n2h 

(Z'  —  + 
r 


Z'  Vn  -  sin  0 


nV  2  , 

sin  9 )  dZ  dZ 


(  B-5 ) 


Changing  the  order  of  integration  : 


/  2  2  2  2  2 
2k yn  -sin  9  »  l  -jkbZ  (a  Z'  sin  9+Bz') 

)  -  /  /  e  sin  k(£-Z)  dZ  e  dz' 

1-R  0  0 


(B-6) 


where 


b  «  h  +  nV 
r  r 


(  B-7  ) 


4 

-jk  T- 
2r 


(  B-8 ) 


■jk  +  i/n2  -  sin29) 
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l  -jkbZ 

/  e  sin  k(A-Z)  dZ 

0 


1 


2 

k  (1-b  ) 


-  jkbi. 

e  +  jb  sin  k£  -  cos  ki 


Neglecting  b^  since 


bi 2  <  < 


gives : 


J l  -jkbZ 

/  e  sin  k(i-Z)  dZ 

0 


1 

k 


-jkbi 

e  +  jb  sin  k l 


cos  k£ 


Substituting  Equation  B-11  into  &-6  gives: 


J 

3 


2  2 
2k\/n  -sin  0 


1-R 


®  1  -jkbi 

J  —  e  +  jb  sin  k£  -  cos  k l 

Ok. 


2  2  2 

exp  (a  z'  sin  0  +  0Z" )  dZ' 


J 


3 


i  2  vn-sin^e  e-jkicos0  ^  +  (jsin  k£  cos0  _  cos  k£ )  j 

J  1  -R  31  J  ^ 

V 

2 

+  3  “  sin  k l  J33 


where 


A 

dz"  =  j -  exp(-w  )  erfc(  j 

2asin0  1 


“  2  2  2 
J  =  /  exp  a  Z  sin  0+0  Z' 

31  0  1 


(B-10) 


(B-11) 


(B-12) 


(B-13) 


( B-14) 
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J32  *  Q  **P  2'2  sin2B  +  62']  32'  -  "]  2os in6  e"”  er<°  <3\/“>  ( B-15) 


J33  -  r0  2'  exp 


p  |a2Z'2  sin29  +  Bz'j  dz" 

- — —  |l-j  \/  erfc  (j  /“)] 


2a  sin  9 


( B-16 ) 


„  3,  „  n  n  (h+i)  .  /  2  .2. 

p  -  Dkx.  —  *  -3k  -  +  y/n  -  sin  9 

r  r 


(  B-17) 


.  2  .2. 

4a  sin  9 


jkr  h  +  &  /n2  -  sin29  2 

2  sin2 9  r  n2 


(  B-18) 


2  2 
4a  sin  9 


2  sin29 


/!  .  2~ 
h  /n  -sin  9 

r  2 

n 


(B-19) 
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Substituting  Equations  B-14,  B-15,  and  B-16  into  B-13  yields: 


J.  = 


n — rs 

vn  -si 


sin  9 


1-R 


- jkicosQ 


2 

n  sin6 


e  1  erfc  ( j  )  +  (jsin  k£cos9  -  cos  ki) 


1  /j2irr  -w  n2sink£  i2r  I —  -w 

~ - \  \ —  e  erfc(  j^)  +  - - - - -  {  1-3vAw  e  erfc^/w)) 

n  sin 0v  n  k  sin  9 

2  .  2.  -jk£cos0  r— -  -w_ 

vn  -sin  9  I  e  h2i rr  1  ,  , — ,  „  „  i 

■«  /  ^ —  e  erf  c(  3^7^)  +  (jsin  ki  cos0-cos  ki)^— 


n  (1-RV) 


sin0 


Vj2irr  -w  _  ,  .  r~.  .  j2sin  ki  r.  .  / —  -w  _ 

d  k  e  erfc(3yw)  +  J - —  {l-Uy/rrw  e  erfc  (3^/w)} 

k  sin  0 


( B-20 ) 


Since  W1  =  W  =  Pg: 


J-  = 


(  B-21 ) 


j2  F  sinki 
e  . 

2 

k  sin  0 
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APPENDIX  C 

INTEGRALS  ENCOUNTERED  IN  VERTICAL-DIPOLE  CALCULATIONS 


In  this  appendix,  the  following  two  integrals  encountered  in  elevated- 
vertical-dipole  calculations  are  evaluated : 


J 

3 


Z 

o 

/  -jkZ  cos0  sin  k (i-Z  +Z) 

Z  -l  e  o 

o 


L 


* 

lift" ) 

e 


(  C-1  ) 
dZ 


J  " 
e 


Z  +£. 
o 


- jkZcosS 


sin  k  U+Z  -Z) 
o 


-P 


1-j  y-rrP  e  erfc  (  j  ^5” ) 


dZ  (  C-2 ) 


The  first  parts  of  the  integrals  are  straightforward.  The  second  part  of  J3 
is  : 


Z  -P 

o  -jkZcosS  _  e 

J  =  /  e  sin  k(£-Z  +Z)  j  v/ttp”  e  erfc  ( j  /p")  dZ  (C-3) 

31  Z  -i  o  e  e 

o 


Allowing  the  same  techniques  as  those  used  in  APPENDIX  B: 


,2k 


-sin  0  c” 


31 


1-R 


/I  /7°r»  e“JkbZ  sin  k  (4-Z_+Z)  dZ 
Jo 


exp  {o^Z'^sinZ0  +  0Z'}  dZ' 


(C-4) 
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where  b,  ct  ,  and  8  are  given  by  Equations  B-7,  B-8,  and  B-9,  respectively. 

-jkb(Z  -A )  -jkbZ 

Z  o  o 

o  - jkbZ  e  -e 

/  e  sin  k(A-Z  +Z)  dZ  =  - 

Z  -A  o  2 

o  k  ( 1-b  ) 


(coskA+jb  sinkA) 


(  C-5) 


- jkb( Z  -A )  -jkbZ 
o  o 

e  -e  ( coskA  +  jb  sinkA ) 


Substituting  Equation  0-5  into  0-4  yields : 


\fFZi 


J  -  j 
31  1-R 


2  -jkb(Z  -A)  -jkbZ 

sin  8  «*  o  o 

•  /  {e  -e  (coskA+jb  sinkA)} 

o 


J31  *  3  1-R 


2  2  2 

x  exp  {a  Z'  sin  9  +  0Z"}  dZ" 


_  /  2  .2.  -jkZ  cosQ  ... 

.2  v  n  -sin  8  J  o  }kAcos6 


'311 


-(coskA  +  jcos8  sinkA)  J  2  -  j  ~  sinkA  J313 


where 


(C-6) 


(  C-7  ) 


311 


exp  (a  Z"  sin  9  +  B2Z')dZ"  *  -j  2ag^'nQ  e  erfc  (j  (0-8) 
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J  " 
31 


Z  +i  -P 

o  -  jkZcosG  _  e 

/  e  sin  k  (A+Z  -Z)  j  e  erfc  (j^/T-)  dZ  (C-15) 

Z  o  e  e 

o 


Rjllowing  the  same  procedure  as  that  used  to  evaluate  Equation  C-3: 


„  /  2  .  2q  -jkZ  cos9  ...  . 

_  2  vn  -sin  9  o  -jk£cos9 


31 


1-R 


311 


( C-16 ) 


-  ( coski- j  cos6  sink£)  J312  +  3  “  sinki.J313 


where 


-  W 

-  2  2  2  /T"  4 

J  *  /  exp  (a  z"  sin  9  +  6  Z" )  »  -j  -  e  erfc  (jy/w”)  (C-17) 

311  o  4  2asin0  4 


84  *  B-jk  —  <  VZ)  *  "jk 


n  ( h+Z+2 )  ~ 

- 2 -  +.  /n2-sin20 


(C-18) 


4  4a2sin29 


-j 


kr 


2  sin  9 


htZ  U  /  2 

-r~*  V— 


sin29 


(019) 
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J  +  J 
31  31 


.  2  \Zn2-sin“0  '  o 


1-R 


2„  -jkZ^cos 0 
e 


-jkZ  cos©  _  ,  -jkicos8  „ 

o  J,.  +  e  J  J, . .  -  2  j _  coski 


311 


311 


312 


(C-20) 


Since: 


2 

and 
J 


WL 


311  311 


W4  =  P3 


312 


J  +  J 
31  31 


/l T  -jkZ  cos0 
2  *n  -sin  0  o  r 

j  - - -  e  J  2  cos  (kicos0)-2  coski 


1-R 


31 1 


J2  2_  -jkZ  cos0  -P 


2  Vn  -sin  0  o 

-  e 

1-R 


-  cos 


ki} 


e  _ 

-  e  erfc  ( j  /P  )  {cos(kicos0) 

asin8  e 

(C-21) 


J31  +  J3l' 


2  >/n2-»in29  'jkZoCOs6  .  feZ 
n2(1-R  )  S  V  * 


-P 


e  6  erfc  ( j  ) 


cost kicosO )-coski 
sin8 
( C-22 ) 
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-jkZ  cosQ 
2  e  ° 

J  +J  '  =  - - -  { cos (k£ cos0  ) -coski} 

3  k  sin  8 

-jkZ  cos6 

_  2e  ° 

1-R 

v 


-P 

x  e  erfc  (  j  JP) 
v  e 


cos(k£cos9)  -  cosk£ 
sin6 


Equation  C-23  can  be  simplified  further  to  yield  : 


- jkZ  cos0 

5-  ® 

J  +  J  "  *  - - -  {cos(  k£cos8  )-cosk£}  -  coskf} 

k  sin  0 

r 5  -p  i 

sin0  vn  -sin  0  rZ~,  e  .  .  rz~ . 

x  1  “  T^r”  - 2 -  VD2irkr  e  erfc  (Dx/Pg) 

v  n 


cos(k£cos9)  -  cosk£ 
sin^0 

2  -pe  ' 
sin  9  e  erfc  ( j  y P  ) 


(C-23) 


( C-24) 
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APPENDIX  D 

FORMULAS  FOR  CALCULATING  FTELD  STRENGTH  OVER  SPHERICAL  EARTH 
WITHIN  THE  RADIO  HORIZON 

The  formulas  presented  in  APPENDIX  A  are  for  calculating  field  strength 
and  gain  over  planar  earth  and  over  spherical  earth  in  the  diffraction  region 
beyond  the  radio  horizon.  The  formulas  presented  in  this  appendix  are  for 
calculating  field  strength  over  spherical  earth  within  the  radio  horizon. 

Within  the  radio  horizon,  the  effects  of  spherical  earth  are  taken  into 
account  by  use  of  the  divergence  factor.  Only  the  direct  and  ground-reflected 
contributions  to  the  fields  are  considered  in  this  formulation,  because  the 
surface  wave  is  negligible  when  the  feed  points  of  the  antennas  are  at  least 
several  wavelengths  above  the  surface  of  the  earth. 

The  quantities  in  the  following  equations  are  described  in  Section  5  and 
in  APPENDIX  A.  The  appropriate  formulas  for  calculating  directive  gain  from 
the  following  equations  are  given  in  TABLE  2  of  Section  5  which  also  provides 
expressions  for  the  radiation  efficiency. 

HORIZONTAL  DIPOLE 


ESD-TR-8Q-1 02 


Appendix  D 


VERTICAL  MONOPOLE 


e"jkRd  A2  +  jB2  e'jkRr  A2-jB2 

E0  Im  R  sin8  +  Rv  adiv  R  sin0 

d  r 


(D-3) 


E ,  =0 


(D-4) 


VERTICAL  MONOPOLE  WITH  RADIAL- WIRE  GROUND  SCREEN 


'jkRd  A2  +  jB  -jkRr  A  - j  B 

Ee  =*  330  Im  v- - sin9  +  Rv  adiv  S - X  (GSCE1  (D"5) 

d  r 


/2  2  2 

ka  -j  yx  +k  l  -jx 

n  sin0  /  e  -e  coski.  J  (xsinO)  dx 

0  J  1 

120it  sinki.  (cos  (k l  cos0)  -  coskiJ 


(D-6) 


E.  »  0 


(D-7) 
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Ee  =  j6°  Zm 


cos(k&  cosQ)  -  coski 
sin0 


“  jk^j 


"3*Rrl 


+  R  a . . 
v  div 


(D-8) 


0 


(D-9) 


INVERTED-L 


E 


9,v 


j30  I 


m 


A4+jB4  D 
sin6  +  v  °div 


sin9 


(D-10) 


(D-11) 


E 


9,h 


A5+jB 

-j30  I  - - -  cos9  sin* 

m  ,2 .  + 
sin 

n 


R 

v 


(D-12) 


A5+jB5 

E.  .  =  j30  I  - - -  cos* 

<M  m  sin  *  ' 


r  -3**, 


•jkR, 


Rh° 


div 


( D  —  1  3  ) 
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-j60  I 


e  cos(k&  cosd>)  -  cosk£ 


m  R . 


■  2, 

san  ip 


(cosa'  sind>  cos6  -  sina'  sin0 ) 

(D-14) 


-R  a..  — 

v  div  R 


cos(k&  cos^)  ~  cos  k l 
sin2^ 


(cos  a'  sin<(i  cos0  +  sin  a'  sin0) 


E,  =»  j60  I  cosa'  cos<}>  — ~ 

<f>  m  R . 


cos  ( ki  co sd)) _ -  cos  k£ 


“jfcRr 

e _  cos  (k&  cosd>' )  -  cos  ki. 


div  R 


sln2\p' 


(D-15) 


SLOPING  LONG-WIRE 


A,  +  j  B, 

E«  =  -13  0  I  — - -  - - —  (cos  a'  sind)  cos0  -  sina'  sin0) 

0  m  R  ,  2  . 

d  sxn  ip 

-  jkR 

J  r  A  +  jB 

-R  a. .  — - -  - - - -  (cos  a'  sin®  cos0  *  sin  a"  sinQ) 

v  div  R  .  2  . , 

r  sin  d» 


(D-16) 


E.  *  j30  I  cosa 


r-j> 

'  cosd) 


d  VjB6 


-T-  +  a,,.. 


V  +jB  1 

(D-17) 
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-jkR^  -jkit^ 

E0  =  30  I  cosa'  cos9  cos(4>-0')  - -  — - - R  a 

cn  m  T  r  n  v 


v  div  R 


sina'  sin6 


-jkRd  -jkXU 
e  1-e 


(D-18  ) 
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E._  =  -30  I  cosa'  sin(9+0') 

92  m 


-jkRd  -jkiU2 
e  1  -e 


-jkRr  -j^U4 

.  _  e  1-e 

+  R.  a  -  - 

n  div  R  u 


(D-21  ) 


E 


0 


'01 


'02 


(D-22) 


E .  = 


>1 


+  E 


92 


(D-23 ) 


TERMINATED  SLOPING  RHOMBIC 


The  total  field  components  are  given  by: 


E0  E01  +  E02  +  E03  +  E04 


and 


E  =»E  +E  +  E  +E 
9  91  92  93  94 


(D-24) 


(D-25 ) 


where  Eg^  Eg2,  E^,  and  E^2  are  given  by  Equations  D-18,  D-19,  D-20,  and 
D-21,  respectively,  and: 
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-jklu 


30  I  e 
m 


-cosa'  cos0  cos(<f>-£3') 


-jkR  -jkiu, 
a  7 

e  1-e 


-  R  a,.  e 
v  div 


-j2kH  cos0  "3kRr  "3kAu3 
1  e  1-e 


(D-26  ) 


+  sina'  sin0 


-jkR  -j  k&U. 

e _ f  1-e _ * 

Rd  U1 


+  R  a  e 
v  div 


•  j2kH  cos0  ~jkRr  ^MtJ3 

1  e _  1-e _ 

Rr  U3 


E  =  30  I  e 

04  m 


Jt  -^kRd  -i*lv2 

cosa'  cos0  cos((j>+0'’)  — -  — - - 

Rd  U2 


-  R  a  .  e 
v  div 


-j2kH  cos0  ~jkRr  .  _3kAU4 
i  e _  1-e 


(D-27 ) 


sina'  sin0 


-jkRd  -jkAtt 
e  1-e  " 


+  R  a  e 
v  div 


-j2kH  cos0  ~jkRr  -jk*U4 
I  e _  1-e 


E  =  -30  I  e 

<p3  m 


-jk£U. 


L  "jkRd  .  -jk£Ul 

cosa'  sin(<|>-0')  — - — - 

Rd  U1 


-j2kH  cos0  "jkRr  .  -3kXU3l 


(D-28) 
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E  =  30  I  e 

4>4  m 


-jktu,  -lkRd  ,  -jM02 

s  cosa'  sin(<f+0')  — -  - — - 

Rd  U2 


+  *h  “div  6 


TERMINATED  HORIZONTAL  RHOMBIC 


a  -jkR  - jkiU  " 

-l2kH.cos0  r  J  4 

J  1  e  1-e 


- j  kiU1  -jkRd  -jkRr 

E01  =  30  2»  -^V -  COS0  COS  (*“Y)  "  Rv  “div  h - 

1  dr 


(D-29  ) 


(D-30) 


-jk£U2  (“  -jkRd  -jkRr‘ 

E  =  .30  I  ^ -  cos0  C03(^Y)  ~~  ~  Rv  «div  h -  (D_31) 

2  dr 


-jUo  f  ”3kRa  e”3kRr 

E93  *  -30  1.  -  *  coa6  coslJ-Y)  -  R„  -  «>- 

1  dr 


■jkiU2  -jki.U  "jkRd 

1*6  J  j  e 

30  I  - - -  e  cos0  cos(<))+Y)  — - - 

m  U2  Rd 


R  a  — — 
v  div  R 


jkRr' 

R 


(D-33) 
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1-e 


sin  ( $-y ) 


e"jkRd 

£"!  +  ®h  “div 

a 


R 

r 


(D-34) 


_jkAU2 

V  =  "30  Zm  IZ15 -  Sin  (<t>+r) 


-  jkR ,  - jkR 

a  r 
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d  r 


(D-35) 
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+  R.  a . . 


div  R 


-jkR. 


(D-37 ) 


4 

E  =*  Z  E 
0  i=*1  0i 


(D-38 ) 


4 

E  -  Z  E  (D-39 ) 
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i 
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HORIZONTALLY  POLARIZED  LOG-PERIODIC  DIPOLE  ARRAY 
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